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ABSTRACT 



o 

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first 
time, this library provides the rate probability density at all temperature grid points 
for convenient implementation in models of stellar phenomena. The recommended rate 
and its associated uncertainties are also included. Currently, uncertainties are absent 
from all other rate libraries, and, although estimates have been attempted in previous 
evaluations and compilations, these are generally not based on rigorous statistical defi- 
nitions. A common standard for deriving uncertainties is clearly warranted. STARLIB 
represents a first step in addressing this deficiency by providing a tabular, up-to-date 



database that supplies not only the rate and its uncertainty but also its distribution. 
Because a majority of rates are lognormally distributed, this allows the construction of 
rate probability densities from the columns of STARLIB. This structure is based on a 
recently suggested Monte Carlo method to calculate reaction rates, where uncertainties 
are rigorously defined. In STARLIB, experimental rates are supplemented with: (i) 
theoretical TALYS rates for reactions for which no experimental input is available, and 
(ii) laboratory and theoretical weak rates. STARLIB includes all types of reactions of 
astrophysical interest to Z = 83, such as (p, 7), (p, a), (a,n), and corresponding reverse 
rates. Strong rates account for thermal target excitations. Here, we summarize our 
Monte Carlo formalism, introduce the library, compare methods of correcting rates for 
stellar environments, and discuss how to implement our library in Monte Carlo nucle- 
osynthesis studies. We also present a method for accessing STARLIB on the Internet 
and outline updated Monte Carlo-based rates. 



Corresponding author: anne.sallaska@nist.gov. Present address: National Institute of Standards and Tech- 
nology, Gaithersburg, MD 20899-8462. 



Introduction 



In the mid 20th century, Fowler and collaborators revolutionized stellar evolution simulations 
by emphasizing standardization of recommen ded reaction rates fo r ease of collaboration and com- 
parison in the stellar modeling community (jCaughlan fc Fowled (|1988l ) and references therein). 
The NACRE collaboration followed, providing refined nu clear physics data, and attempted to 
estimate reaction rate uncertainties never before addressed (jAngulo et al.lll999l ). However, the sta- 
tistical meaning of their "lower" and "upper" rate limits is unclear, and the need to establish a 
co nsistent defin i tion is evident. The turn of the century brought forth an additional evaluation 
by llliadis et al.l (J200ll ) , which introduced radioactive target nuclei and recommended standard res- 
onance str engths to use as normalizat ions. However, no rigorously determined uncertainties were 
presented. iDescouvemont et al.l ([20041 ) then presented an evaluation of reaction rates important to 
Big Bang nucleosynthesis and attempted to quantif y uncertainties in a more stati stically meaning- 
ful manner. Similar evaluations were published by (jAdelberger et al.lll998l . l201ll ) that focused on 
solar fusion cross sections and their uncertainties. 

Another paradigm shift is currently upon us. Although stellar reaction rate libraries embrace 
experimental rates over theoretical ones and cull their experimental data from the above evalua- 



tions, none of the l ibraries in c lude any estimation of uncertainties (JIN A REACLIB: ICvburt et al 



(J2010l ); BRUSLIB: IXu et al.l (J2013I )). This information is essential for realistic simulations of the 
energy generation and nucleosynthesis that occurs in various stellar phenomena. In the past, re- 
porting rates without uncertainties was accepted because the additional computing time to allow 
for uncertainties was exorbitant and because it was unclear how to proceed in a statistically rigor- 
ous way. This, however, can lead to erroneous conclusions drawn from the output of stellar models, 
as it is not uncommon for the uncertainties in many astrophysically-relevant reaction rates to span 
orders of magnitude. To produce realistic nucleosynthesis results and confidence levels, it is im- 
perative for modelers to be able to access a library that includes not only meaningful uncertainties 
but also probability density functions for each reaction rate at a given temperature. The advent of 
readily- available high-powered computing allows not only for applying a Monte Carlo technique to 
evaluate reaction rates but also for incorporating the reaction rate probability densities into stellar 
models. This paper presents a next-generation reaction rate library, STARLIB, built specifically 
to begin addressing these problems. 

STARLIB is a tabular, stellar reaction rate library that includes all types of incident particles, 
including neutrons, protons, a particles, and 7 rays. Target nuclei range between Z = 1 — 83. Ex- 
perimental data are used when available, and, if not, theory is considered. The str ucture of STAR- 



LIB rests on a new meth od to quantitatively define reaction rate uncertainties (JLongland et al. 



2010bl : llliadis et al. 



2010c 



bJJaD , although only a fraction of the library as yet contains quantitative 



uncertainties. This method uses experimentally determined nuclear physics quantities (resonance 
strengths and energies, S-factors, partial widths, etc.) as inputs to a Monte-Carlo algorithm. Each 
quantity and its associated uncertainty is sampled according to its physically-motivated probability 
density function. In 2010, experimental Monte Carlo rates for 62 charged-particle nuclear reactions 



on A = 14 — 40 target nuclei were evaluated (jlliadis et al.ll2010d1. The library pre sented here in- 
cludes updates to 7 Monte-Carlo rates and 1 entirely new rate. ISallaska et al. (20121) is an overview 
of preliminary work for STARLIB. In a recently published work (jlliadis et al.ll201ll ). a preliminary 
version of STARLIB was employed_j. 

This work is organized as follows. We begin with a discussion of the thermonuclear reaction 
rates in Sec.JH including analytical evaluation techniques (Sec. l2.ip . Monte Carlo methods (Sec. l2.2D . 
theoretical rates based on the TALYS code (Sec, 12.3] ). and modification of laboratory rates for stellar 
environments (Sec. 12.4]) . We present our prescription for building STARLIB in Sec. [3l and imple- 
mentation of STARLIB in Monte Carlo nucleosynthesis studies is discussed in Sec. |U Appendix lAl 
details the formatting of STARLIB, and Appendix [B] specifies how to access the STARLIB library 
and our Monte Carlo rate calculator on the Int ernet at Istarlib . physics . unc . edu[ Appendix O 
details updated Monte Carlo rates since 2010 (jlliadis et al.ll2010bl ). The complete library, with 
references, is provided on our website and on the Astrophysical Journal website. 



2. Charged-Particle Reaction Rates 

For a charged particle reaction in a stellar plasma at thermal equilibrium, the reaction rate, 
Na(ctv), for two interactin g particles with a Maxwell-Boltzm ann distribution of velocities at a 
temperature T is given by (jlliadis! 120071 ; JRolfs h Rodnevlll988l ): 



N A (av) 




/ Ea{E)e- E ' m ' dE, 
Jo 



(1) 



where /x is the reduced mass of the particles in the entrance channel, E is the center-of-mass energy 
of the reaction, and a is the reaction cross section. The rate may be simplified for resonant (non- 
interfering and i nterferi ng) and non resonant reaction mechanisms. Details of each are discussed in 
Longland et al.l (|2010bl ) and references therein. 



2.1. Analytical Methodology for Total Reaction Rates 

When applying the above formalism to compute the total thermonuclear reaction rate, two 
main questions arise: how does one account for uncertainties, and how does one incorporate ex- 
perimental or theoretical upper limits of partial widths and resonance strengths in a meaningful 
way? Only in special circumstances are analytical error propagation methods applicable for reac- 
tion rates, for example, when uncertainties in resonance energies are small (i.e., errors less than a 
few keV) and when upper limits on strengths or partial widths can be disregarded. In such cases, 



L In the preliminary 2011 version, the Hauser-Feshbach rates were adopted from iRauscher fc Thielemannl ( 20001 ) 



the uncertainties in both the strength and energy may be propagated using standard techniques 
of derivatives. However, in most cases there are many different contributions to the total reaction 
rate, and analytic methods quickly begin to fail. 

Fowler and collaborators ( Caughlan fc Fowlerlll98§) provided no un certainty informatio n what- 
soever. Although the NACRE collaboration (jAngulo et all LL999T) and llliadis et al.1 (|200ll ) did at- 
tempt to address this void, their "upper" and "lower" limits were not established from rigorous 
statistical methods. Often what is reported are limits where the authors have guessed the largest 
sources of uncertainty, and, although this undertaking is a laudable step in the right direction, 
it has no statistical merit. There is also the perpetually troublesome problem of how to handle 
upper limits of nuclear physics input quantities (resonance strengths, partial widths, etc.) in the 
calculation of the total reaction rate. The Monte Carlo formalism allows inclusion of upper limits 
in a natural, statistically precise way. 

The interpretation and the meaning of "upper" and "lower" limits can also be ambiguous. 
Because rate probability density functions are usually not reported, these limits have been treated 
as sharp boundaries in most nucleosynthesis simulations, if uncertainties were taken into account 
at all. Several stud i es att e mpted to quantif y the actual r e action rate probability density. While 
Krauss Sz Chaboyerl (120031 ): ICoc et al.l (12002 ); llzzard et al.l (120071 ) used unphysical, unif orm distri 



butipns, a few modelers proceeded a step further and assumed Gaussian distribution s (jCoc et al. 



20041 ) and even lognormal distributions for rates (IHix et al.l 120031 : iParikh et al.l 120081 ) . One of the 
main thrusts of the present work is to begin supplying not only rigorously determined values for 
recommended rates but also to provide reliable probability densities for a given rate at any tem- 
perature. 



2.2. Monte Carlo Rates 



To address the persistent probl e m of q uant ifying reaction rate uncertainties in a statisti- 
cally meaningful way. lLongland et al.l (J2010bl ) and llliadis et al.l (|2010d lbllal) detailed a Monte Carlo 
method founded on physically motivated distributions assigned to each nuclear physics input quan- 
tity. Here we present an overview of this method, which motivates the structure of the library 
presented in this work. 

To calculate the thermonuclear reaction rate, each input parameter must be described by a 
probability density function. The distributions employed in the present work are summarized in 
Table [TJ Once distributions are assigned for each input parameter, a random value can be sam- 
pled for each, and the reaction rate is computed from Eq. [T] (or its subsequent simplifications) at 
a given temperature using standard Monte Carlo sampling techniques. The sampling continues 
until the user's desired precision is reached, with a minimum of 5,000 samples in order to achieve 
reproducibility within a few per cent. All Monte C arlo calculations reported in this work are per- 
formed using the code RatesMC (llliadis et al.ll2010bl ). The method allows for a statistically rigorous 



Table 1: Summary of probability density functions used to estimate reaction rates. 



Parameter 



Distribution 



Resonance energies 
Resonance strengths 
Partial widths 
Nonresonant S-factors 
Upper limits (r, ory) 
Interference 



Gaussian 

lognormal 

lognormal 

lognormal 

Porter-Thomas 

binary 



definition of the low rate, median (recommended) rate, and high rate as the 16th, 50th, and 84th 
percentiles of the cumulative rate distribution, respectively. The values chosen here represent a 
coverage probability of 
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Fig. 1. — (color online) Three-dimensional reaction rate probability density function, normalized 
to the median (recommended) rate, for 22 Na(p,7) 23 Mg (see Appendix IC.6J) . Notice how the distri- 
bution becomes narrower with increasing temperature, i.e., the reaction rate uncertainty decreases 
because the nuclear physics input becomes more reliable at higher energies. 

Although the "upper limit" and "lower limit" values of rates presented in the literature are 
sometimes interpreted by nuclear astrophysicists as sharp boundaries, it should be clear from the 
above discussion that sharp reaction rate boundaries have no statistical meaning. In fact, the values 
of the bounds simply depend on the rate probability density function and the associated desired 
coverage probability. To emphasize this point more clearly, we show in Fig. [1] a 3-dimensional reac- 
tion rate probability density function, for the 22 Na(p, 7) 23 Mg reaction (Appendix IC6[) . Although 



we chose (arbitrarily) to associate the 16th and 84th percentiles with the "low" and "high" rates, 
this plot clearly illustrates that the probability density is not zero outside of these bounds. Notice 
that the sharp peak at high temperatures indicates a smaller uncertainty band, whereas, at low 
temperatures, the probability density is much broader because of more uncertain nuclear physics 
input parameters. 

Every experiment performed will be subject to a cutoff at higher energies. In other words, 
if the energy window of effective stellar burning at a given temperature is only partially covered 
by data, then the experimental rate will be smaller than the actual rate. Two questions then 
arise: (i) what is the actual effective stellar energy window? (ii) what is the highest temperature 
at which experimental reactio n rat es are still based on data? These questions are addressed in 



detail in iNewton et al.l (J2008J) and llliadis et al.l (j2010d ). Methods of extrapolating rates to high 



temperatures will be discussed in Sec. I2.4L 



The evaluation of llliadis et al.l (|2010d ) contains 62 reaction rates determined with the Monte 
Carlo method for A = 14 to 40 target nuclei. In this work, we present a new Monte Carlo 
rate for 38 Ar(p, 7) 39 K and updates for the reactions 18 0(p,j) 19 F, 18 0(p, q) 15 N, 22 Ne(p, 7) 23 Na, 
22 Ne(a,7) 26 Mg, 22 Ne(a,n) 25 Mg, 22 Na(p,7) 23 Mg, and 29 P(p,7) 30 S. These results are discussed in 
Appendix O 



2.2.1. Lognormal Reaction Rates 

The reaction rate probability density function (the Monte Carlo output) is determined by the 
interplay of the probability densities for all input quantities. Depending on which contribution 
dominates the total rate, the rate probability density may assume many different shapes. Experi- 
ence has shown, however, that in the majority of cases (but certainly not all) the rate probability 
density is well approximated by a lognormal function. This observation is important for two rea- 
sons, as will be discussed in more detail in Sec. [3j (i) it allows for a convenient implementation of 
the Monte Carlo rate probability density in reaction rate libraries, and (ii) it allows us to adopt a 
lognormal rate probability density for all other reactions for which Monte Carlo rates are not yet 
available. 

The lognormal distribution is a close cousin of the Gaussian distribution and allows for an 
asymmetric sampling of a manifestly positive quantity. Consider a continuous random variable y 
that is Gaussian distributed. If y = ln(x), then x follows a lognormal distribution. In other words, a 
lognormal distribution is simply a distribution where the natural logarithm of a continuous random 
variable follows a Gaussian distribution. It is parameterized by two variables, /x and a: 

f( x ) = ^—- e -^ nx -^ 2 ^ 2 ^ for 0<x<oo. (2) 



However, the distribution is no longer symmetric, and the parameters do not represent the mean 



and square root of the variance, as in the case of the Gaussian distribution. A corollary to the 
central limit theorem illustrates the power of this distribution: the product of a large number of 
distributions of any shape tend toward a lognormal distribution, whereas the sum follows a Gaussian 
distribution. 

The mean, median, and respective uncertainties are given as: 

mean = E[x] = e» +a2/2 , V[x) = e 2t * +a2 (e° 2 - 1) (3) 



median = x me( i = e M , f.u. (68% coverage) = e CT , (4) 

where E[x] is the expectation value and V[x] is the variance (square of the standard deviation) of 
a normalized distribution, f.u. denotes the factor uncertainty. Given values for the expectation 
value and variance, Eq. [3] can be inverted to derive the lognormal parameters, \x and a: 



^ '?&) a = Jln(l + ^ 

EfxlV' V V Efxl 



M = ln(E[x])--ln(l + ^), ff = Jln(l + ^ll) (5) 



For a manifestly positive data set {z{\, Eq. U is derived from the geometric mean and geometric 
standard deviation: 



ix g = $zi ■ z 2 ■ ■ ■ z n = e M , a g = exp[(l/n^(lnzj//i 5 ) 2 ] = e CT . (6) 

i=i 

If the data are lognormal, the quantity e^ is given by the median value. Then, the lower and upper 
bounds are defined as: 



•Plow — l^g/^g — ^ ) %high — P'g^g — ^ i vv 

for a coverage probability of 68%. From either relation, one may determine the factor uncertainty 
such that, 



f.u. = e^+7e^ = e^/e^ a = e° . (8) 

These quantities are crucially important for our Monte Carlo method (Sec. I2.2|) and for using our 
library to calculate rate probability density functions (Sec. [3]). 

We will now discuss the reasons for the observation that most experimental reaction rates 
follow a lognormal probability density. This is illustrated in Fig. [2] for various cases in one sample 
reaction. Panel (a) displays the Monte Carlo reaction rate probability density obtained by assuming 
that only the (non-resonant) direct capture process contributes to the total rate. Since we assumed 
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Fig. 2. — (color online) Monte Carlo reaction rate probability density functions (red histograms) 
for hypothetical input for the 22 Na(p, 7) 23 Mg reaction. The distributions have been obtained using 
5,000 samples at T = 0.3 GK. The solid curves are the lognormal approximations to the rates, (a) 
direct capture (DC) only; (b) single narrow resonance (NR) at Er = 200 keV with large energy 
uncertainty; (c) single NR with large resonance strength uncertainty; (d) single NR with both small 
uncertainties in resonance energy and strength; (e) multiple narrow resonances: Er = 204 ± 1 keV, 
274 ± 1 keV, 434 ± 1 keV, 480 ± 2 keV, 583 ± 1 keV with ujj = 6 ± 1 meV, 39 ± 8 meV, 170 ± 20 
meV, 93 ± 36 meV, 590 ± 70 meV, respectively; and (f) upper limit to a particle partial width 
at the same resonance energy as (a,b,c). The rate distribution in (f) deviates from a lognormal 
distribution since the reduced width input is based on a Porter-Thomas distribution. 



a lognormal distribution for the S-factor input, the reaction rate output distribution is also lognor- 
mal. Shown in panel (b) is the rate probability density for a single narrow resonance with a large 
uncertainty in the resonance energy and small uncertainty in the resonance strength. Because we 
assume that Er is distributed in a G aussian fashion, the B o ltzmann fact or e - E R/ kT that enters 
into the reaction rate formalism (see lLongland et al.l l2010bl ; llliadid 120071 ) will be, by definition, 
lognormal, as is evident in the figure. Panel (c) illustrates the opposite limit case wherein the un- 
certainty in the resonance strength is large for a single narrow resonance with a small uncertainty 
for the reson ance energy. The reaction rate is dir ectly proportional to the lognormally-distributed 
strength (see lLongland et al.l l2010bl ; llliadid 120071 ) and will thus also follow a lognormal distribu- 
tion. Panel (d) shows the reaction rate probability density function for a single narrow resonance 
with small uncertainties in both the resonance energy and strength. In this case, the rate is still 
lognormally distributed but with a small spread parameter, thus appearing more Gaussian. This is 
similar to panel (e), which presents the reaction rate probability density for multiple contributing 
resonances. Here, five resonances were chosen within the Gamow window, with relatively small 
uncertainties. As per the central limit theorem, the sum of these resonances appears Gaussian, 
which can be well-approximated by a lognormal distribution. Panel (f) shows the reaction rate 
distribution for a single resonance with an upper limit on the particle partial width (and, by ex- 
tension, strength). Reduced widths of upper limits are assumed to be Porter-Thomas distributed, 
not lognormal, hence the deviation from the lognormal approximation. Ot her special cases where 
the rate is not lognormally distributed are discussed in llliadis et al.l (|2010d ). 



2.3. Theoretical Reaction Rates: TALYS 



In the following sections, we will discuss and employ theoretical estimates of reaction rates. 
These are important for a number of issues, including: (i) rate estimates for reactions where no 
experimental nuclear physics information exists; (ii) estimation of stellar rates from experimental 
rates, and extrapolation of stellar rates to temperatures at which no experimental rates exist; (hi) 
calculation of stellar enhancement factors, stellar rate ground state fractions, and nuclear partition 
functions (Sec. E2 



In this work, we adopt theoretical estimates for reaction rates based on the statistical (Hauser 



Feshbach) model of nuclear reactions, computed using the code TALYS (JKoning et al. 



2004 



2008 



Goriely et al.ll2008al ). TALYS is a modern nuclear reaction code, which includes many state-of- 
the-art nuclear models to cover the main reaction mechanisms. Recent ly TALYS was updat ed to 
estimate nuclear reaction rates of particular relevance to astrophysics (jGoriely et al.ll2008al ). All 
available experime ntal information on nuclear masses (including the latest release of the 2012 atomic 
mass eval uation bvlWang et al.ll2012r ). deformations, and low-lying states spectra (from the RIPL3 
database. ICapote et al.ll2009l ) are considered. Otherwise, various local and global models are used 
to describe the nuclear structure properties, optical potentials, level densities, 7-ray strengths, and 
fission properties. 



10 



In the special case of radiative neutron captures, Maxwellian-averaged cross sections have been 
measured mostly at 30 keV with relatively hi gh accuracy for ab out 281 stable or long-li ved target 



nuclei . See, in particular, the compilation of iBao et al.l (|2000l ) and references found in IXu et al 



(J2013I ). For those reactions, TALYS (laboratory) Maxwellian-averaged cross sections are normalized 
using the experimental value at 30 keV to ensure an accurate description in the whole temperature 
range, the Hauser-Feshbach being known to provide a proper description of the energy dependence 
of the reaction cross section. Figure [3] displays a comparison of experimental and unnormalized 
TALYS Maxwellian-averaged laboratory (n, 7) cross sections at 30 keV for 281 target nuclei between 
Li and Bi. It is apparent that for medium and heavy target nuclei (A > 50) most of the radiative 
capture rates agree within 40%, although deviations up to a factor of 3 can be found in some 
cases. This r epresents an improv ement in comparison with the former calculation based on the 
MOST code (JArnould et al.l 120071 . see in particular Fig. 39) where globally an 85% accuracy was 
obtained. For light nuclei, larger deviations are found either because the resonance contribution is 
overestimated, or the direct capture contribution has been disregarded. 



A 

D 
v 



A 


v 



10 


•• 








1 


• 

• • 

• 


• 
• 

• 


• 

• 

* 
• 


• • 

• • • 

•• • 

• « - 




• % • .* • 
...» 
.• • • 

• 
• • 


0.1 


• 






- 



50 100 150 

Mass number A 



200 



Fig. 3. — (Color online) Comparison between experimental and TALYS Maxwellian-averaged (lab- 
oratory) (n, 7) cross sections at 30 keV for 281 nuclei between Li and Bi. Experimental data are 
taken fr om the compila tion of IBao et al.l (J2000l ). as well as from more recent measurements com- 
piled in lXu et al.l (|2013l ). Note that for the heavier target nuclei (A > 50), most experimental and 
theoretical values agree within a factor of 2. 



For neutron rich and neutron deficient unstable nuclei, the nuclear reaction rates are esti- 
mated using microscopic models , in particular nuclear structure properties based on the HFB-21 
mass model (IGorielv et al.ll2010T). the Ha rtree-Fock-Bogolvuobov (HFB) plus combinatorial nuclear 
level density model (IGorielv et al.ll2008bl). and the HFB plus Quasi- Random Phase Approximation 
(QRPA) model for the dipole strength (JGoriely et al.ll2004l ). While the particle- induced reaction 
rates are directly computed with the T ALYS code, t he reverse photodisintegration rates are derived 
from the detailed balance relation (see llliadis! 120071 ) . 



11 



Finally, we comment briefly on the predictive power of TALYS rates. The uncertainties in- 
volved are mainly of two origins: (i) the description of the reaction mechanism, i.e., the model 
of formation and de-excitation of the compound nucleus, including a possible pre-equilibrium and 
direct capture contribution, and (ii) the evaluation of the nuclear quantities entering the calcu- 
lation of the transmission coefficients for each entrance or exit channel. Our estimate of factor 
uncertainties for the adopted TALYS rates will be discussed in Sec. 13.41 

TALYS not only provides reaction rate predictions for target nuclei in their ground state 
("laboratory" reaction rates), but also reaction rates that include contributions from thermally 
populated excited target states ("stellar" rates). These topics, including theoretical estimates of 
stellar enhancement factors, stellar rate ground state fractions, and partition functions, will be 
discussed in the next section. 



2.4. Modification of Laboratory Rates for Stellar Environment 

Before the reaction rates discussed above can be used in stellar model calculations, they must 
be corrected for a number of effects. Almost all measurements necessarily determine the reaction 
rate with the target in its ground state, x L = Na{o~v), referred to as "laboratory rate" or "ground 
state rate". However, in stellar environments, excited levels in the target nucleus may be populated 
that can participate in the stellar burning, giving rise to a "stellar rate", x* = Na(o~v)* . This issue 
is particularly important at elevated stellar temperatures. Because there is significant confusion in 
the community regarding how to address this correction, we present a discussion of the different 
methods. 

Consider a reaction + 1 — > 2 + 3, with and 3 denoting the target and final (residual) nucleus, 
respectively. Species 1 and 2 are assumed to be light particles (p, n or a). The laboratory reaction 
rate is obtained by summing over all relevant transitions to excited states, v, in the final nucleus 3: 



x L 



N A (av) = J2 N A(™) 9 - s -^ (9) 



For a non-degenerate plasma in thermodynamic equilibrium, the ratio of the number density iVj M 
of nuclei i in excited state /z and the total number density iVj of nuclei i is given by a Boltzmann 
distribution: 

p _Ni !L g ili e- E ^ kT g lf ,e- E ^ kT 

with g^ = (2jj M + 1), jj M , and Ei^ the statistical weight, spin, and excitation energy, respectively, 
of state fi in nucleus i. T is the plasma temperature. The sum over jjl in the denominator includes 
the ground state and is referred to as the partition function, Gj, of nucleus i. The quantity 
Gf = Gi/gig.s. is called normalized partition function, since gi g . s . refers to the statistical weight of 
the ground state in nucleus i. When experimental information on excited states is not available, 
the partition functions can be obtained from nuclear level density estimates. In this work, we adopt 
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the results computed with the combinatorial model of lGorielv et al.l (|2008bl ) as used in the nuclear 
reaction code TALYS (Sec. 12, 3| ). The stel lar rate is th en obtained by appropriately averaging over 
excited states, /i, in the target nucleus (jlhadis! 120071 ): 



N A (av)* = ^P ^N A {av)^ 



E,g^e- Eo » /kT ZyN A (<Jv)^'- 



(11) 



From these expression we can derive a number of useful quantities. The ratio of stellar to 
laboratory reaction rates, referred to as stellar enhancement factor, is defined by: 



SEF 



N A (av)* _ ^E^o^-^/^E.NAiav 



U-+V 



N A (av) 



Z u N a (<jv)9-°- 



(12) 



Given a theoretical estimate for SEF, a stellar rate can be derived from a laboratory rate via 
N A {o~v)* = SEF N/\(o~v) - Numerical examples for deriving stellar rates from laboratory rates are 
presented in llliadisl (J2007I ). In the present work, we compute the stellar enhancement factors using 
the nuclear reaction code TALYS (Sec. \2.3\i . It is frequently overlooked that the numerical value 
of the stellar "enhancement" factor can be smaller than unity. This may happen, for example, if a 
significant fraction of target nuclei, 0, reside in excited states, and if, at the same time, the reaction 
rates involving these excited states for some reason (e.g., angular momentum, parity, or isospin 
selection rules) are much smaller compared to the ground state rate. For the same reason, one 
should not conclude that a value of SEF = 1 implies a negligible rate contribution from excited 
target states: the interplay of significant excited target state population, Pq^, and small reaction 
rates from excited target states can indeed give rise to a stellar enhancement factor near unity. 

The fractional contribution of the la boratory rate to the s tellar rate, referred to as stellar rate 



ground state fraction, can be defined as (JRauscher et al.ll201ll ): 



GSF 



Pog. s N A {av) 

N A (av}* 



■-E,N A (av) 



303 ■ 

E^E^H^ ~ Go SEF 

a v 



(13) 



The range of possible values amounts to < GSF < 1. For the limiting value of GSF = 1, the 
stellar rate is equal to the laboratory rate, implying SEF = 1. It is apparent that the stellar rate 
ground state fraction contains more information than the stellar enhancement factor. The interplay 
of population, Po^, and reaction rates, N A {av) tJj ^ fV , involving excited target states, //, may give rise 
to a significant overall contribution of excited target states to the total rate (GSF < 1), although 
the corresponding stellar "enhancement" factor may not be affected (SEF = 1). 

Numerical values of SEF and GSF versus temperature for ~ 100 charged-particle induced 
reactions on A < 40 targets, including all reactions with experimental rates estimated using the 
Monte Carlo method, are shown in Fig. 01 With few exceptions, values of SEF range from 0.5 to 
1.8, while GSF exceeds a value of 0.3. 
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Fig. 4. — Values of (top) SEF and (bottom) GSF versus stellar temperature for ~ 100 charged- 
particle induced reactions on A < 40 targets, including all reactions with experimental rates esti- 
mated using the Monte Carlo method. Only those data points are shown for which either value 
deviates from unity (20% of cases). 
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As described in Sec. 12.21 an experimental laboratory rate is typically determined by data over 
a wide range of temperatures. In some instances, for example, neutron captures, an experimental 
rate is only available at a single or a few temperature values only. For the sake of simplicity, 
we will assume in the following discussion that an experimental laboratory rate, x^ x , has been 
measured at a single temperature, T m . Ideally, a theoretical laboratory rate, x^ H , should reproduce 
the experimental value. However, this is rarely the case in practice (Sec. I2.2f) . At this point we 
are faced with a number of questions. Given an experimental value for a laboratory rate at some 
temperature, T m , what is the stellar rate at that very same temperature? What is the stellar rate 
at temperatures not covered by experimental data? And what uncertainties should be adopted for 
the stellar rate? 

To explain the different methods applied so far, it is instructive to chose a simple example. 
Suppose that only a ground state, -Eo^.s. = keV, and one excited state, Eq\ = 50 keV, can 
participate in the stellar burning. For simplicity, we will set the statistical weights, g^^, equal to 
unity. For a temperature of T m = 1 GK, we find kT = 86.173 keV, Go = 1.56, Po g . s . = 0.64, 
and P01 = 0.36, that is, the population probabilities for ground and excited state amount to 64% 
and 36%, respectively. Furthermore, assume that the theoretical ground state (laboratory) and 
excited state rates are xif H = 40.0 cm 3 s" 1 mol" 1 and x\ u = 30.0 cm 3 s _1 mol -1 , respectively. 
Consequently, the stellar rate becomes x^ H = 0.64 x 40.0 + 0.36 x 30.0 = 51.4 cm 3 s" 1 mol" 1 . This 
value is entirely based on theory, and the idea is to scale this value, using a suitable method, if 
experimental information is available for the laboratory rate. Let us assume that the experimental 
laboratory rate amounts to x^ x = 20.0 cm 3 s" 1 mol" 1 . Two methods are applied in practice, 
which are depicted in Fig. [5j 

Method A: In our example, the theoretical and experimental rates are different at T m , as shown 
in the top left of the figure. In the first step, the theoretical laboratory rate at all temperatures, 
Xrpjj (solid blue line), is normalized to the experimental rate, x EX , at T m . This procedure yields 
the dashed blue line in the middle left panel. In the second step, the normalized laboratory rate 
is multiplied by the stellar enhancement factor, which in our case, for example, at T m = 1 GK, 
amounts to {SEF) Tm = {x* TH /x^ H ) Tm = 51.4/40.0 = 1.29 (see Eq. Q2J. The adopted stellar rate 
for this method, x*^, shown as a dashed red line on the lower left, is given by: 

X EX \ (qpTP\ — (~L \ X TH X TH f-> A \ 

-j— (SEF) - (x EX ) Tm (14) 

x TH j Tm x TH yx TH ) Tm 

For our numerical example at T m = 1 GK, we find a value of (x^)T m = 40.0 x (20/40) x 1.29 = 
25.8 cm 3 s _1 mol" 1 . This method is conceptually simple, in the sense that a theoretical stellar 
enhancement factor is applied to a theoretical laboratory rate that is normalized to an experimental 
rate. Its advantage derives from the fact that the quantities SEF = x^ H /x^ H and x^ H /(x^ H )T m 
are ratios of two theoretical model predictions. For example, if a theoretical laboratory rate deviates 
from the experimental one, it is assumed that the theoretical stellar rate is affected by a similar 
deviation, hence lessening the systematic effect on the ratio, SEF. The disadvantage is that no 
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Fig. 5. — Comparison of methods for estimating stellar rates when an experimental laboratory 
rate is known. (Left panels) x^ H denotes a theoretical laboratory (i.e., ground state) rate (blue 
solid line), while Xg X is an experimental rate (data point). The theoretical laboratory rate is first 
normalized to the experimental rate (dashed blue line in middle left) and then multiplied by the 
stellar enhancement factor (SEF), to yield the adopted stellar rate (dashed red line in lower left). 
(Right panels) x^ H denotes a theoretical stellar rate (solid red line), and the fractional ground and 
excited state contributions are displayed as blue and green solid lines, respectively (top right). The 
fractional ground state rate is first normalized to the experimental rate, giving the dashed blue line 
(middle right), and then added to the fractional excited state contribution, to yield the adopted 
stellar rate (dashed red line in lower right). 
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information on the fractional contribution of the laboratory rate to the stellar rate enters in this 
method. 

Method B: The starting point is the theoretical stellar rate (red solid line on the upper right) 
and the fractional contribution of the laboratory (ground state) rate (solid blue line; the solid green 
line shows the fractional contribution of the excited state). For our numerical example, we find for 
the stellar rate ground state fraction at T m = 1 GK a value of (GSF)x m = {Pog.s.x^ H /x^p H )T m = 
0.64 x 40.0/51.4 = 0.50. That is, the ground and excited state each contribute 50% (51.4/2 = 25.7 
cm 3 s _1 mol -1 ) to the total stellar rate at this temperature. In the first step, only the ground 
state rate fraction, (GSF)x^ H , is normalized to the experimental laboratory rate at T m , yielding 
the dashed blue line in the middle right panel. At 1 GK, for example, we find a normalized value 
of {GSF) Tm (x* TH ) Tm ( x ^ x /x^ H ) Tm = 0.50 x 51.4 x (20/40) = 12.9 cm 3 s" 1 mol" 1 . In the second 
step, the normalized rate (dashed blue line) is added to the (unnormalized) excited state fractional 
contribution to the total rate, (1 — GSF)x^ H , shown as green solid line. This sum represents the 
adopted rate for this method, x*^ D , displayed as dashed red line on the lower right, given by: 
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In our numerical example we obtain at 1 GK a value of {x* AD )T m = 0.50 x 51.4 x (20/40) + (1 — 
0.50) x 51.4 = 38.6 cm 3 s _1 mol^ 1 . The advantage of this method is that it accounts for the stellar 
rate ground state fraction. A disadvantage is that it normalizes only the theoretical ground state 
rate to the experimental rate, while the theoretical excited state rate is adopted at face value. In 
other words, if experimental and theoretical laboratory rates deviate from each other, necessitating 
a normalization of the theoretical laboratory rate, then a correlated systematic deviation in the 
excited state fractional contribution to the total stellar rate remains unaccounted for. Furthermore, 
in the extreme case of GSF <C 1, the experimental laboratory rate becomes irrelevant and provides 
no constraint whatsoever on the adopted stellar rate. 



Both methods give identical results if laboratory experimental and theoretical rates agree at 
T m . The same holds for situations where GSF = 1, implying SEF = 1. In both cases, the adopted 
and theoretical stellar rates are equal, x* AD = x^ H . Differences in the two methods arise when 
GSF < 1, as demonstrated i n the numerical example above. Meth od A has been adopted in most 
previous investigations (e.g.. lAngulo et al.lll999i ; IHiadis et al.ll200ll ) and is al so used in the present 



work. Method B has been advocated recently by Rauscher and collaborators (JRauscher et al.ll2011 



Rauscherll2012al lbl). In this context it should be remembered that the higher the temperature, the 
more likely are contributions from excited target states, i.e., GSF < 1. However, at the same time 
reverse reactions become more important an d, in those ca ses where rate equilibrium between a pair 
of forward and reverse reactions is achieved (jlliadisl 120071 ). the actual reaction rates do not impact 
the nucleosynthesis and the distinction between SEF and GSF, or Method A and Method B, is 
unimportant. 



For Method A, the uncertainties in the adopted stellar rate at a temperature of T m , according 
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to Eq. [TH arise from: (i) the experimental rate, Xg X , and (ii) the theoretical quantity SEF. If 
the adopted stellar rate needs to be estimated at temperatures for which no experimental rate 
is available, then an additional uncertainty arises from (iii) the temperature dependence of the 
theoretical ratio, x^ H /(x^ H )T m - For Method B (see Eq. [T5]) . the uncertainties in the adopted 
stellar rate at T m arise from: (i) the normalized partition function, Gq] (ii) the experimental rate, 
Xex'i ( m ) the absolute uncertainty of the theoretical laboratory rate, x^ H ; and (iv) the absolute 
uncertainty of the theoretical stellar rate, x^ H . At temperatures for which no experimental rate is 
available, the uncertainty of the adopted stellar rate has an additional contribution from (v) the 
temperature dependence of the theoretical ratio, x^ H /(x^ H )T rn - 

It should be obvious to the reader that all of the uncertainties estimated from theory are 
difficult to quantify and that either method makes some questionable assumptions. In our opinion, 
it would be worthwhile to pursue another procedure in the future, which we call "Method C". 
Here, the starting point is the basic assumption that nuclear data provide a constraint for the 
nuclear reaction model predictions, no matter the actual fraction of the ground state contribution. 
More precisely, instead of normalizing a theoretical prediction to an experimental reaction rate, the 
nuclear model should be fine-tuned by adjusting nuclear physics ingredients of the Hauser-Feshbach 
calculation, for example, nuclear level densities and 7-ray strength functions (see Sec. l2.3p . in order 
to reproduce both the absolute magnitude and energy dependence of the experimental cross section. 
Once the nuclear reaction model prediction reproduces the nuclear physics data, i.e., the ingredients 
for the calculation of the experimental laboratory rate, the stellar rate predicted by that very same 
model should be more trustworthy. However, this procedure is time consuming. It requires a 
detailed look at each reaction, taking also into account other constraints, for example, the s-wave 
spacing at the neutron binding energy and the Giant Dipole Resonance strength, including the 
corresponding uncertainties. For this reason, Method C has rarely been appli ed for reaction rate 
libraries, although it is com mon practice in some other fields of nuclear science (jCapote et al.ll2009l ; 
Koning &, Rochmanll2009l ). 



Finally, it must be pointed out that at elevated densities the bare nucleus reaction rates dis- 
cus sed above hay e to be corrected for electron screening effects. For more information on this topic, 



sec 



Iliadia (|2007l ). and references therein. 



For stellar model calculations, it is also necessary to include the corresponding reverse rate for 
each forward rate. This rate, Na((tv)Z can be calculated from the forward stellar rate, Na(o-v}% 
using the reciprocity theorem (jlliadisl 120071 ) . From the present results for rates and probability 
density functions of a forward reaction, it is straightforward to calculate, after proper corrections 
for thermal target excitations have been taken into account, the rates and probability density 
functions of the corresponding reverse reaction. We may write the reverse rate as, 



N A (av};=a(e cQ )N A (av)}, 



(16) 



with a equal to a constant factor (see llliadis! 120071 ) and c = — II.6O5/T9. The Q-value is expected 
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to follow a Gaussian probability density function, where the value of Q and the uncertainty AQ 
reported in the literature represent the expectation value and square root of the variance, respec- 
tively (Sec. I2.2p . Therefore, e^ will be represented by a lognormal distribution. In most cases the 
forward rate, NA{crv)*r, is also described by a lognormal distribution, with location and spread pa- 
rameters of \x and o that are tabulated here. Consequently, the reverse reaction rate, Na(ctv}*, will 
follow a lognormal distribution as well, wi th location and spread p arameters of /i r = In a + \i + cQ 
and of = a 2 + c 2 (AQ) 2 , respectively (see lLongland et all ()2010bl ) or the discussion of combining 
lognormal distributions in Sec. l3.3p , The reverse rate factor uncertainties of all Monte Carlo based 
rates are derived in this work from the lognormal spread parameter, according to f.u. r = e°" r . 



3. A New Nuclear Rate Library for Stellar Models: STARLIB 

The ideas presented above are the foundation of a new kind of stellar reaction rate library, 
which we ca ll STARLIB. Other rate libraries co mmonly used fo r stellar models, for example, JINA 
REACLIB (jCvburt et all bold ) or BRUSLIB dXu et al.ll20ia ). provide reaction and decay rates 
versus temperature for all stable and unstable nuclei. In contrast, STARLIB contains informa- 
tion not only on rates versus temperature but also on rate uncertainties and probability density 
functions. With this additional information, more realistic nucleosynthesis simulations can be per- 
formed than previously possible. In particular, STARLIB begins to enable rigorous quantification 
of how improvements in determining stellar nuclear reaction rates impact the structure, evolution, 
and observational consequences of stellar phenomena. 

The new information on rate uncertainties and probability densities can be presented by pro- 
viding only three quantities: temperature, rate, and factor uncertainty. To better understand 
how this can be accomplished, recall the arguments from Sec. 12.2.11 that most reaction rates can 
be described by a lognormal probability density function. The lognormal function has only two 
parameters, fi and a. Applying some of its properties to a reaction rate, x, we find that: 

(i) the median rate, x me( i, is related to the location parameter fi via (combining Eqs.[J]and[6]): 



•Emed 



y^low-Khigh 



(17) 



(ii) the factor uncertainty, f.u., of the rate (for a coverage probability of 
spread parameter a via (combining Eqs.[f3[71 and [8]): 



is related to the 



f.u. = e a = ^x high /xi c 



(18) 



(iii) the high and low rates (for a coverage probability of 68%) are given by (rewriting Eq. [7] 
using Eq. 0]) : 
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•Klow 6 /C %med/ J ' -U- %high 6 6 %med J -W. v.-'-"'/ 

Therefore, by providing temperature, recommended rate (x me d) and factor uncertainty (f.u.), not 
only can the low rate, x; OUI , and high rate, Xhigh, be obtained easily but, in addition, the lognormal 
parameters can be conveniently computed from Eq. [4] or Eqs. [17] and [18J and hence the lognormal 
probability density is determined at any given temperature grid point. We assume a lognormal 
distribution for all reaction and decay rates included in STARLIB. 

Next, the issue of how to present this informa tion in a data library must be addressed. For 
example, the JINA REACLIB (ICvburt et al.l 120101) pr ovides rates versus temperature as an ana- 
lytical fit formula, whereas BRUSLIB (jXu et al.ll2013l ) presents temperature and rates in tabular 
(i.e., two-column) format. While the factor uncertainty could, in principle, be presented as an 
analytical formula, we chose instead the simpler solution of a tabular format. Thus, STARLIB 
consists of three columns, listing, for each reaction or decay, temperature, recommended rate, and 
factor uncertainty. 

It should be clear that a general-purpose nuclear reaction and decay library must encompass 
tens of thousands of nuclear interactions. On the other hand , experiment a l Mont e Carlo rates are 
available so far for only 63 reactions (i.e., those published in llliadis et al.l (J2010d ) and those listed 
in Appendix [C]). In order for the reader to better grasp what kind of information is provided in 
STARLIB, we will discuss the step-by-step construction of STARLIB in the following subsections. 
O ur procedure for trea ting and evaluating the nuclear physics input is described in detail in Sec. 2 
of llliadis et al.l (j2010bl ). Table [2] presents an overview of the nuclear rate sources used in STARLIB, 
including the number of rates and the assumption for uncertainties. STARLIB formatting and 
reference labels are discussed in Appendix [Al 



Table 2: Summary of STARLIB rates and assumption for uncertainties. Numb er indicates forward 
rate only. See Table H] for comprehensive list of rate sources. JINA REACLIB (jCyburt et al.ll2010l ) 
supplements STARLIB for rates not available from sources listed below. 



Rate Evaluation 


Number 


Uncertainties 


Monte Carlo (experimental) 


63 


individual 


NACRE 


24 


individual 


Big Bang 


9 


individual 


/^-decays (experimental) 


2225 


individual 


/^-decays (theoretical) 


5556 


individual 


TALYS (theory) 


47,721 


factor of 10 


Neutron capture 


281 


individual 


26 Al 7-ray transitions 


17 


individual 
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3.1. Starting point: Experimental Monte Carlo Based Reaction Rates 

The three-column structure of STARLIB is motivated by the experimental Monte C arlo based 



thermo nuclear reaction rates, as explained above. These rates consist of those published in llliadis et al 



(J2010d ) and those evaluated in the present work (Appendix [C]). A list of reactions for which Monte 
Carlo based rates are available is given in Table El Notice that in this case the recommended rate, 
listed in column 2 of STARLIB, is directly obtained from the Monte Carlo procedure as the "me- 
dian" rate (Sec. I2.2| ). In particular, no assumption is made regarding the rate probability density 
following a lognormal distribution or not, though assumptions have been made on the input pa- 
rameter distributions. On the other hand, the factor uncertainty, listed in column 3 of STARLIB, 
is explicitly derived under the assumption that the rate distribution is lognormal. 

All experimental "laboratory" Monte Carlo rates are corrected for thermal target excitations, 
according to Method A (see Sec. 12.41 and Eq. [T4"j) . before they are introduced into the rate library. 
The necessary stellar enhancement factors and the extrapolation to higher temperatures are based 
on TALYS calculations (Sec. I2.3|) . For the uncertainties, we have to distingu ish the temperature 



ranges below and above the matching temperature, T matc h (see Sec. 12.21 and iNewton et al.1 12008 



Iliadis et al.ll2010d ). For T < T matc h, we adopt the uncertainties of the experimental Monte Carlo 
rates for the stellar rates. This disregards any additional uncertainties arising from the calculated 
stellar enhancement factors, SEF = x^ H /x^ H (Sec. I2.4p . and, hence, we may underestimate the 
actual uncertainty. However, inspection of the computed SEF values below T = T matc /j (Fig. dj 
top) reveal that they are close to unity, and thus the additional uncertainty is likely to be small. For 
T > T maic /j, we adopt for the stellar rates the uncertainties of the experimental Monte Carlo rate 
at T matc h- This procedure disregards the additional uncertainty from the ratio x^ H /(x^ H )T match , 
which is difficult to quantify at present, as explained in Sec. 12. 4[ From their stellar rates, the rates 
of the corresponding reverse reactions are computed (see Sec. 12.4ft . which are also inserted in the 
new library. 



3.2. Other Experimental Reaction Rates 

Experimental Monte Carlo based reaction rates are available for 63 reactions (Table [3]). How- 
ever, for many other reactions for which experimental nuclear physics information exists, Monte 
Carlo rates have not been computed yet. Thus, the question arises of how to implement such rates 
published in the literature. In the simplest case, when only the experimental reaction rate versus 
temperature has been published and no estimate of a rate uncertainty is provided, the published 
rate was introduced as the recommended rate into our new library, and the factor uncertainty was 
set equal to 10. This represents our best estimate given the absence of information, and the user 
is encouraged to change this if a more appropriate value is suspected. 

A num ber of publica t ions d o pr ovide reaction rates includi ng estimates of uncertainties. For 



example, in lAngulo et al.l (J1999I ) and iDescouvemont et al.l (J2004I ). these uncertainties are expressed 
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Table 3: Reactions with ava ilable experimental M onte Carlo based reaction rates. Most rates and 
Q values are adopted from llliadis et al.l (|2010d ); those that have been evaluated after 2010 are 
indicated by an asterisk and bolded; for more details on the latter reactions, see Appendix ICl 



Reaction a 



Q b (keV) 



Reaction a 



Q b (keV) 



14 C(p, 7 ) 15 N 


10207.42±0.00 


25 Mg(p, 7 ) 26 Al m 


6078.15±0.05 


14 C(a, 7 ) 18 


6226. 3±0. 6 


26 Mg(p, 7 ) 27 Al 


8271.05±0.12 


14 N(a, 7 ) 18 F 


4414.6±0.5 


23 Al(p, 7 ) 24 Si 


3304±27 


15 N(a, 7 ) 19 F 


4013.74±0.07 


24 Al(p, 7 ) 25 Si 


3408±10 


15 0(a,7) 19 Ne 


3529.1±0.6 


25 Al(p, 7 ) 26 Si 


5513.7±0.5 


16 0(p, 7 ) 17 F 


600.27±0.25 


26 A13( P , 7 ) 27 Si 


7462.96±0.16 


16 O(a, 7 ) 20 Ne 


4729.85±0.00 


27 Al(p, 7 ) 28 Si 


11585. 11±0.12 


17 0(p, 7 ) 18 F 


5606.5±0.5 


27 Al(p,a) 24 Mg 


1600.96±0.12 


17 0(p,a) 14 N 


1191.82±0.11 


26 Si(p, 7 ) 27 P 


861±27 


ls O(p, 7 ) 19 F* 


7994.8±0.6 C 


27 Si(p, 7 ) 28 P 


2063±3 


18 0(p,a) 15 N* 


3981.09±0.62 


28 Si(p, 7 ) 29 P 


2748.8±0.6 


18 0(a, 7 ) 22 Ne 


9668.1±0.6 


29 Si(p, 7 ) 30 P 


5594.5±0.3 


17 F(p, 7 ) 18 Ne 


3923.5±0.4 


30 Si(p, 7 ) 31 P 


7296.93±0.19 


18 F(p, 7 ) 19 Ne 


6411. 2±0. 6 


27 P(P, 7 ) 28 S 


2460±30 


18 F(p,a) 15 


2882.15±0.73 


29 P(p, 7 ) 30 S* 


4398.72 ± 3.06 c 


19 Ne(p, 7 ) 20 Na 


2193±7 


31 P( P , 7 ) 32 S 


8863.78±0.21 


20 Ne(p, 7 ) 21 Na 


2431.69±0.14 


31 P(p,a) 28 Si 


1915.97±0.18 


20 Ne(a, 7 ) 24 Mg 


9316.55±0.01 


30 S(p, 7 ) 31 Cl 


290±50 


21 Ne(p, 7 ) 22 Na 


6739.6±0.4 


31 S(p, 7 ) 32 Cl 


1574±7 


21 Na(p, 7 ) 22 Mg 


5504.18±0.34 


32 S( P , 7 ) 33 C1 


2276.7±0.4 


22 Ne(p, 7 ) 23 Na* 


8794.11 ±0.02 c 


31 Cl(p, 7 ) 32 Ar 


2420±50 


22 Ne(a, 7 ) 26 Mg* 


10614.78 ± 0.03 


32 Cl(p, 7 ) 33 Ar 


3343±7 


22 Ne(a,n) 25 Mg* 


-478.29 ± 0.04 


35 Cl(p, 7 ) 36 Ar 


8506.97±0.05 


22 Na(p, 7 ) 23 Mg* 


7580.717 ± 0.708 c 


35 Cl(p,a) 32 S 


1866.21±0.13 


23 Na(p, 7 ) 24 Mg 


11692. 68±0. 01 


34 Ar(p, 7 ) 35 K 


84.5±0.7 


23 Na(p,a) 20 Ne 


2376.13±0.00 


35 Ar(p, 7 ) 36 K 


1668±8 


22 Mg(p, 7 ) 23 Al 


122±19 


36 Ar(p, 7 ) 37 K 


1857.63±0.09 


23 Mg(p, 7 ) 24 Al 


1872±3 


38 Ar(p, 7 ) 39 K* 


6381.43±0.24 c 


24 Mg(p, 7 ) 25 Al 


2271. 6±0. 5 


35 K(p, 7 ) 36 Ca 


2556±40 


24 Mg(a, 7 ) 28 Si 


9984.14±0.01 


39 Ca(p, 7 ) 40 Sc 


538±3 


25 Mg(p, 7 ) 26 Al* 


6306.45±0.05 


40 Ca(p, 7 ) 41 Sc 


1085.09±0.08 


25 Mg(p, 7 ) 26 AF 


6306.45±0.05 







"The ground, isomeric, and total states of 26 A1 are denoted by the superscripts g, m, and t, respectively. 
^Reaction Q values from llliadis et alj ( 2010(3 ). unless otherwise noted. 



c Reaction Q value from W ang et al.l (J2012J) 
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in terms of "lower limits" and "upper limits". Since we assume for all rates in STARLIB that 
the rate probability densities are lognormally distributed, the recommended rate and its factor 
uncertainty are obtained from the geometric mean and geometric standard deviation of the "upper 
limit" and "lower limit" values, according to Eqs. 1171 and 1181 respectively. Notice that Eg. 1181 only 
applies for a coverage probability of 68%, whereas the literature rates have been estimated without 
any reference to probability densities. Therefore, our adopted geometric mean value may differ 
from the reported recommended ratq 2 ], especially if the literature rates have very large uncertainties 
and, in addition, are non-symmetric on a logarithmic scale (i.e., are not lognormal; see Sec. I2.2.T]) . 
However, this difference is found to be very small compared to the large reported rate uncertainties. 
Each of these experimental rates has been corrected for thermal target excitations, according to 
Method A, as explained in Sec. 12.41 (see Eq. ll4p . The necessary stellar enhancement factors, SEF, 
are computed using the TALYS code. The additional uncertainty introduced by SEF has been 
disregarded (see Sec. 13. ip . 



3.3. Experimental Laboratory f3 Decay Rates 

STARLIB lists laboratory /3-decay constants, A = 1n2/ti/2, with t±/2 as the half-life, in units 
of s^ 1 . For a given decay, the same decay constant is listed at all 60 temperature grid points, as 
laboratory decay rates are temperature- independent. In many cases, a /3-decay competes with a 
/3-delayed particle decay. This distinction has not always been made in previous rate libraries but 
is important because these processes represent different links in a reaction network. Consider, for 
example, the nucleus 29 S. Not only does it /3 decay to the ground state of 29 P, it also j3 decays 
to unbound states of 29 P with comparable probability. These excited states, in turn, may decay 
by proton-emission, which results in the final nucleus 28 Si. Consequently, the j3 decay 29 S — s- 29 P 
and the /3-delayed proton decay 29 S — )- 28 Si compete with each other and must be accounted for as 
separate links in nucleosynthesis calculations. 

The partial decay constant, X p , for a given link (i.e., either /3 decay or /3-delayed particle decay) 
is related to the partial half- life, £? /2 ) an d the branching ratio, Bp, via: 



A p = t^ = — ^, (20) 

h/2 J l/2 

where t P x ,„ denotes the total laboratory half-life. Since /3-decay half-lives are manifestly positive 
quantities, we assume that they are associated with a lognormal probability density. It may be 
argued that half-lives are more properly described by a Poissonian distribution, especially if the 
uncertainty is determined by counting statistics. However, almost all of the measured laboratory 



2 We emphasize that any such differences between recommended rates simply reflect our choice of a probability 
density shape and do not imply any indication of quantitative improvement. 
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half-lives have relatively small uncertainties, and thus the Poissonian is well approximated by a 
Gaussian (for sample sizes > 15), which in turn is closely approximated by a lognormal distribution. 
Therefore, we can derive the lognormal parameters fi and a for the partial half- life, i?/ 2 ; from the 
reported recommended value (i.e., the expectation value, E[x]) and its uncertainty (i.e., the square- 
root of the variance, yV^c]), according to Eq. [5j 

To obtain the probability density for the c orresponding partial decay constant, An, we use a 
well-known property of lognormal distributions (JLongland et al.ll2010bl : ICrow &; Shimizulll988l ): for 
two independent, lognormally distributed, random variables x\ and X2, with location and spread 
parameters of /ii, o\ and fi2, o~2, the product ax^x^ (where a > 0) is also lognormally distributed, 
with location and spread parameters of//' = In a + /3i/ii + /32M2 and a' 2 = fi\o\ + {i\o\, respectively. 
From /3i = 0, xi = 1, /?2 = — 1, and X2 = £?/ 2 , we find that the partial decay constant is also 
lognormally distributed, with the median value and factor uncertainty of: 



\ p 



Jin (In 2)-//] 



f.U. 



(21) 



where the lognormal parameters \i and a for the partial half-life are related to E[x] and V[x] 
according to Eq. The latter two quantities are assumed to be equal to the mean value of the 
reported or calculated partial half-life and the square of the associated uncertainty, respectively. 
Therefore, the probability density for the /3-decay constants listed in STARLIB is obtained in 
exactly the same manner as for reactions. That is, the lognormal location parameter // is equal to 
the natural logarithm of the listed partial decay constant, while the lognormal spread parameter a' is 
equal to the natural logarithm of the factor uncertainty. For example, if the reported partial half- life 
amounts to t p , 2 = 10 ± 1 s, we obtain from the above expressions A p = 0.0697 s" 1 and f.u. = 1.105. 
The lognormal parameters of the partial decay constant are // = —2.66 and a' = 0.0998. 



Measured /3-decay half lives and branching r atios for all nuclid e s up to 121 Te are included 
in STARLIB. Most values have been adopted from lAudi et al.l (|2003l ); I Wang et al.l (J2012l ). except 
when newer information was available. Theoretical /3-decay rates were taken from other sources 
(see Table HD. 

For use at elevated temperatures and densities, the /3-decay half lives in STARLIB need to 
be replaced by stellar /3-decay rates. A number of possible complications arise in this regard. 
First, weak interaction rates depend ,in general, on both temperature and density. Therefore, 
libraries that are available to the community (e.g., REACLIB) only include laboratory /3-decay rates, 
and we follow the same strategy in the format of STARLIB. In actual stellar model calculations, 
all laboratory weak rates need to be replaced by temperature and density dependent /3-decay 
ra tes, which are usually contained in a separate input file (based, e.g., on the results presented 



m 



Fuller et al,lll982l ; lOda et al,lll994l ; lLanganke k, Martmez-Pineddl200ll ). Second, when a stellar 



modeler uses any of these "stellar" weak rates from a theoretical model, the question arises of what 
uncertainty to assign to the estimated values. At present, there is no simple answer. The only 
compilation of /3-decay rates that includes a coherent treatment of the temperature and density 
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dependence concerns the work of iTakahashi fc Yokoil (|1987l ). A speci fic detai l ed an alysis of the 
uncertainties related to unknown transition log ft values is given in iGorielvl (| 19991 ) (see Table 
1). The latter work concerns nuclei close to the valley of ^-stability of relevance for s-process 
nucleosynthesis. For nuclei far away from stability, for which no experimental data exist, the 
uncertainties related to the population of excited states are usually smaller compared to those 
associated with the adopted theoretical model (as inferred from different model predictions), and 
more future work on those models is required before reliable uncertainties can be derived. For other 
nuclei, it remains difficult to estimate the uncertainties related to the population of excited states if 
all the transitions and their relative strengths cannot be determined experimentally. Clearly, more 
work is needed in the future to estimate reliable uncertainties of "stellar" weak interaction rates. 



3.4. Theoretical Reaction Rates 

For the vast majority of reaction rates in STARLIB, no experimental cross section information 
exists, and theoretical estimates using the code TALYS have been adopted (see Sec. 12. 3p . These 
theoretical rates have been obtained for target nuclei in the region Z = 3—83. As already mentioned, 
the latest experimental information on Q-values and level structure is adopted in the TALYS 
calculations. Only when no experimental information was available for these input quantities did 
the calculations resort to nuclear mass and structure models (Sec. 12.31) . 



Reliable uncertainties for theoretical reaction rates are difficult to assess. Various claims have 
been made in the literature ("on average within a factor of 2"), which may have been too op- 
timistic. Previously, uncertainties have been systematically evaluated for each target and each 
re action channel from the u se of different sets of nuclear input models (see, in particular, Fig. 22 
of lArnould &: Gorielyil2003l ). A similar approach could be followed to estimate the uncertainties 
affecting the TALYS rates. However, for the present version of STARLIB, we restrict ourselves to 
recommend a factor of 10 uncertainty for any given reaction rate for which no experimental cross 
section information exists. Of course, the theoretical rates are labeled (see Table 0|), and thus it 
is a simple matter for the user to modify the reaction rate uncertainty value. More work will be 
required in the future to refine the estimated uncertainties of purely theoretical reaction rates. 

Since we assume lognormal probability densities for all rates in STARLIB, Eqs. [T7lfT9l ap- 
ply for theoretical rates as well, i.e., the probability densities are easily obtained from the listed 
recommended rate and its factor uncertainty at each temperature grid point. 



3.5. Experimental Neutron Capture Rates 



Exp erimental neutr on capture rates are compiled in lBao et al.l (|2000l ) (see also references pro- 
vided in IXu et al,l (|2013l ). and the KADoNiS project at |www. kadonis .orgp . Most of these experi- 
mental results provide the Maxwellian-averaged cross sections at only a few energies, or sometimes 



25 



even a single energy, near kT = 30 keV, i.e., close to canonical s-process conditions. For a total of 
281 neutron capture reactions involving stable or long-lived target nuclei, the experimental labora- 
tory rates are used to estimate the corresponding stellar rates at all temperatures of astrophysical 
interest, according to Method A in Sec. 12.41 (see Eq. [T4")) . The stellar enhancement factors and 
the extrapolation to temperatures outside the measured region are based on TALYS calculations 
(Sec.ESD- 

For the uncertainties of the adopted stellar rates we took, at all temperatures, only the un- 
certainties of the experimental rates into account, as described in Sec. 13. 1\ although the radiative 
neutron capture rates for those 281 nuclei are also affected by the additional uncertainties related 
to the prediction of the SEF. Such uncertainties can be estimated by a system atic modification of 
the nuclear input models, as discussed above (see also lArnould fc Gorielyll2003l ). They remain how- 
ever relatively small and are neglected in the present version of STARLIB. Again, this procedure 
underestimates the actual uncertainties, and more work is required in the future (see Method C in 
Sec.GG 



Probability densities for the neutron capture rates listed in STARLIB are obtained, as for all 
other reactions, using Eqs. [TTlfl"9l 



3.6. Al /3-Decays and 7-Ray Transitions 

The nuclide 26 A1 has an isomeric state (E x = 228 keV; J w = + ) that does not always achieve 
thermal equilibrium with the ground state (J n = 5 + ) in a stellar plasma. Therefore, the ground 
and isomeric state need to be treated separately in nuclear reaction networks, d epending on the 
tempe ratures involved. Stellar modelers generally follow the recommendation of IWard fc Fowler 
(| 19801 ): the ground and isomeric states in 26 A1 will be in thermal equilibrium at temperatures above 
~ 400 MK and can be treated as one species in the network. Below this temperature, the ground 
and isomeric states are not in thermal equilibrium and should be adopted as distinct species. 



This e xtreme assumptio n based on a sharp temperature boundary was questioned by lRunkle et al 



(J200ll ) and lCoc et al.1 (|1999l ). who studied the equilibration of ground and isomeric state via 7-ray 
excitations involving higher-lying levels of 26 A1. They showed that it is sufficient to consider just 
three excited levels for the equilibration: E x = 417 keV (J 71- = 3 + ), 1058 keV (J 77 = 1 + ), and 
2070 keV (J n = 2 + ). STARLIB includes all levels discussed above, and in a nucleosynthesis sim- 
ulation the user may select only a single species (i.e., 26 A1 in thermal equilibrium, presumably 
because temperatures are very high), two species (i.e., the ground a nd isomer, assuming that the 



equilibration can be disregarded; see, for example. Illiadis et al.ll201ll ). or all five species (i.e., if the 



equilibration needs to be taken into account explicitly). The safest method of implementation in a 
reaction network is to consider all five species such that the equilibration of 26 A1 is always explicitly 
incorporated. However, processing time will increase significantly. 



The necessary 7-ray transition and /3-decay constants are listed in Appendix A of Illiadis et al. 
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(|201ll ). and all of these transition rates are included in STARLIB. A factor uncertain ty of 10 was 



estim ated for those transition rates that have been derived from nuclear theory (see llliadis et al 



20111 . for details). 



3.7. Other Rates 

For other rates not available from the above so urces, STARLIB is supplemented by rates 



from the JIN A REACLIB library (jCvburt et al.ll2010l ). Their reaction rates, provided as a seven 



parameter fit, were converted into a tabular format, and a factor of 10 uncertainty was assumed 
for each rate, in absence of more reliable information. 



4. Use of STARLIB in Monte Carlo Nucleosynthesis Studies 

We will briefly comment on how to randomly sample reaction and decay rates in a Monte Carlo 
nucleosynthesis calculation, since there seems to be no consensus on a procedure in the literature. 
We anticipate that such studies will become increasingly important in the future, allowin g for more 



realist ic elemental and isotopic abundance estimates. We will follow here the ideas of lLongland 



(J2012J), to which the reader is referred for more details. 



Recall that STARLIB provides for each nuclear interaction (nuclear reaction, weak interaction) 
a lognormal rate probability density function at 60 temperature grid points between 1 MK and 10 
GK: the lognormal parameter [i can be approximated by the natural logarithm of the recommended 
rate listed in column 2 (Eq. [T7l) , and the lognormal parameter a is obtained from the natural log- 
arithm of the rate factor uncertainty listed in column 3 (Eq. [T8l) , With these two parameters, 
the lognormal probability density is precisely defined (see Eq. [2]) . In a Monte Carlo nucleosynthe- 
sis simulation, all rates can be simultaneously sampled according to their individual probability 
densities, or the rates for a subset of links can be sampled, depending on the application. 



For a lognormal probability density, the samples are drawn according to the function (j Evans et al 



20001 ): 



Xl = e ^+Pi- = e » ( e *)« = x med (f-u-) Pi (22) 

where Xi denotes the sampled rate, and pi is a standard normal deviate (i.e., a Gaussian deviate 
with a mean of zero and a standard deviation of unity). See also Eqs. [T7] and [THJ For a random 
sample of pi = 0, the recommended rate (e^) is obtained. It is apparent from the right hand side of 
Eq. [22]that the first term (e^) is listed in column 2 of STARLIB (i.e., the recommended rate), while 
the base (e a ) is listed in column 3 of STARLIB (i.e., the factor uncertainty). Thus, it becomes 
obvious how the structure of STARLIB is conveniently tailored for Monte Carlo nucleosynthesis 
simulations. 
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In practical terms, we find it useful to sample the rates according to Eq. [22] only once at 
the beginning of each network calculation. That is, in the beginning of each network run (for 
example, post-processing studies using constant or varying temperatures and densities, or multi- 
zone calculations), the probability factor pi is sampled for each reaction, i, in the network, where 
for a given reaction pi has the same value at all temperatures. Notice that this procedure does 
not imply that, for a given reaction, the rate variation is constant with temperature. Rather, the 
lognormal parameter a depends on the temperature and, therefore, randomly sampling according to 
Eq. [25] automatically takes into account the temperature dependence of the rate factor uncertainty. 
Furthermore, the parameters Pi could be saved for each sample reaction network run in order to 
study correlations between rates and abundances. 

Note that the procedure described above is dif ferent from Monte Carlo variations in previous 
studies. Simulations of big bang nucleosynthesis (JKrauss &: Romanellil ll99Q : ISmith et al.l Il993l ) 
employed Gaussian rather than lognormal rate probability densities. Their assumption is only 
approximately correct when the rate uncertainties are relatively small (say, less than 20%) for 
all interactions in the network. Lo gnormal rate probability d ensities were used in simulations 



of cla ssical novae ([Smith et al.l [2002]) and Type I x-ray bursts ([Roberts et al.l 120061 ; iParikh et al. 



2008j). In these latter works, the assumed enhancement factor multiplying each rate, i.e., the factor 
e api , instead of pi in Eq. (22] was sampled randomly. That procedure is applicable only for the 
special case that the rate uncertainty is constant with temperature, for example, if rigorous rate 
uncertainties are not available and global factor uncertainties (e.g., "factor of 10") need to be used 
in the simulations instead. On the other hand, our method described above, of sampling the rate 
according to Eq. [22] using the information provided by STARLIB, is more generally applicable, i.e., 
if the rates of a given reaction have large uncertainties or if they vary with temperature. For more 
in yolved methods o f rand omly sampling rates in nucleosynthesis simulations, the reader is referred 



to 



Londand et al.l (J2012h . 



An example is shown in Fig. [6] displaying results of a Monte Carlo post-processing reaction 
network simulation for the evolution of a 25 Mq star with initial solar composition from the end of 
hydrogen burning to the end of helium burnin g. The tim e evolu tion of temperature and density in 
the stellar core is adopted from the models of iThe et al.l (|2000l ). All rates are sampled randomly, 
simultaneously, and independently according to Eq. [22j As explained above, this is achieved by 
drawing random variables pi for each reaction at the beginning of each sample network calculation 
and by using the lognormal parameters fi and a provided in STARLIB. The results are obtained for 
1500 reaction network samples. The main part of the figure displays the 26 Mg abundance versus 
the random variable pi for the 22 Ne(a,7) 26 Mg reaction rate. There is a noticeable slope in the 
data scatter, indicating that the 26 Mg abundance at the end of helium burning varies as a result of 
changing the 22 Ne(a,7) 26 Mg reaction rate. The inset shows the resulting distribution of final 26 Mg 
abundances, which is given by a projection of the scatter plot onto the y-axis. The location and 
width of this distribution can be used to determine a mean value and an uncertainty for the final 
26 Mg abundance. The scatter plot also indicates that the 22 Ne(a,7) 26 Mg reaction rate is not the 
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Fig. 6. — (color online) Results of Monte Carlo post-processing nucleosynthesis study for the 
evolution of a 25 M Q star of initial solar metallicity from the end of core hydrogen burning to the 
end of cor e heliu m burning. Temperature and density evolutions are extracted from the models of 
The et al.l (|2000l ). The results are obtained for 1500 reaction network samples. (Inset) Distribution 
of 26 Mg abundance at the end of helium burning, where the spread in the values provides a realistic 
estimate for the abundance uncertainty; (Main) final 26 Mg abundance versus random variable pi 
for the 22 Ne(a,7) 26 Mg reaction rate (see Eq. [22]). The impact of the 22 Ne(a,7) 26 Mg rate on the 
final 26 Mg abundance is apparent. 



only important factor determining the 26 Mg abundance (otherwise there would be less scatter in 
the y-direction) . 

In the example above we assumed that all rates are independent. However, correlations cer- 
tainly exist between reaction rates. For example, the rates of the competing 22 Ne(a,7) 26 Mg and 
22 Ne(a,n) 25 Mg reactions are computed using the same values of particle separation energies and 
partial widths since the same compound levels are populated. Therefore, these two rates are indeed 
correlated. Such input parameter correlations may indeed be important for the nucleosynthesis, 
and it would be of interest to study their impact. As of yet, they are not included in the STARLIB 
format. Furthermore, although STARLIB contains probability densities for all reverse rates, these 
should in general not be sampled independently because forward and reverse rates are correlated 
according to detailed balance. A proper Monte Carlo reaction network study would involve the 
following steps: (i) sampling over the rate of the forward reaction, as described above; (ii) sampling 
over the reaction Q-value, which is described by a Gaussian probability density with a mean of 
Q and a standard deviation of AQ; (iii) calculating the resulting reverse rate from Eqs. 3.44 and 
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3.45 of llliadia (|2007l ). and Eq. [16] for the rate of the forward reaction that was sampled under 
step (i). Disregarding the correlation between forward and reverse rate may greatly exaggerate the 
Monte Carlo derived abundance uncertainties at elevated temperatures, where reverse rates become 
important 







For the reverse rate sampling method described above, an estimate for the Q- value uncertainty 
is needed. Since this information is not provided in the present version of STARLIB, the user needs 
to supplement STARLIB with other sources. Future releases will include the Q-value uncertainties. 



5. Conclusions 

In this work, we have presented STARLIB, a next-generation reaction-rate library for nuclear 
astrophysics. This library fills a void in the community and urges a paradigm shift for stellar 
modelers to incorporate not only rate uncertainties but also their probability densities, as no cur- 
rent library provides this information. We have outlined the Monte Carlo procedure that enables 
assignments of rigorous uncertainties for experimental reaction rates. Monte Carlo rates are sup- 
plemented by a variety of other types to facilitate an all-purpose nuclear reaction library, including 
weak rates and theoretical rates for which no experimental information is known. Methods to cor- 
rect laboratory rates for a stellar environment have been discussed, and an application of STARLIB 
in a Monte Carlo nucleosynthesis study was demonstrated. 

This work was supported in part by the National Science Foundation under award number 
AST-1008355 and by the Department of Energy under grant number DE-FG02-97ER41041. We 
would also like to thank Richard Longland, Jordi Jose, Alain Coc, and Matthew Buckner. In 
addition, we thank Arjan Koning and Stephane Hilaire for their valuable contribution in writing 
the TALYS code and making it publicly available. Our collaborators at ASU acknowledge partial 
support from NASA. 

The following appendix details STARLIB's formatting, its website, and new and updated 
Monte Carlo rates. 



3 For the example shown in Fig. 6, the effects of the correlations discussed here are small. First, the main 
uncertainties in the 22 Ne+a rates at s-process temperatures are determined by the contributions from unobserved 
resonances, which are uncorrelated. Uncertainties from correlated quantities, such as excitation energies, properties 
of observed resonances, etc., are small. Second, s-process temperatures are typically < 300 MK, which is too low for 
a significant abundance flow via reverse reactions to occur. 
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A. STARLIB Format 

For each nuclear interaction, STARLIB lists a single-line header, followed by three columns: 
temperature (in GK), recommended rate (cm 3 s _1 mol -1 for reactions; s~ 1 for photodisintegration 
and /3-decay constants), and factor uncertainty of the rate. The temperature grid consists of 60 
values between 1 MK and 10 GK. STARLIB reference labels are listed in Table HI 

Table 4:: Labels and references in STARLIB 

Label Reference 

REACTION RATES BASED ON EXPERIMENT 

cf88 Caughlan & Fowler, At. Data Nucl. Data Tab. 40, 283 (1988) [CF88] 

de04 Descouvemont et al, At. Data Nucl. Data Tab. 88, 203 (2004) 

mclO Iliadis et al, Nucl. Phys. A 841, 31 (2010) [Monte Carlo rates] 

mcl3 Present work [updated Monte Carlo rates]; see Appendix IU1 

nacr Angulo et al, Nucl. Phys. A 656, 3 (1999) [NACRE] 

taex Bao et al. At. Data Nucl. Data Tab. 76, 70 (2000), and updates; present work 

[extrapolated using TALYS]; see Sections 12.31 and 13.51 
REACTION RATES BASED ON HAUSER-FESHBACH THEORY 
taly Present work [calculated using TALYS]; see Sections 12.31 and 13.41 
DECAY RATES BASED ON EXPERIMENT 
au03w Audi et al., NPA 729, 3 (2003) 
auecw Audi et al., NPA 729, 3 (2003)1 
bet-w /3~-decay 
bet+w /3 + -decay 

ecw pep reaction; 3He — > t decay 1 
ru09w Rugel et al., Phys. Rev. Lett. 103, 072502 (2009) 
wcl2w Tuli, wallet cards, National Nuclear Data Center (2012) 
willw Wietfeldt and Greene, RMP 83, 1173 (2011) 
DECAY RATES BASED ON THEORY 

bkmow Klapdor, Metzinger k Oda, At. Data Nucl. Data Tab. 31, 81 (1984) [/3~-decay] 
btykw Takahashi, Yamada & Kondo, At. Data Nucl. Data Tab. 12, 101 (1973) [/3+-decay] 
illlg Iliadis et al., ApJS 193, 16 (2011) [7-ray transitions] 

ka88w Kajino et al, Nucl. Phys. A 480, 175 (1988) [/3-decay of excited 26 A1 levels] 
mo92w Moeller et al. (1992) [/3~-decay] 

mo03w Moeller, Pfeiffer and Kratz, Phys. Rev. C 67, 055802 (2003) 
INDIVIDUAL RATES, MAINLY BASED ON EXPERIMENT 
an06 Ando et al, Phys. Rev. C 74, 025809 (2006) 
arl2 Arnold et al, Phys. Rev. C 85, 044605 (2012) 



'Rates need to be multiplied by the electron density, pY e . 
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Table 4 - continued 



Label Reference 



bb92 Rauscher et al., ApJ 429, 499 (1994) 

beOl Beaumel et al, Phys. Lett. B 514, 226 (2001) 

bul2 Buckner et al, submitted to Phys. Rev. C (2012) 

cy08 Cybert & Davids, Phys. Rev. C 78, 064614 (2008) 

fu90 Fukugita & Kajino, Phys. Rev. D 42, 4251 (1990) 

halO Hammache et al, Phys. Rev. C 82, 065803 (2010) 

illl Iliadis et al, Astrophys. J. Suppl. 193, 16 (2011) 

im05 Imbriani et al, Eur. Phys. J. A 25, 455 (2005) 

ku02 Kunz et al, Astrophys. J. 567, 643 (2002) 

lol2 Longland et al, Phys. Rev. C 85, 065809 (2012) 

maf o Malaney & Fowler, Astrophys. J. 345, L5 (1989) 

nk06 Nagai et al, Phys. Rev. C 74, 025804 (2006) 

re98 Rehm et al, Phys. Rev. Lett. 80, 676 (1998) 

rolf C. Rolfs and collaborators 

ta03 Tang et al, Phys. Rev. C 67, 015804 (2003) 

wies M. Wiescher and collaborators 



The header consists of one line of 64 fields, indicating (i) the type of interaction, (ii) the interact- 
ing nuclei, (iii) a reference label, (iv) an interaction specific label, and (v) the recommended Q- value. 



(i) Type of interaction (left-aligned in fields 1-2): these follow the JINA REACLIB chapter num- 
bers, specifically their "format 2, with chapters 9-11" version (see lREACLIBll2012l ). They are in 
principle redundant since the type of interaction can always be reconstructed from the interaction 
label (see below). However, they have been included in STARLIB as convenient labels for some 
users. The chapter numbers are (e^ stand for interacting nuclei): 



1 ei — > e 2 

2 ei — > e 2 + e 3 

3 ei — > e 2 + e 3 + e 4 

4 ei + e 2 — > e 3 

5 ei + e 2 — > e 3 + e 4 

6 ei + e 2 — > e 3 + e 4 + e 5 

7 ei + e 2 — > e 3 + e 4 + e 5 + e 6 

8 e 4 + e 2 + e 3 — > e 4 

9 e 4 + e 2 + e 3 — > e 4 + e 5 

10 e 4 + e 2 + e 3 + e 4 — > e 5 + e 6 

11 e 4 — > e 2 + e 3 + e 4 + e 5 
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(ii) Interaction label (fields 6-35): all species in the entrance channel are listed first (lightest to 
heaviest), followed by species in the exit channel; each species in the interaction occupies 5 fields 
(right-aligned), with the element symbol given first, followed by the mass number; exceptions are 
neutrons ("n"), protons ("p"), and tritons ("t"); levels in 26 A1 (see Sec. I3.6[) have special labels: 
E x = keV (ground state; 5+), 228 keV (isomeric state; 0+), 417 keV (3 + ), 1058 keV (1+) and 
2070 keV (2+) are denoted by "al-6", "al*6", "alOl", "al02", and "al03", respectively; the label 
"al26" is used to denote 26 A1 with all levels in thermal equilibrium (presumably at sufficiently high 
temperatures). 

(iii) Reference label (right- aligned in fields 44'4 r O : these labels, which usually have a length of 
four characters, provide information on the source of the rates; a complete list of STARLIB ref- 
erence labels is given in Table 01 Note that the characters "ec" in fields 46-47 indicate electron 
capture, and, therefore, these decay rates need to be multiplied by pY e , where p and Y e denote the 
density and electron mole fraction, respectively. 



(iv) Interaction specific label (field ^8): indicates a weak interaction ("w"), a 7-ray transition 
("g"), or if the reaction rate has been calculated from a forward rate using the principle of detailed 
balance ("v"). 



(v) The Q-value (fields 53-64): indicating the energy re lease in the n u clear interaction, in units o f 
MeV; values are obtained from the recommendations of llliadisl (|2007l ) ; iThompson &: Iliadisl ()1999l ) 
or, if no experimental information is available, from nuclear models (see Sec. I2.3[) . 



A number of examples are listed below: 



1 


ol5 nl5 




1 


al-6 alOl 




4 


P P 


d 


4 


P P 


d 


5 


he4 nl4 


P 



ol7 



au03w 


2 


75420e+00 


illlg 


-4 


16800E-01 


nacr 


1 


44222E+00 


ecw 


1 


44206e+00 


mclOv 


-1 


19182e+00 



The first example denotes the /3-decay 15 0— » 15 N, "au03" indicates that the decay rate was 
adopted from lAudi et al.1 (120031 ). the "w" label indicates a weak interaction, and the total energy 
release (not corrected for neutrino losses) amounts to 2.7542 MeV. The second example illustrates 
the 7-ray transition from the gro und state of 2 6 A1 ( 5 + ) to the third excited state at E x =417 keV 
(3 + ), listing the source of the rate llliadis et al.l (|201ll ) and the label "g" for a 7-ray transitions, and 
providing the energy for the tran sition (-0.4168 MeV ) . The third example refers to the p + p — > d 
reaction, with the source listed as lAngulo et al.l (|1999l ). and an energy release of 1.44222 MeV. The 
fourth example illustrates the p + e+p + v — > d reaction, which depends on the electron density in 
the plasma and thus is distinct from the previous example. In this case, the label "ec" indicates that 
the reaction rate needs to be multiplied by pY e in a reaction network calculation, similar to a weak 
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interaction (hence the label "w"). The fifth and last example denotes the 14 N(g, p ) 1 O re action. 



The source of the reaction rate is the Monte Carlo based evaluation of llliadis et al.l (|2010d ). with 
the label "v" (reverse) indicating that this rate has been obtained from the rate of the forward 
reaction, 17 0(p, a) 14 N, using detailed balance. The reaction Q-value amounts to -1.19182 MeV 
(endothermic reaction) . 

Below we summarize the meaning of each of the three columns that follow the single- line header: 

(i) Column 1: Temperature, in GK (10 9 K). 

(ii) Column 2: Recommended rate, either the thermonuclear reaction rate in units of cm 3 mol^ 1 
s" 1 , or the decay constant of a photodisintegration or a /3-decay in units of s" 1 . 

All rates in STARLIB are assumed to be lognormally distributed (see Sec. I2.2.ip . however, only 
for the Monte Carlo based rates do the listed rate and factor uncertainty have a precise statistical 
meaning. See Sees. l2~?2l and 151 for a detailed discussion. The lognormal probability density is given 
by Eq. m 

f(x) = —^=-e- {lnK -rf /{2a2) for 0<x<oo, (Al) 



where x is the reaction rate. The lognormal parameters are related to the median, low, and high 
rates via X[ ow = e M_cr , x me d = e^, and Xhigh = e M+cr , respectively, for a coverage probability of 68%. 
Thus, the natural logarithm of the recommended rate provides the lognormal parameter \x. 

(iii) Column 3: Factor uncertainty, f.u. One may multiply and divide the recommended rate 
(Column 2) by the factor uncertainty to calculate the high and low rates. The factor uncertainty 
is related to the shape parameter, a, by f.u. = e a . For small values of a, below 0.1, the lognormal 
distribution appears more Gaussian (i.e., symmetric) in shape and deviates from symmetry as a 
increases in value. The natural logarithm of the factor uncertainty yields the lognormal parameter 
a . 



Once both the lognormal parameters ix and a have been determined by the natural logarithms 
of the values listed in Column 2 and 3, respectively, the probability density function may be 
calculated using Eq. [21 



B. STARLIB Website 

The STARLIB library is available at: 
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starlib . physics . unc . edu 



where one may display an individual rate, download the entire library, or download a user-specified 
subset of the library. The website allows ease of dissemination into the community and a portal 
for real-time updating of new rates. Also included is a reaction rate calculator, where one may use 
Monte-Carlo simulation code RatesMC to calculate an experimental reaction rate. The following sec- 
tions outline each major function of the site. Details on STARLIB, its construction, hyperlinked ref- 



erences, and information on updates may be accessed at starlib . physics . unc . edu/Details . php 



Step-by-step instructions for the website are shown in lSallaska &: Iliadisl (120121 ) 



B.l. Library 

The entire STARLIB library can be downloaded from: 



starlib . physics . unc . edu/RateLib . php 



As we update new versions of STARLIB, archived versions will also be available for comparison 
purposes. Because some stellar models may not require the full STARLIB library, one may also 
download a subset of the library for use in a reduced network, after uploading a user-specified list 
of nuclides. The formatting of the list may have the chemical name before or after the mass number 
(for example, "na22" or "22na") and is case insensitive. The site will accept either left or right 
aligned lists. 

Individual rates may also be displayed with the current or archived versions of the library. 
There are three available menus from which to choose. The left menu requires all interacting nuclei 
to be entered (i.e., electrons or 7-rays should be omitted), up to three initial and three final species. 
This menu also allows a choice of special reactions, for example, those involving four interacting 
nuclei or the p + e + p — > d reaction. 

The middle menu allows input in the form of a target nuclide and reaction. The reaction may 
be formatted as "x,y" or "(x,y)". A complete list of formatting options is shown in Table [5j 

The right menu is for /3-decays and 7-ray transitions, where parent and daughter nuclei may 
be entered. There is also a button that will display examples of the various options for each menu; 
multiple clicks cycle through examples. 

Upon submission of any of the three menu options, the individual rate will be presented 
along with its uncertainty, its source (including a link to the source description and Digital Online 
Identification), Q-value, and the library version. 
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Table 5: Formatting options for user input on the STARLIB website. 



Nuclide 


Format 


26 


Al 


Reactions 


proton 


P 


hi lh 


al26 


al*6 


x,y 


deuteron 


d 


h2 2h 


ai-6 


alOl 


(x,y) 


triton 


t 


h3 3h 


al02 


al03 




a-particle 


a 


he4 4he 








neutron 


n 


nl In 








general 


- 


na22 22na 









B.2. Calculator for Experimental Monte Carlo Based Reaction Rates 

In addition to ac cessing the STARLIB library, the Monte Carlo code RatesMC (written by 
Richard Longland; see lLongland et al,ll2010bl ) is available to users to run on our website at: 



starlib . physics . unc . edu/RateCalc . php 



We currently employ an Intel Xeon Processor X5650 server with dual 6-core, hyper-threaded, 2.66 
GHz processors (maximum of 24 simultaneous simulations of RatesMC). 

The first set of menu items on the Rate Calculator webpage allows access to a repository of 
current input files for experimental Monte Carlo rates. The same formatting applies as outlined 
in the previous subsection, as well as displaying examples via the "?" buttons. For the reaction 
25 Mg(p, 7) 26 A1, options for choosing the isomeric state, ground state, or both combined are avail- 
able. This database is available freely to all users. However, to run simulations, we require users 
to register their email addresses. 

We recommend that the user edit an existing input file from our online database. Once the 
input file is uploaded or edited, the Monte Carlo simulation may be started by logging in with 
a registered email address. Our server employs a PBS queue system, and, after submission, the 
website displays the job ID, status ("R" for running, "Q" for in the queue), the run time, simulation 
progress (in percent), and estimated time remaining. The site employs several shell scripts, one 
of which determines if any resonances are being integrated. Depending on the circumstance, the 
jobs are filtered into different queues to allow jobs without integrated resonances to be completed 
immediately and instantaneously. A Monte Carlo simulation requires significantly more CPU time 
if resonance contributions to the total rate have to be determined by numerical integration. As 
discussed in Sec. 12. 2\ a (p, 7) reaction with one integrated resonance will take approximately 45 
minutes on our current server. Simulations may be aborted by supplying the job ID. 

When the simulation is complete, an email alert will be sent to the address provided, including 
the job ID and a link to access the output files. Files may also be accessed from: 



starlib . physics . unc . edu/RateCalc . php#output 
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or by scrolling to the bottom of the Rate Calculator page to log in. Files are delineated by the 
reaction name, the job ID, and the appropriate file extensions (".in", ".out", ".hist", ".latex"). For 
ease of bookkeeping, both the input file (".in") and output file (".out") are in the data repository. 
The ".hist" file includes the binned reaction rate probability density function for each of the 60 
temperature values, and the ".latex" file contains the reaction rate table in Latex format. Files 
may be sorted by name, size, or date and may also be deleted by specifying the job ID. 



C. Updated Experimental Monte Carlo Based Rates 

Here we provide updated Monte Carlo based thermonuclear rates for the reactions 18 0(p, 7) 19 F, 
18 0(p,a) 15 N, 22 Ne(p, 7 ) 23 Na, 22 Ne(a, 7 ) 26 Mg, 22 Ne fa,n) 25 Mg, 22 Nafa,7) 2 3 Mg, 29 P(p,7 ) 30 S, an d 
new rate 38 Ar(p, 7) 39 K. Similarly, to the procedure in lLongland et al.l (J2010bl ); llliadis et al.l (J2010a lb|jal). 
we provide for each reaction (i) a rate table, (ii) the complete input file used to calculate the rate; 
(iii) the figure with the rate probability densities at six sample temperatures; and (iv) a figure 
comparing the new Monte Carlo based rate with the latest previously reported rate. Some of 



the r eaction rates listed here have been published elsewhere (jLongland et al.l 12012 : iBuckner et al 



2012j), but the information according to (ii), (iii), and (iv) has not been reported elsewhere. Details 
are given below. Energies are in the center-of-mass system, unless otherwise noted. 

(i) Reaction rate table: 



T 

Low rate 
Median 
rate 
High rate 

/' 
a 
A-D 







Temperature in GK. 

0.16 quantile of the cumulative reaction rate distribution. 

0.50 quantile of the cumulative reaction rate distribution. 

0.84 quantile of the cumulative reaction rate distribution. 
All reaction rates are in units of cm 3 s _1 mol -1 . 

Parameter determining the location of the lognormal reaction rate probability density 
function. 

Parameter determining the width of the lognormal reaction rate probability density func- 
tion. 
Anderson-Darling test statistic, t* AD , indicating how well the Monte Carlo reaction rates 



are ap proximated by a lognormal distribution (see lLongland et al.l l2010bl : llliadis et al. 



2010c 



Oft. 



Values given in parenthesis are not obtained from the Monte Carlo method but are found 
from extrapolation to elevated temperatures (see Sec. 12. 2p . 



(ii) Nuclear data input table: an example input file and the corresponding detailed explanations 
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are reproduced from Appendix A of llliadis et al,l (|2010bl ). and are repeated here for convenience. 



03 


1 


04 


8 


05 


2 


06 


1.0078 


07 


16.999 


08 


4.0026 


09 


0.5 


10 


1.5 


11 


0.0 


12 


5261.3 


13 


4015.3 


14 


1.25 


15 


3 


16 


****** 


17 


1.0 


18 


5000 


19 






01 17X(p,a)14Y 

Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 

Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures; =NT for rate output at selected temperatures 

21 Nonresonant Contribution 

22 S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

23 3.1e2 -2.1e-l 4.5e-6 0.4 1200.0 

24 0.0 0.0 0.0 0.0 0.0 

26 Resonant Contribution 

27 Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg.Gx in (eV) ! ! 

28 Note: if Er<0, 

29 Ecm DEcm wg 

30 -5.2 0.6 

31 477.3 1.1 

32 510.5 2.0 6.0e-l 

34 Upper Limits of Resonances 

35 Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 

36 Note: . . .PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 

37 Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 
38-3.4 0.7 1 7.9e-3 0.0 1 0.0045 15.0 2.9 1 0.0 0.33 0.07 1 0.0 1 
39 73.5 0.4 1 4.5e-7 1.8e-7 1 0.0 150.0 0.0 1 0.010 0.21 0.02 1 0.0 1 

41 Interference between Resonances [numerical integration only] 

42 Note: + for positive, - for negative interference; +- if interference sign is unknown 

43 Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 

44 +- 

45-2.3 0.6 1 6.2e-3 0.0 1 0.0045 17.0 1.1 1 0.0 0.71 0.07 1 

46 89.0 0.8 1 2.9e-9 0.8e-9 1 0.0 230. 09. 1 0.0 0.66 0.03 1 

48 Comments : 



Dwg 


Jr 


Gl 


DG1 


LI 


G2 


DG2 


L2 


G3 


DG3 


L3 


Exf 


Int 





1 


0.022 


0.010 


2 


30.8 


2.6 





0.77 


0.21 


1 


0.0 


1 





4 


140 


32 


1 


110 


12 


3 


0.13 


0.02 


1 


0.0 


1 


1.8e-l 
































0.0 






PT 






Exf 

0.0 
0.0 
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49 1. Narrow resonance information is adopted from Peterson et al . 1973. 



Row 01 




Row 02 




Row 03 




Row 04 




Row 05 




Row 06 




Row 07 




Row 08 




Row 09 




Row 10 




Row 11 




Row 12 




Row 13 




Row 14 




Row 15 




Row 16 




Row 17 




Row 18 




Row 19 




Row 20 




Row 21-22 


Row 23 





Reaction label. 

Separator. 

Projectile charge. 

Target charge. 

Charge of exit particle; it refers to the channel other than the entrance and the 7-ray 

channel; = if only two channels are open and radiative capture is the only possible 

reaction. 

Projectile mass in atomic mass units. 

Target mass in atomic mass units. 

Mass of exit particle in atomic mass units. 

Projectile spin. 

Target spin. 

Spin of exit particle. 

Separation energy of incident particle in keV. 

Separation energy of exit particle in keV. 

Radius parameter in fm, used for calculating penetration factor. 

Label of 7-ray channel; channel 1 refers to the incident particle, channel 2 to the emitted 

quantum, and channel 3 to the spectator quantum. 

Separator. 

Minimum energy cutoff (in keV) for numerical integration of rates. 

Number of random samples. 

Flag for temperature output; = outputs results at all temperatures. 

Separator. 

Comments. 

Input for nonresonant contribution; S, S, S are the parameters S(0), S'(0), S"(0) of the 

astrophysical S-factor in units of keV b, b, b/keV, respectively; f racErr is the fractional 

uncertainty, yM/[x]/^E[x] , of the effective S-factor; Cutoff Energy labels the energy E cuto g 

(in keV) at which the S-factor is cut off at higher energies; it is related to the cutoff 

temperature by the expression 

3/2 
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1QQ9F°/^ / / 72 72 MqMi 
iy - y2jfc/ cutoff/ V Z Z 1 Mn+M 



Row 24: 
Row 25: 
Row 26-29: 



9,cutoff — ^--^ -^cutoff/ y ^0^1 M0+M1 

where -E cu toff is m MeV and all other quantities have the same meaning as in Sec. (2) 

Input for a second nonresonant contribution, if needed. 

Separator. 

Comments. 
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Row 30: Input for resonance contribution, one input row for each resonance; Ecm, DEcm: resonance 

energy and la uncertainty; wg, Dwg: resonance strength, wy, and associated uncertainty; 
Jr: resonance spin; G1,DG1,L1: incident particle partial width, uncertainty, orbital angu- 
lar momentum quantum number; for a subthreshold resonance {E r < 0) the dimensionless 
reduced width is listed instead of the entrance channel partial width; G2,DG2,L2: partial 
width, uncertainty, angular momentum quantum number (or multipolarity for 7-rays) of 
emitted quantum; G3,DG3,L3: partial width, uncertainty, angular momentum quantum 
number (or multipolarity for 7-rays) of spectator quantum; Exf : excitation energy of level 
in residual nucleus that is populated in primary transition; Int: = for analytical rate 
calculation; = 1 for numerical rate calculation; for E r < the rate contribution is always 
computed numerically; when the resonance strength is entered the rate contribution is 
always computed analytically, regardless of the flag value; Ecm, DEcm, Exf are in units of 
keV, while resonance strengths and partial widths are in eV. 

Row 31: Input for second resonance; in this example, the rate contribution is calculated from partial 

widths and the rate is chosen to be integrated numerically. 

Row 32: Input for third resonance; in this example, the rate contribution is calculated from the 

resonance strength and the rate is necessarily computed analytically. 

Row 33: Separator. 

Row 34-37: Comments. 

Row 38: Input for resonance contribution when only a partial (or reduced) width upper limit is 

available for at least one reaction channel; the number of upper limit channels must be 
less than the number of open channels; the meaning of the input quantities is the same as 
for row 30, except that (i) resonance strengths are not allowed as input, and (ii) the mean 
value, PT, for the Porter Thomas distribution of dimensionless reduced widths is entered 
for each upper limit channel; upper limit channels are identified by a non-zero value for 
the partial (or reduced) width and a zero value for the corresponding uncertainty; in 
this example, an upper limit for the dimensionless reduced proton width is entered (since 
E r < 0). 

Row 39: Input for second resonance; in this example, an upper limit for the dimensionless reduced 

a-particle width is entered. 

Row 40: Separator. 

Row 41-43: Comments. 

Row 44: Flag for interference between resonances of same spin and parity; +: positive interfer- 

ence; -: negative interference; +- : unknown interference sign (a binary probability density 
function is then used for the random sampling. 

Row 45-46: Input for two interfering resonances; the meaning of the input quantities is the same as 
for row 38, except that the rate contribution is always computed numerically. 

Row 47: Separator. 

Row 48-49: Comments. 
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(iii) Rate probability density figure: This plot consists of six panels showing the Monte Carlo reaction 
rate probability density function (in red) at sample temperatures of T = 0.03, 0.06, 0.1, 0.3, 0.6, and 
1.0 GK. The lognormal approximation to the rate is shown in black. Note that it does not represent a 
least-squares fits to the red histogram, but the lognormal parameters \x and a are computed directly 
from the expectation value and the variance of the logarithm of the rate distribution, according 
to Eq. For a re action, each panel disp lays the temperature, T, and the Anderson-Darling test 
statistic, A-D (see lLongland et al.ll2010bl ). 



(iv) Comparison to literature figure: This plot displays reaction rate ratios as a function of tem- 
perature in two panels, using the following notation: 



NA{o-v) Cyhigh 

NA(o-v) c>med 

Na{o-v) c j ow 

NA(o-v) Pyhi g h 

NA(o-v) p>med 

Na(o-v) P)1ow 



Present high Monte Carlo reaction rate 
Present median Monte Carlo reaction rate 
Present low Monte Carlo reaction rate 



Previous low Monte Carlo reaction rate (Iliadis et al. 



Previous high Monte Carlo reaction rate (llliadis et al.ll2010q) 
Previous median Monte Carlo reaction r ate (llliadis et al.ll2010a ) 



2010 







Top graph: 



upper solid 
lower solid 
upper dashed line 
lower dashed line 



N A {o-v) cMgh /N A {o-v) C:Tned 

NA{o-v) c j ow /N A {o-v) Ctmed 

NA(o-v) pMgh /N A (o-v) P:med 

NA(o-v} p j ow /N A (o-v}^ med 



Bottom graph: 



upper thin solid line NA{o-v) cdl ig h /N A {o-v) p , med 
thick solid line NA{o-v)^ med /N A {o-v) p ^ rned 

lower thin solid line NA(o-v) c j ow /N A (o-v) p , med 
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Table 8:: Total thermonuclear reaction rates for 18 0(p,7) 19 F. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


0.010 


2.92xl0" 24 


4.97xl0" 24 


8.79xl0~ 24 


-5.364xl0+ 01 


5.49xl0^ 01 


7.33x10+°° 


0.011 


2.58xl0~ 23 


4.22xl0" 23 


7.34xl0~ 23 


-5.149xl0+ 01 


5.18xl0~ 01 


1.12xl0+ 01 


0.012 


1.72xl0~ 22 


2.68xl0~ 22 


4.47xl0" 22 


-4.964xl0+ 01 


4.78xl0~ 01 


1.63xl0+ 01 


0.013 


9.30xl0~ 22 


1.37xl0~ 21 


2.18xl0" 21 


-4.800 xl0+ 01 


4.28xl0~ 01 


2.33xl0+ 01 


0.014 


4.22xl0~ 21 


5.97xl0~ 21 


8.99xl0" 21 


-4.653 xl0+ 01 


3.83xl0- 01 


2.60xl0+ 01 


0.015 


1.71xl0~ 2 ° 


2.32xl0~ 2 ° 


3.31xl0~ 2 ° 


-4.519xl0+ 01 


3.40X10" 01 


2.10xl0+ 01 


0.016 


6.18xl0~ 2 ° 


8.15xl0~ 2 ° 


l.llxlO" 19 


-4.393 xl0+ 01 


3.01xl0~ 01 


1.47xl0+ 01 


0.018 


6.28xl0~ 19 


8.03xl0~ 19 


1.04xl0~ 18 


-4.166xl0+ 01 


2.54xl0~ 01 


6.81x10+°° 


0.020 


4.82xl0~ 18 


6.04xl0~ 18 


7.72xl0" 18 


-3.964xl0+ 01 


2.40X10" 01 


6.73x10+°° 


0.025 


3.23xl0" 16 


4.04xl0~ 16 


5.17xl0~ 16 


-3.543 xl0+ 01 


2.44X10" 01 


1.47xl0+ 01 


0.030 


8.73xl0~ 15 


1.09xl0~ 14 


1.40xl0" 14 


-3.212xl0+ 01 


3.29xl0~ 01 


2.41xl0+° 2 


0.040 


1.21xl0~ 12 


1.55xl0~ 12 


2.34xl0~ 12 


-2.701 xl0+ 01 


7.52xl0~ 01 


1.06xl0+° 3 


0.050 


9.35xl0- n 


1.23xl0~ 10 


2.06xl0~ 10 


-2.259 xl0+ 01 


8.58xl0- 01 


1.07xl0+° 3 


0.060 


9.55xl0~° 9 


1.27xl0~° 8 


1.82xl0~° 8 


-1.809xl0+ 01 


5.78X10" 01 


6.00xl0+° 2 


0.070 


3.57xl0-° 7 


4.76X10" 07 


6.50xl0~° 7 


-1.452xl0+ 01 


4.14X10" 01 


2.16xl0+° 2 


0.080 


5.50xl0~° 6 


7.30xl0~° 6 


9.89xl0~° 6 


-1.181xl0+ 01 


3.46X10" 01 


5.98xl0+ 01 


0.090 


4.55xl0~° 5 


6.05xl0"° 5 


8.14xl0~° 5 


-9.706x10+°° 


3.18xl0~ 01 


1.74xl0+ 01 


0.100 


2.44xl0-° 4 


3.24xl0~° 4 


4.35 xlO" 04 


-8.031x10+°° 


3.05X10" 01 


6.39x10+°° 


0.110 


9.50xl0~° 4 


1.26xl0"° 3 


1.69xl0~° 3 


-6.673x10+°° 


2.99xl0~ 01 


2.81x10+°° 


0.120 


2.92xl0~° 3 


3.88xl0~° 3 


5.19xl0"° 3 


-5.551x10+°° 


2.96xl0~ 01 


1.39x10+°° 


0.130 


7.48 xl0~° 3 


9.95xl0"° 3 


1.33xl0~° 2 


-4.610x10+°° 


2.94xl0~ 01 


8.23xl0~ 01 


0.140 


1.66xl0~° 2 


2.21xl0~° 2 


2.96 xl0"° 2 


-3.812x10+°° 


2.93X10" 01 


5.51X10" 01 


0.150 


3.30xl0~° 2 


4.39xl0~° 2 


5.87xl0~° 2 


-3.127x10+°° 


2.92xl0~ 01 


4.28xl0~ 01 


0.160 


5.97xl0~° 2 


7.94xl0~° 2 


1.06xl0~ 01 


-2.533x10+°° 


2.91xl0~ 01 


3.57xl0~ 01 


0.180 


1.58xl0~ 01 


2.10xl0~ 01 


2.81 xl0~ 01 


-1.560x10+°° 


2.90xl0~ 01 


3.07xl0- 01 


0.200 


3.39X10" 01 


4.51xl0~ 01 


6.03xl0~ 01 


-7.970xl0- Q1 


2.90xl0~ 01 


2.95xl0~ 01 


0.250 


1.27x10+°° 


1.69x10+°° 


2.27x10+°° 


5.265xl0~ 01 


2.89xl0~ 01 


2.87xl0~ 01 


0.300 


2.93x10+°° 


3.90x10+°° 


5.21x10+°° 


1.360x10+°° 


2.89xl0- 01 


2.80X10" 01 


0.350 


5.15x10+°° 


6.85x10+°° 


9.13x10+°° 


1.923x10+°° 


2.87X10" 01 


2.71X10" 01 


0.400 


7.70x10+°° 


1.02xl0+ 01 


1.36xl0+ 01 


2.322x10+°° 


2.85xl0~ 01 


2.69xl0~ 01 


0.450 


1.04xl0+ 01 


1.37xl0+ 01 


1.82xl0+ 01 


2.616x10+°° 


2.83xl0~ 01 


2.73xl0~ 01 


0.500 


1.30xl0+ 01 


1.71 xl0+ 01 


2.27xl0+ 01 


2.842x10+°° 


2.79xl0~ 01 


2.91xl0~ 01 


0.600 


1.84xl0+ 01 


2.39xl0+ 01 


3.12xl0+ 01 


3.175x10+°° 


2.66xl0~ 01 


3.86xl0" 01 


0.700 


2.46 xl0+ 01 


3.12xl0+ 01 


3.99 xl0+ 01 


3.445x10+°° 


2.43xl0~ 01 


7.54xl0~ 01 


0.800 


3.36xl0+ 01 


4.14xl0+ 01 


5.13xl0+ 01 


3.727x10+°° 


2.14xl0~ 01 


1.30x10+°° 


0.900 


4.76xl0+ 01 


5.71 xl0+ 01 


6.94xl0+ 01 


4.050x10+°° 


1.91xl0~ 01 


3.36x10+°° 


1.000 


6.92 xl0+ 01 


8.17xl0+ 01 


9.82 xl0+ 01 


4.412x10+°° 


1.79xl0" 01 


1.13xl0+ 01 
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Table 8 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


1.250 


1.72xl0+ 02 


2.00xl0+ 02 


2.38xl0+° 2 


5.312x10+°° 


1.73xl0~ 01 


3.37xl0+ 01 


1.500 


3.63xl0 +02 


4.21 xl0+ 02 


4.98 xl0+ 02 


6.053x10+°° 


1.68xl0" 01 


2.35xl0+ 01 


1.750 


6.40 xl0 +02 


7.43xl0 +02 


8.73 xl0 +02 


6.619x10+°° 


1.63xl0~ 01 


1.33xl0+ 01 


2.000 


9.92 xl0+ 02 


1.15xl0+ 03 


1.34xl0+ 03 


7.053x10+°° 


1.58xl0~ 01 


8.27x10+°° 


2.500 


1.84xl0 +03 


2.13xl0 +03 


2.49 xl0 +03 


7.668x10+°° 


1.52X10" 01 


4.13x10+°° 


3.000 


2.80xl0 +03 


3.23xl0+ 03 


3.77xl0 +03 


8.086x10+°° 


1.53xl0~ 01 


7.12x10+°° 


3.500 


3.78xl0 +03 


4.37xl0 +03 


5.13xl0+ 03 


8.389x10+°° 


1.55xl0~ 01 


1.05xl0+ 01 


4.000 


4.76xl0 +03 


5.51 xl0 +03 


6.51 xl0 +03 


8.625x10+°° 


1.63xl0~ 01 


2.03xl0+ 01 


5.000 


6.60xl0 +03 


7.73xl0+° 3 


9.35 xl0 +03 


8.968x10+°° 


1.80xl0- 01 


2.77xl0+ 01 


6.000 


(8.73xl0+ 03 ) 


(1.06xl0+ 04 ) 


(1.28xl0+ 04 ) 


(9.265x10+°°) 


(1.90xl0- 01 ) 




7.000 


(1.17xl0+ 04 ) 


(1.41 xl0+ 04 ) 


(1.71 xl0+ 04 ) 


(9.555x10+°°) 


(1.90xl0~ 01 ) 




8.000 


(1.45xl0+ 04 ) 


(1.76xl0+ 04 ) 


(2.13xl0+ 04 ) 


(9.774x10+°°) 


(1.90xl0^ 01 ) 




9.000 


(1.72xl0+ 04 ) 


(2.08 xl0+ 04 ) 


(2.51 xl0+ 04 ) 


(9.942x10+°°) 


(1.90X10" 01 ) 




10.000 


(2.03xl0+ 04 ) 


(2.46 xl0+ 04 ) 


(2.97xl0+ 04 ) 


(1.011xl0+ 01 ) 


(1.90xl0~ 01 ) 





C.l. 18 Q(p, 7 ) 19 F 



Comments: This rate was adopted from lBuckner et al.l (J2012l ). which updated the upper and 
lower limits on the E^ = 95 keV resonance strength. In this work, a new upper limit of luj < 7.8 x 
10~ 9 eV (90% CL) was determined. A lognormal probability density function was constructed from 
the improved strength constraints, and the expectation value and standard deviation of this function 
were inp ut into the reactio n rate calculation. A new Q-value (7993.5994 ± 0.0011 keV) was also 
ad opted ( Wang et al.1120121 ) . Proton energies were calculated from the excitation energies presented 



m 



et al.1 
19951' 



Tillev et al.l (119951) fo r all r esonances except the E|^ = 150.8 ± 0.1 keV reson ance measured 



111 



direct ly by iBecker et al.l (j 19951 ) . The direct captur e expansion was ad opted from iBuckner et al 



(120121) . and it was determin ed us ing the code TEDCA (IKrauss et al.lll992l ) and parameters presented 



Iliadis &: Wiescherl (J2004I ). InlBuckner et al.l (120121'). a zero-scattering potential was used based 
on the argument presented in llliadis <fe Wiescherl (|2004l ) . The S-factor fit agreed with the TEDCA 
calculations up to an energy of E£™ off = 2500 keV. For temperatures T > 5 GK, the rate was 



extrapolated with the BRUSLIB Hauser-Feshbach rates (JGorielv et al.ll2008al ) 



180(p,g)19F 

1 ! Zproj 

8 ! Ztarget 

2 ! Zexitparticle (=0 when only 2 channels open) 
1.0078 ! Aproj 

17.999 ! Atarget 
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4.0026 

0.5 

0.0 

0.0 

7993.60 

4013.80 

1.25 

2 



Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 






1.0 
10000 





Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures; =NT for rate output at selected temperatures 



Non-Resonant Contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

7.06e0 2.98e-3 -5.20e-7 0.5 2500.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 



Ecm 


DEcm 


wg 


Dwg 


Jr 


Gl 


DG1 


LI 


G2 


DG2 


L2 


G3 


DG3 


L3 


Exf 


Int 


20.4 


0.7 








2.5 


2.3e-19 


0.5e-19 


2 


2.3 


1.0 


1 


2.5e3 


1.0e3 


3 


0.0 


1 


90.4 


3.0 


5.3e-8 


9.0e-6 
































0.0 





142.8 


0.1 








0.5 


1.67e-l 


0.12e-l 





0.72 


0.15 


1 


1.23e2 


0.24e2 


1 


0.0 


1 


205.4 


1.0 


5.0e-6 


1.0e-6 
































0.0 





260.7 


2.6 


3.7e-5 


0.5e-5 
































0.0 





316.4 


1.3 








2.5 


1.9e-2 


0.3e-2 


2 


0.78 


0.34 


1 


47.0 


19.0 


3 


0.0 


1 


589.9 


1.7 


1.0e-2 


0.2e-2 
































0.0 





598.3 


1.2 








1.5 


1.4e2 


0.7e2 


1 


0.71 


0.39 


1 


2.0e3 


0.1e3 


2 


0.0 


1 


799.6 


1.6 








0.5 


24.6e3 


1.4e3 





2.5 


0.4 


1 


20. e3 


1.0e3 


1 


0.0 


1 


933.1 


2.8 








1.5 


76.0 


7.0 


1 


0.34 


0.06 


1 


3.5e3 


0.3e3 


2 


0.0 


1 


1106.1 


4.0 


0.29 


0.03 
































0.0 





1173.4 


1.5 








0.5 


. 38e3 


0.03e3 





1.4 


1.0 


1 


5.4e3 


. 38e3 


1 


0.0 


1 


1324.4 


2.1 


0.08 


0.01 
































0.0 





1327.4 


1.2 








0.5 


0.22e3 


0.02e3 





3.4 


1.7 


1 


4.7e3 


0.4e3 


1 


0.0 


1 


1542.8 


2.1 


0.025 


0.005 
































0.0 





1572.4 


3.0 


0.041 


0.010 
































0.0 





1581.4 


4.0 


0.06 


0.01 
































0.0 





1592.4 


3.0 


0.025 


0.004 
































0.0 





1673.9 


1.6 








1.5 


2.0e3 


0.6e3 


2 


1.0 


0.4 


1 


1.4e3 


0.4e3 


1 


0.0 


1 


1826.4 


1.2 


2.8 


0.7 
































0.0 





1880.4 


1.9 


0.13 


0.04 
































0.0 





1893.4 


3.0 








0.5 


ll.e3 


3.0e3 





0.36 


0.20 


1 


18.0e3 


5.4e3 


1 


0.0 


1 



Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 

Note: . . .PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 

! 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
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Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf 

!+- 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Almost the same input data as in Angulo et al. 1999 (NACRE) are adopted here; however, we obtain the resonant 
rate contributions by numerical integration, whenever possible. 

2. Proton width for Er=20.4 keV from Champagne and Pitt 1986 and La Cognata et al. 2008, while total and 
radiative widths are from a private communication quoted in Wiescher et al . 1980. 

3. Er=90.4 keV wg expectation value and variance from Buckner et al . 2012. 

4. Direct capture S-f actor adopted from Buckner et al . 2012. 

5. Above T=5 GK the rate is extrapolated using the MOST Hauser-Feshbach rate. 
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18 0(p,y) 19 F 



T9 = 0.03 
A-D = 240.9 



T9 = 0.06 
A-D = 599.E 



1 e-05 2e-05 



CO 



CO 
O 



T9 = 0.1 
A-D = 6.392 



0.005 0.010 0.015 0.020 



d 


1 1 T9 = 0.6 


£ 


A-D = 0.3858 


o - 


ilIUI 


o 


||l 


o 
o - 


,fit 


o 


I'll, 


o 
o - 


l 1 


o 


1 "I 




5 

o - 


1 1 




o 


1 ""liiiiii 


o 
o 

o - 


/ Hk, 




Reaction rate (cm 3 mol 1 s 1 ) 
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18 0(p,Y) 19 F 




0.01 



0.10 



1.00 



10.00 



co 
cc 

CD 
CO 

cc 



to 

CD 

cc 



0.01 




0.10 

Temperature (GK) 



1.00 



10.00 
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Table 9:: Total thermonuclear reaction rates for 18 0(p, 



o 



1 15 



N. 



T (GK) Low rate Median rate High rate lognormal \x lognormal a 



A-D 



0.010 


7.18xl0" 21 


8.66xl0~ 21 


1.04xl0" 20 


-4.620 xl0+ 01 


1.87x10" 


-01 


3.07X10" 01 


0.011 


6.56xl0" 2 ° 


8.00xl0~ 20 


9.76xl0"" 2 ° 


-4.397xl0+ 01 


1.97x10" 


-01 


3.15xl0" 01 


0.012 


4.66xl0~ 19 


5.71xl0~ 19 


6.97xl0" 19 


-4.201 xl0+ 01 


2.03x10" 


-01 


5.47xl0~ 01 


0.013 


2.66xl0~ 18 


3.32xl0~ 18 


4.10xl0~ 18 


-4.025 xl0+ 01 


2.16x10" 


-01 


1.80xl0~ 01 


0.014 


1.31xl0~ 17 


1.63xl0~ 17 


2.04xl0~ 17 


-3.865 xl0+ 01 


2.26x10" 


-01 


1.09x10+°° 


0.015 


5.60xl0~ 17 


6.98xl0~ 17 


8.89xl0~ 17 


-3.719xl0+ 01 


2.38x10" 


-01 


2.59x10+°° 


0.016 


2.13xl0~ 16 


2.68xl0~ 16 


3.44xl0" 16 


-3.585 xl0+ 01 


2.44x10" 


-01 


1.76x10+°° 


0.018 


2.37xl0~ 15 


3.08xl0~ 15 


3.95xl0" 15 


-3.342 xl0+ 01 


2.55x10" 


-01 


4.88xl0~ 01 


0.020 


1.96xl0~ 14 


2.52xl0~ 14 


3.34xl0" 14 


-3.130xl0+ 01 


2.70x10" 


-01 


1.51x10+°° 


0.025 


1.39xl0~ 12 


1.83xl0" 12 


2.41xl0"" 12 


-2.702 xl0+ 01 


2.76x10" 


-01 


9.46X10" 01 


0.030 


3.69X10" 11 


4.84xl0~ n 


6.40xl0~ n 


-2.375 xl0+ 01 


2.76x10" 


-01 


8.57xl0- 01 


0.040 


4.51xl0~° 9 


5.86xl0~° 9 


7.67xl0"° 9 


-1.895xl0+ 01 


2.67x10" 


-01 


9.51xl0~ 01 


0.050 


1.52xl0~° 7 


1.92xl0~° 7 


2.47xl0""° 7 


-1.546xl0+ 01 


2.47x10" 


-01 


1.03x10+°° 


0.060 


3.91 xl0~° 6 


4.48xl0~° 6 


5.30xl0""° 6 


-1.230xl0+ 01 


1.54x10" 


-01 


1.16xl0+ 01 


0.070 


9.27xl0~° 5 


1.00 xlO" 04 


1.08xl0""° 4 


-9.208 xl0+ 00 


8.06x10" 


-02 


2.33x10+°° 


0.080 


1.25xl0~° 3 


1.34xl0~° 3 


1.44xl0~° 3 


-6.612x10+°° 


7.07x10" 


-02 


5.24X10" 01 


0.090 


l.OOxlO" 02 


1.07xl0~° 2 


1.15xl0""° 2 


-4.535x10+°° 


7.01x10" 


-02 


2.53xl0~ 01 


0.100 


5.29xl0~° 2 


5.67xl0~° 2 


6.09xl0"° 2 


-2.870x10+°° 


7.09x10" 


-02 


2.24xl0- 01 


0.110 


2.05xl0~ 01 


2.20xl0~ 01 


2.36xl0~ 01 


-1.515x10+°° 


7.11x10" 


-02 


3.97X10" 01 


0.120 


6.29xl0~ 01 


6.73X10" 01 


7.24xl0~ 01 


-3.944xl0~ 01 


7.12x10" 


-02 


4.06xl0" 01 


0.130 


1.61 xl0+ 00 


1.73xl0+ 00 


1.86x10+°° 


5.462xl0~ 01 


7.13x10" 


-02 


2.71xl0" 01 


0.140 


3.58x10+°° 


3.84x10+°° 


4.13x10+°° 


1.345x10+°° 


7.11x10" 


-02 


3.75xl0- 01 


0.150 


7.11x10+°° 


7.61x10+°° 


8.19x10+°° 


2.031x10+°° 


7.12x10" 


-02 


3.61xl0~ 01 


0.160 


1.29xl0+ 01 


1.38xl0+ 01 


1.48xl0+ 01 


2.626x10+°° 


7.12x10" 


-02 


3.26xl0~ 01 


0.180 


3.43xl0+ 01 


3.67xl0+ 01 


3.95xl0+ 01 


3.605x10+°° 


7.10x10" 


-02 


4.34xl0~ 01 


0.200 


7.42 xl0+ 01 


7.94xl0+ 01 


8.53xl0+ 01 


4.375x10+°° 


7.05x10" 


-02 


3.73X10" 01 


0.250 


2.88xl0+° 2 


3.08xl0+° 2 


3.31xl0+° 2 


5.732x10+°° 


6.99x10" 


-02 


5.87xl0~ 01 


0.300 


7.09xl0+° 2 


7.62xl0+° 2 


8.19xl0+° 2 


6.636x10+°° 


7.24x10" 


-02 


2.06X10" 01 


0.350 


1.41xl0+° 3 


1.53xl0+° 3 


1.67xl0+° 3 


7.336x10+°° 


8.48x10" 


-02 


1.05x10+°° 


0.400 


2.58xl0+° 3 


2.86xl0+° 3 


3.23xl0+° 3 


7.968x10+°° 


1.13x10" 


-01 


7.92x10+°° 


0.450 


4.73xl0+° 3 


5.47xl0+° 3 


6.36xl0+° 3 


8.612x10+°° 


1.48x10" 


-01 


4.23x10+°° 


0.500 


8.94xl0+° 3 


1.06xl0+° 4 


1.27xl0+° 4 


9.272x10+°° 


1.77x10" 


-01 


4.04x10+°° 


0.600 


3.08xl0+° 4 


3.73xl0+° 4 


4.61xl0+° 4 


1.054xl0+ 01 


2.06x10" 


-01 


1.94x10+°° 


0.700 


9.01 xl0+ 04 


l.llxl0+° 5 


1.37xl0+° 5 


1.162xl0+ 01 


2.14x10" 


-01 


1.29x10+°° 


0.800 


2.17xl0+ 05 


2.71xl0+° 5 


3.36xl0+° 5 


1.251xl0+ 01 


2.17x10" 


-01 


5.57xl0" 01 


0.900 


4.53xl0+ 05 


5.56xl0+° 5 


6.91xl0+° 5 


1.323xl0+ 01 


2.10x10" 


-01 


1.95x10+°° 


1.000 


8.27xl0+ 05 


1.02xl0+° 6 


1.25xl0+° 6 


1.383xl0+ 01 


2.04x10" 


-01 


7.73xl0~ 01 
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Table 9 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a 


A-D 


1.250 


2.56xl0+ 06 


3.08 xl0+ 06 


3.71 xl0+ 06 


1.494xl0+ 01 


1.89xl0" 01 


2.33x10+°° 


1.500 


5.49xl0+° 6 


6.50xl0 +06 


7.81 xl0+ 06 


1.569xl0+ 01 


1.77xl0~ 01 


2.42x10+°° 


1.750 


9.44xl0+ 06 


l.llxl0 +07 


1.31 xl0 +07 


1.623xl0+ 01 


1.65xl0~ 01 


2.24x10+°° 


2.000 


1.40xl0+ 07 


1.63xl0 +07 


1.92xl0+° 7 


1.661xl0+ 01 


1.61xl0~ 01 


1.33x10+°° 


2.500 


2.38xl0 +07 


2.75xl0 +07 


3.23xl0 +07 


1.714xl0+ 01 


1.55xl0" 01 


3.67x10+°° 


3.000 


3.30xl0+ 07 


3.79xl0 +07 


4.42 xl0 +07 


1.746xl0+ 01 


1.48xl0~ 01 


2.20x10+°° 


3.500 


4.09xl0+° 7 


4.70xl0+° 7 


5.44xl0+° 7 


1.767xl0+ 01 


1.45xl0~ 01 


1.46x10+°° 


4.000 


4.72 xl0+ 07 


5.43xl0 +07 


6.28xl0 +07 


1.781xl0+ 01 


1.44xl0~ 01 


9.12xl0" 01 


5.000 


5.68xl0 +07 


6.49 xl0 +07 


7.47xl0+° 7 


1.799xl0+ 01 


1.37xl0~ 01 


9.16X10" 01 


6.000 


6.40xl0+ 07 


7.26xl0+° 7 


8.30xl0+° 7 


1.810xl0+ 01 


1.32xl0~ 01 


1.34x10+°° 


7.000 


6.98xl0+ 07 


7.86 xl0 +07 


8.97xl0 +07 


1.819xl0+ 01 


1.25xl0~ 01 


3.57x10+°° 


8.000 


7.62 xl0+ 07 


8.52xl0 +07 


9.56xl0 +07 


1.826xl0+ 01 


1.15xl0~ 01 


2.65x10+°° 


9.000 


8.29xl0+ 07 


9.17xl0 +07 


1.02 xlO+ 08 


1.834xl0+ 01 


1.06xl0~ 01 


3.06x10+°° 


10.000 


8.98xl0+ 07 


9.89 xl0+ 07 


1.09xl0+ 08 


1.841 xl0+ 01 


9.83xl0" 02 


2.30x10+°° 



C.2. i8 0(p,a) 15 N 



Comments: In the previous input file for 18 0(p, a) 15 N (see llliadis et al.ll2010bl ). only two inter- 
fering resonances were taken into account, and the interference sign was assumed to be unknown. 
Here, we take interference between three 1/2+ states into account at Er = 143, 609, and 807 keV. 
The interference sign between the 609 and 807 keV resonances is known to be (-) (i.e., destruc- 
tive outside the resonances and constructive between) . This was determined by calcula ting the 
cross section for each sign and comparing with the data from lLorenz-Wirzba et al.l (11979!) . Ener - 
gies and partial w i dths w ere adopted from the weighted average of Yagil (1962T ^ Mak et al.l ( 19781 ); 
La Cognata et al.l ( 2010al ). using the method discussed in lWietfeldt Sz Greene! ( 20 111 ) to expand un- 
certainties by adding an unknown systematic uncertainty in quadrature. For the 143 and 609 keV 
resonances, the (-) interference sig n is known acco rding to the yield curve shown in iBecker et al 



(|l995l ); lLorenz-Wirzba et al.l (J19791 ) (strength from lBecker et al.lll995l ). Because the RatesMC code 
can only include interference from two resonances, we ran the code twice with two input files. 
Input file #1: the 143 and 609 keV resonances interfering, with the 807 resonance in the resonant 
contribution section (shown in the input file below). Input file ^2: the 609 and 807 keV resonances 
interfering, with the 143 keV resonance in the resonant contribution section. Between temperatures 
of 0.07 and 0.16 GK, the rates deviate from one another by less than a percent. Therefore, the 
rate tabulated here is from input file #1 for T < 0.06 GK, while above this threshold, the rate is 
from input file #2. Only one input file is shown for clarity (input file #1). To construct input file 
#2, move the "!" comment from the 143 keV resonance to the 807 keV resonance in the resonant 
contribution section and from the 807 keV resonance to the 143 keV resonance in the interference 
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section. 

Additional changes from the previous input file include moving the 90-keV resonance from 
the upper limit sec tion to the resonant c ontrib ution section with a strength of 1.6 ± 0.5 eV di- 
rectly measured by lLorenz-Wirzba et al.1 (|1979l ). In addition, re sonance ene r gies le s s than 2 MeV 
were a d opted from the 1 8 Q(p, 7 ) 19 F reaction i nput f ile (see IC.1|) . ISellin et al.1 (|1969l ) ; lOrihara et al. 
(j 19731 ) ; lAlmanza et al.l (j 19751 ) ; iMurillo et al.l (|1979l ) provided values for the highe r resonance en- 
ergies and partial widths, and Q- values and masses are from IWang et al.1 (120121 ). Our adopted 
strengths for the 20 and 90 keV resonances agree with lLa Cognata et al.l (|2010bl ). 
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17.999 
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Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 






1.0 
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Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures; =NT for rate output at selected temperatures 



Non-Resonant Contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 



i i 



icm 


DEcm 


wg 


Dwg 


Jr 


Gl 


DG1 


LI 


G2 


DG2 


L2 


G3 


DG3 


L3 


Exf 


Int 


20.4 
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1172.2 1.5 

1326.2 1.2 

1672.7 1.6 

1825.2 1.2 

1892.0 3.0 

2167.0 3.0 

2237.0 3.0 

2259.0 3.0 

2313.0 4.0 

2501.5 1.4 

2619.5 1.7 

2768 .5 2.6 

2864 .9 2.0 

2980.2 2.6 

3291.0 7.0 

3355.0 25.0 

3455.0 3.5 

3507.0 5.0 

3545 .0 7.0 

3608.0 12.0 

3658 .0 4.0 

4045 .0 20 . 

4141.0 8.0 

4227.0 12. 

4527.0 7.0 

4582.0 10.0 

4585.0 25.0 

4785.0 10.0 

4865 .0 30 . 

4945.0 25.0 

4985.0 50.0 

5095.0 75.0 

5322 .0 8.0 

5365.0 25.0 

5737.0 11.0 

6045 .0 20 . 

6105.0 21.0 

6335.0 20.0 

6705.0 20.0 

6745 .0 50 . 

6785.0 20.0 

6925.0 30.0 

7365.0 20.0 

7405.0 30.0 

7775.0 21.0 

8205.0 40.0 

8235.0 30.0 

8285.0 20.0 

9055.0 40.0 

9165.0 40.0 
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Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 
Note: ...PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 
Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 
204.2 1.0 2.5 7.7e-4 2.0e-4 2 0.0 0.8e3 0.0 3 0.010 0.0 0.0 0.0 0.0 1 

Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf 



142.8 0.1 


0.5 


1.67e-l 


0.12e-l 
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1.23e2 


0.24e2 
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0.15 
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0.0 
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! 807 . 3.0 


0.5 


2 . 46e4 


0.14e4 
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2.2e4 


4.0e3 
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0.0 
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Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Up to Er=2 MeV, input data are taken from the same source as for the 180(p,g)19F reaction (see previous 
table), except for comments below; for higher energies the partial widths are adopted from Sellin et al . 1969, 
Orihara et al . 1973, Almanza et al . 1975 and Murillo et al . 1979. 

2. Q-values and masses from Audi and Meng 2011. 

3. For 90 keV resonance, directly measured strength adopted from Lorenz-Wirzba et al . (1979). 

4. Resonance energy for 143 keV from Becker et al. 1995. 

5. For 20 and 90 keV resonances, strengths agree with La Cognata et al . (2010b). 

6. For energies and partial widths of broad 609 keV and 807 keV resonances, we adopted the weighted average from 
Yagi 1962, Mak et al . 1978, and La Cognata et al . 2010a, using the "additive variable method"; see Sec. 8.2.; 
note that the results of Lorenz-Wirzba et al . 1979 are explicitly included in the values reported by La Cognata 
et al. 2010a. 

7. Interferences between 143, 609, and 807 keV resonances: we assumed that the two-level interference is 
destructive (constructive) outside (between) (i) 143-609 keV (see Becker et al. 1995 and Lorenz-Wirzba et al. 
1979); (ii) 906-807 keV (see Lorenz-Wirzba et al . 1979 and La Cognata et al . 2010a). 
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8 0(p, <x) 15 N 
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Table 10:: Total thermonuclear reaction rates for 

22 Ne(p, 7 ) 23 Na. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a 


A-D 


0.010 


4.21xl0~ 25 


6.76xl0" 25 


1.07xl0~ 24 


-5.566 xl0+ 01 


4.68xl0" 01 


4.60xl0" 01 


0.011 


1.60xl0 -23 


2.46xl0~ 23 


3.71 xlO" 23 


-5.206 xl0+ 01 


4.27xl0~ 01 


4.88xl0~ 01 


0.012 


3.25xl0~ 22 


4.85xl0~ 22 


7.12xl0~ 22 


-4.908 xl0+ 01 


4.00xl0~ 01 


4.80xl0~ 01 


0.013 


4.07xl0~ 21 


6.00xl0~ 21 


8.64xl0" 21 


-4.657xl0+ 01 


3.83xl0- 01 


3.83xl0- 01 


0.014 


3.52xl0~ 2 ° 


5.10xl0~ 2 ° 


7.35 xl0~ 20 


-4.442 xl0+ 01 


3.73xl0~ 01 


2.28xl0~ 01 


0.015 


2.25xl0~ 19 


3.24xl0~ 19 


4.70xl0~ 19 


-4.257xl0+ 01 


3.70xl0~ 01 


2.06xl0~ 01 


0.016 


1.13xl0~ 18 


1.63xl0~ 18 


2.37xl0~ 18 


-4.096 xl0+ 01 


3.69xl0~ 01 


2.80xl0~ 01 


0.018 


1.64xl0~ 17 


2.37xl0~ 17 


3.44xl0~ 17 


-3.828 xl0+ 01 


3.76xl0~ 01 


3.33xl0~ 01 


0.020 


1.35X10" 16 


1.98xl0~ 16 


2.91xl0~ 16 


-3.616xl0+ 01 


3.87xl0- 01 


4.12xl0~ 01 


0.025 


5.71xl0~ 15 


8.68xl0~ 15 


1.30xl0" 14 


-3.239 xl0+ 01 


4.19xl0~ 01 


5.82xl0~ 01 


0.030 


6.55xl0~ 14 


1.02xl0~ 13 


1.57xl0~ 13 


-2.992 xl0+ 01 


4.47xl0~ 01 


6.44xl0~ 01 


0.040 


1.25xl0~ 12 


2.04xl0~ 12 


3.27xl0~ 12 


-2.693 xl0+ 01 


4.87xl0~ 01 


6.52xl0~ 01 


0.050 


6.93xl0" 12 


1.16xl0 -11 


1.89xl0 _u 


-2.519xl0+ 01 


5.08xl0~ 01 


5.88xl0~ 01 


0.060 


2.27X10" 11 


3.74xl0~ n 


6.08X10" 11 


-2.401 xl0+ 01 


4.97xl0~ 01 


3.37xl0~ 01 


0.070 


7.06X10" 11 


1.05xl0~ 10 


1.59xl0~ 10 


-2.297xl0+ 01 


4.13xl0~ 01 


1.83x10+°° 


0.080 


2.85xl0~ 10 


3.84xl0~ 10 


5.18xl0~ 10 


-2.168xl0+ 01 


3.00xl0~ 01 


2.97xl0~ 01 


0.090 


1.32xl0~° 9 


1.89xl0-° 9 


2.67xl0"° 9 


-2.010xl0+ 01 


3.35xl0~ 01 


9.47x10+°° 


0.100 


5.52xl0~° 9 


9.22xl0~ 09 


1.56xl0~° 8 


-1.850xl0+ 01 


4.66xl0~ 01 


3.13xl0+ 01 


0.110 


2.08xl0~° 8 


4.02xl0~ 08 


7.67xl0"° 8 


-1.704xl0+ 01 


5.73xl0~ 01 


4.61 xl0+ 01 


0.120 


6.51xl0~° 8 


1.45xl0~ 07 


3.00xl0~° 7 


-1.578xl0+ 01 


6.62xl0~ 01 


5.75xl0+ 01 


0.130 


1.80xl0~° 7 


4.48xl0~ 07 


9.60xl0~° 7 


-1.468xl0+ 01 


7.23xl0~ 01 


6.86xl0+ 01 


0.140 


4.42xl0~° 7 


1.18xl0~ 06 


2.61 xlO" 06 


-1.373xl0+ 01 


7.65xl0~ 01 


6.94xl0+ 01 


0.150 


9.89xl0~° 7 


2.76xl0~ 06 


6.27xl0"° 6 


-1.288xl0+ 01 


7.92xl0~ 01 


7.35 xl0+ 01 


0.160 


2.01xl0~° 6 


5.90xl0~ 06 


1.33xl0"° 5 


-1.214xl0+ 01 


8.08xl0~ 01 


7.61 xl0+ 01 


0.180 


7.34xl0~° 6 


2.09xl0~ 05 


4.73 xlO" 05 


-1.086xl0+ 01 


7.91 xl0~ 01 


7.67xl0+ 01 


0.200 


2.68xl0~° 5 


6.30xl0~ 05 


1.33xl0"° 4 


-9.715x10+°° 


6.81xl0~ 01 


7.40xl0+ 01 


0.250 


9.43 xl0~° 4 


1.18X10" 03 


1.57xl0~° 3 


-6.722x10+°° 


2.30xl0~ 01 


2.81 xl0+ 01 


0.300 


1.77xl0-° 2 


1.98xl0-° 2 


2.23xl0"° 2 


-3.919x10+°° 


1.17X10" 01 


9.06X10" 01 


0.350 


1.56xl0~ 01 


1.73xl0~ 01 


1.92xl0" 01 


-1.751x10+°° 


1.06xl0~ 01 


1.06x10+°° 


0.400 


8.19xl0~ 01 


9.03xl0~ 01 


9.97xl0~ 01 


-1.006xl0~ 01 


9.93xl0~ 02 


7.22xl0~ 01 


0.450 


3.00xl0 +00 


3.28xl0+ 00 


3.60x10+°° 


1.190x10+°° 


9.33xl0~ 02 


4.60xl0~ 01 


0.500 


8.50x10+°° 


9.27x10+°° 


1.01xl0+ 01 


2.228x10+°° 


8.81 xl0~ 02 


3.56xl0~ 01 


0.600 


4.13xl0+ 01 


4.47xl0+ 01 


4.83 xl0+ 01 


3.800x10+°° 


8.04xl0~ 02 


2.68xl0~ 01 


0.700 


1.30xl0+° 2 


1.41xl0+° 2 


1.51xl0+° 2 


4.945x10+°° 


7.61 xl0~ 02 


2.87xl0~ 01 


0.800 


3.14xl0+° 2 


3.38xl0+° 2 


3.64xl0+° 2 


5.824x10+°° 


7.46xl0~ 02 


4.03xl0~ 01 


0.900 


6.32 xl0+ 02 


6.81xl0+° 2 


7.34 xlO+ 02 


6.524x10+°° 


7.51 xlO" 02 


5.80xl0" 01 
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Table 10 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


1.000 


1.12xl0+ 03 


1.21 xl0+ 03 


1.30xl0+° 3 


7.100x10+°° 


7.68xl0" 02 


7.91 xl0~ 01 


1.250 


3.28xl0 +03 


3.54xl0+ 03 


3.86 xl0 +03 


8.176x10+°° 


8.35xl0"° 2 


2.51x10+°° 


1.500 


6.90xl0 +03 


7.49xl0 +03 


8.22 xl0 +03 


8.926x10+°° 


9.01 xl0"° 2 


4.82x10+°° 


1.750 


1.20 xl0+ 04 


1.30xl0+ 04 


1.44xl0+ 04 


9.479x10+°° 


9.45xl0~° 2 


6.75x10+°° 


2.000 


1.82 xl0+ 04 


1.99xl0+ 04 


2.20xl0 +04 


9.901x10+°° 


9.68xl0~° 2 


7.77x10+°° 


2.500 


3.31 xl0 +04 


3.60 xl0+ 04 


3.99 xl0+ 04 


1.050xl0+ 01 


9.67xl0~° 2 


8.73x10+°° 


3.000 


4.93xl0+ 04 


5.35xl0+° 4 


5.90 xl0 +04 


1.089xl0+ 01 


9.37xl0~° 2 


8.78x10+°° 


3.500 


6.51 xl0 +04 


7.05 xlO+° 4 


7.74xl0+ 04 


1.117xl0+ 01 


8.98xl0~° 2 


8.67x10+°° 


4.000 


7.96 xl0 +04 


8.59xl0+° 4 


9.39 xl0 +04 


1.137xl0+ 01 


8.59X10" 02 


8.44x10+°° 


5.000 


(1.07xl0+ 05 ) 


(1.16xl0+ 05 ) 


(1.27xl0+ 05 ) 


(1.166xl0+ 01 ) 


(8.56xl0"° 2 ) 




6.000 


(1.31xl0+° 5 ) 


(1.43xl0+° 5 ) 


(1.56xl0+° 5 ) 


(1.187xl0+ 01 ) 


(8.56xl0~° 2 ) 




7.000 


(1.51xl0+ 05 ) 


(1.65xl0+ 05 ) 


(1.79xl0+ 05 ) 


(1.201xl0+ 01 ) 


(8.56xl0"° 2 ) 




8.000 


(1.68xl0+ 05 ) 


(1.83xl0+ 05 ) 


(2.00xl0+° 5 ) 


(1.212xl0+ 01 ) 


(8.56xl0"° 2 ) 




9.000 


(1.81xl0+ 05 ) 


(1.98xl0+° 5 ) 


(2.15xl0+° 5 ) 


(1.219xl0+ 01 ) 


(8.56xl0"° 2 ) 




10.000 


(1.96xl0+° 5 ) 


(2.13xl0+ 05 ) 


(2.32xl0+ 05 ) 


(1.227xl0+ 01 ) 


(8.56xl0" 02 ) 





C.3. 22 Ne(p, 7 ) 23 Na 



Comments: The input file from llliadis et al.l (|2010bl ) has been updated with a new normal- 
ization of all previously measured resonance strength s, based on an absolut e measurement of the 
E p = 458 keV resonance (ujj = 0.524 ± 0.051 eV) (ILongland et al.ll2010al ). Uncertainties from 
the previous strengths and the uncertainty of the new measurement were added in quadrature. 
Particle partial widths for the two subthreshold resonances at E p = 35 and 151 keV have not been 
renormalized. Values for upper limits remain unchanged. For T > 4.0 GK, rates were extrapolated 
using Hauser-Feshbach calculations. 



22Ne(p,g)23Na 

1 ! Zproj 

10 ! Ztarget 

! Zexitparticle (=0 when only 2 channels open) 

1.0078 ! Aproj 

21.991 ! Atarget 

! Aexitparticle (=0 when only 2 channels open) 

0.5 ! Jpro j 

0.0 ! Jtarget 

0.0 ! Jexitparticle (=0 when only 2 channels open) 

8794.11 ! projectile separation energy (keV) 

0.0 ! exit particle separation energy (=0 when only 2 channels open) 

1.25 ! Radius parameter RO (fm) 
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2 ! Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 



1.0 

5000 





Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures; =NT for rate output at selected temperatures 



Non-resonant contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

6.2el 0.0 0.0 0.4 1500.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 



i i 



Ecm DEcm wg Dwg Jr Gl DG1 LI G2 



35.4 


0.7 








0.5 


3 


150.9 


2.0 








3.5 


2 


417.0 


0.8 


0.076 


0.019 








458.2 


0.8 


0.524 


0.051 








601.9 


0.3 


0.035 


0.012 








610.3 


0.3 


3.3 


0.5 








631.4 


0.4 


0.41 


0.12 








693.3 


0.7 


0.15 


0.05 








813.7 


0.2 


8.2 


3.0 








857.3 


0.5 


2.1 


1.1 








861.1 


1.0 


1.22 


0.59 








879.5 


1.0 


0.93 


0.31 








888.2 


0.3 


0.41 


0.12 








906.3 


1.0 


7.0 


2.4 








937.91 


0.07 


0.47 


0.13 








960.9 


0.5 


2.8 


0.9 








1021.0 


0.4 


1.0 


0.3 








1040.7 


1.0 


2.50 


0.69 








1055.2 


0.5 


2.50 


0.69 








1096.1 


0.6 


1.7 


0.5 








1122.2 


0.6 


0.70 


0.19 








1208.5 


0.6 


1.28 


0.37 








1221.9 


0.4 


12.2 


1.7 








1254.3 


0.6 


0.23 


0.06 








1275.9 


0.6 


3.20 


0.87 








1281.1 


0.5 


1.40 


0.38 








1290.3 


0.2 


0.93 


0.25 








1319.9 


0.5 


0.76 


0.19 








1331.0 


0.5 


1.63 


0.44 








1374 . 7 


0.2 


3.1 


0.9 








1436 . 7 


0.3 


5.2 


1.3 








1448 . 8 


1.4 


1.3 


0.4 








1486 . 6 


0.6 


1.63 


0.44 








1523.1 


0.6 


5.8 


1.8 








1543.7 


0.7 


1.22 


0.31 









DG2 L2 G3 DG3 
.le-15 1.2e-15 
.3e-9 9.2e-10 3 



L3 



Exf 


Int 














2.2 


1.0 


1 











0.0 


1 


0.02 


0.01 


1 











0.0 


1 




















0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 























0.0 






57 



1551.1 


0.7 


5.2 


1.3 
































0.0 





1558.9 


0.7 


3.49 


0.94 
































0.0 





1645.6 


1.0 


7.6 


1.9 
































0.0 





1653.7 


1.2 


2.04 


0.56 
































0.0 





1683.8 


0.7 


2.97 


0.81 
































0.0 





1706.8 


0.7 


3.84 


1.06 
































0.0 





1712.8 


0.7 


0.58 


0.18 
































0.0 





1724.1 


0.7 


2.56 


0.69 
































0.0 





1739.1 


0.7 


0.87 


0.25 
































0.0 





1754.2 


0.9 


6.4 


1.9 
































0.0 





1779.5 


0.8 


1.34 


0.37 
































0.0 





1821.8 


0.8 


4.37 


1.19 
































0.0 






Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 

Note: . . .PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 



Ecm DE< 


:m Jr 


Gl 


DG1 


LI 


PT 


G2 DG2 


L2 PT 




G3 


DG3 


L3 


PT 


Exf 


Int 




27.9 


3.0 


4.5 


5.2e-26 


0.0 


5 


0.0045 0.1 


0.01 


1 

















0.0 





177.9 


2.0 


0.5 


2.6e-6 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 





247.9 


1.0 


3.5 


3.3e-8 


0.0 


4 


0.0045 0.04 


0.02 


1 

















0.0 


1 


277.9 


3.0 


0.5 


2.2e-6 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 





308.9 


3.0 


0.5 


2.2e-6 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 





318.9 


3.0 


0.5 


3.0e-6 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 





352.9 


5.0 


0.5 


6.0e-4 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 





376.9 


3.0 


0.5 


6.0e-4 


0.0 





0.0045 0.1 


0.01 


1 

















0.0 






Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

L2 PT G3 DG3 L3 PT Exf 



Ecm DEcm Jr 

!+- 

0.0 0.0 


Gl 


DG1 


LI PT 


G2 


DG2 


0.0 


0.0 


0.0 





0.0 


0.0 0.0 


0.0 


0.0 


0.0 





0.0 



0.0 0.0 0.0 0.0 
0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Information for Er>400 keV from Endt 1990 (strengths normalized relative to Er=458 keV) . 

2. Er=178, 278-377 keV: s-wave resonances (Jp=l/2+) assumed for upper limit; value of Gg=0.1+-0.01 eV is a guess 
(inconsequential since Gp<<Gg) ; Gp upper limit values calculated from strength upper limits of Goerres et al. 
1982. 

3. Er=28 keV: h-wave resonance (Jp=9/2-) assumed for upper limit; Gg=0.1+-0.01 eV is a guess (inconsequential). 

4. Er=151 keV: contrary to Hale et al . 2001, we adopt C2S=0.0011 (see Hale's Ph.D. thesis). 

5. Er=248 keV: g-wave resonance (Jp=7/2+) assumed for upper limit. 

6. Direct capture S-f actor adopted from Goerres et al . 1983, with uncertainty estimate from Hale et al. 2001. 

7. Levels at Ex=8862, 8894 and 9000 keV (Powers et al . 1971) have been disregarded. 

8. All strengths have been renormalized to the resonance at 458 keV (Longland 2010a), with the exception of 
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upper limits and resonances at 35 and 151 keV. 



2 Ne(p, y) 23 Na 



o 




T9 
A-D 


= 0.06 
= 0.3373 


yiiuj|UL_ 


ii i in 



1.5e-10 2.09-10 



T9 = 0.3 
A-D = 0.9059 



- 










T9 = 0.6 
A-D = 0.2683 

IkiLi, 




Reaction rate (cm 3 mol 1 s 1 ) 
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2 Ne(p,7) 23 Ne 





Temperature (GK) 
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Table 11:: Total thermonuclear reaction rates for 
22 Ne(a, 7 ) 26 Mg. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal /i 


lognormal a 


A-D 


0.010 


1.05xl0" 77 


2.14xl0" 77 


4.52xl0~ 77 


-1.765xl0+° 2 


7.42xl0" 01 


3.20xl0" 01 


0.011 


3.99xl0" 74 


7.28xl0" 74 


1.34xl0" 73 


-1.684xl0 +02 


6.15xl0" 01 


8.94xl0"° 2 


0.012 


3.69xl0" 71 


6.34xl0" 71 


1.07xl0" 7 ° 


-1.617xl0+ 02 


5.34xl0" 01 


1.60xl0" 01 


0.013 


1.15xl0" 68 


1.90xl0" 68 


3.09 xlO" 68 


-1.559xl0+ 02 


4.92xl0" 01 


1.99xl0" 01 


0.014 


1.55xl0" 66 


2.52xl0" 66 


4.04xl0" 66 


-1.511xl0+ 02 


4.80xl0" 01 


1.92xl0" 01 


0.015 


1.06xl0" 64 


1.73xl0" 64 


2.79 xlO" 64 


-1.468xl0+ 02 


4.90xl0" 01 


2.81xl0" 01 


0.016 


4.11 xlO" 63 


6.96xl0" 63 


1.14xl0" 62 


-1.431 xl0+ 02 


5.13xl0" 01 


4.12xl0" 01 


0.018 


1.80xl0" 6 ° 


3.26xl0" 6 ° 


5.63 xlO" 60 


-1.370xl0+ 02 


5.75xl0" 01 


7.20xl0" 01 


0.020 


2.24xl0" 58 


4.34xl0" 58 


8.04xl0" 58 


-1.321xl0+ 02 


6.43xl0" 01 


9.51xl0" 01 


0.025 


1.54xl0" 54 


3.14xl0" 54 


6.30 xlO" 54 


-1.232xl0+ 02 


7.13xl0" 01 


6.38xl0" 01 


0.030 


2.82xl0" 5 ° 


3.35xl0" 49 


1.30 xlO" 48 


-1.121xl0+ 02 


1.87xl0+ 00 


9.43xl0+ 01 


0.040 


1.81 xlO" 42 


2.31 xlO" 41 


8.91 xlO" 41 


-9.413xl0+ 01 


2.14x10+°° 


1.31xl0+° 2 


0.050 


8.51 xlO" 38 


1.08 xlO" 36 


4.17xl0" 36 


-8.338xl0+ 01 


2.15x10+°° 


1.32xl0+° 2 


0.060 


1.05xl0" 34 


1.34xl0 -33 


5.14xl0~ 33 


-7.624xl0+ 01 


2.08x10+°° 


1.25xl0+° 2 


0.070 


1.95xl0" 32 


2.12xl0" 31 


8.04xl0~ 31 


-7.104xl0+ 01 


1.79x10+°° 


7.79xl0+ 01 


0.080 


2.76xl0" 3 ° 


1.14xl0" 29 


3.67xl0~ 29 


-6.679xl0+ 01 


1.33x10+°° 


3.60xl0+ 01 


0.090 


1.76xl0" 28 


6.30xl0" 28 


1.35xl0" 27 


-6.289xl0+ 01 


1.15x10+°° 


1.26xl0+° 2 


0.100 


4.79xl0" 27 


2.28xl0" 26 


6.55xl0~ 26 


-5.931 xlO+ 01 


1.35x10+°° 


9.20xl0+ 01 


0.110 


8.17xl0" 26 


5.95xl0~ 25 


1.86xl0" 24 


-5.616xl0+ 01 


1.55x10+°° 


1.14xl0+° 2 


0.120 


1.11 xlO -24 


9.63xl0" 24 


3.07xl0~ 23 


-5.343xl0+ 01 


1.64x10+°° 


1.34xl0+° 2 


0.130 


1.23xl0" 23 


1.03xl0~ 22 


3.28xl0" 22 


-5.102xl0+ 01 


1.57x10+°° 


1.22xl0+° 2 


0.140 


1.38xl0" 22 


8.23xl0~ 22 


2.50xl0" 21 


-4.883xl0+ 01 


1.36x10+°° 


8.62xl0+ 01 


0.150 


1.53xl0" 21 


5.57xl0" 21 


1.51xl0" 20 


-4.679xl0+ 01 


1.10x10+°° 


5.66xl0+ 01 


0.160 


1.41xl0" 2 ° 


3.79xl0~ 2 ° 


8.10xl0" 20 


-4.484xl0+ 01 


8.63xl0" 01 


4.85 xlO+ 01 


0.180 


8.05xl0" 19 


1.54xl0~ 18 


2.84xl0" 18 


-4.102xl0+ 01 


6.29xl0" 01 


8.06xl0" 01 


0.200 


3.41xl0" 17 


5.43xl0" 17 


9.60xl0" 17 


-3.740xl0+ 01 


5.19xl0" 01 


2.31 xlO+ 01 


0.250 


5.88 xlO" 14 


7.56xl0" 14 


l.OOxlO -13 


-3.019xl0+ 01 


2.78xl0" 01 


1.22xl0+ 01 


0.300 


9.32xl0" 12 


1.13xl0 _11 


1.38xl0" n 


-2.520xl0+ 01 


1.96xl0" 01 


1.18x10+°° 


0.350 


3.46xl0" 10 


4.08xl0" 10 


4.86xl0" 10 


-2.162xl0+ 01 


1.69xl0" 01 


5.63xl0" 01 


0.400 


5.11xl0"° 9 


5.95xl0~° 9 


6.98 xlO" 09 


-1.894xl0+ 01 


1.56xl0" 01 


6.32xl0" 01 


0.450 


4.09xl0"° 8 


4.72 xlO" 08 


5.50xl0"° 8 


-1.686xl0+ 01 


1.47xl0" 01 


6.44xl0" 01 


0.500 


2.13xl0"° 7 


2.44xl0" 07 


2.82xl0" 07 


-1.522xl0+ 01 


1.41xl0" 01 


6.51xl0" 01 


0.600 


2.47xl0"° 6 


2.79xl0" 06 


3.20xl0~ 06 


-1.278xl0 +01 


1.32xl0" 01 


6.04xl0" 01 


0.700 


1.39xl0"° 5 


1.57xl0" 05 


1.78xl0~ 05 


-1.106xl0+ 01 


1.25xl0" 01 


5.52xl0" 01 


0.800 


5.15xl0"° 5 


5.77xl0" 05 


6.51 xlO" 05 


-9.758x10+°° 


1.18xl0" 01 


4.75xl0" 01 


0.900 


1.48 xlO" 04 


1.66xl0~ 04 


1.88X10" 04 


-8.701x10+°° 


1.19xl0" 01 


1.95x10+°° 
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Table 11 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a 


A-D 


1.000 


3.65 xlO" 04 


4.11X10" 04 


4.73 xlO" 04 


-7.788 xl0+ 00 


1.35xl0~ 01 


1.02xl0+ 01 


1.250 


2.33xl0~ 03 


2.77xl0~° 3 


3.43 xlO" 03 


-5.867x10+°° 


2.02xl0~ 01 


2.19xl0+ 01 


1.500 


1.45xl0" 02 ) 


l.TOxlO -02 ) 


2.21xl0~ 02 ) 


(-4.024x10+°°) 


;2.12xl0" 01 ) 




1.750 


7.64xl0~° 2 ) 


;9.45xl0" 02 ) 


1.17xl(T 01 ) 


(-2.360x10+°°) 


;2.12xl0~ 01 ) 




2.000 


3.00xl0" 01 ) 


;3.70xl0~ 01 ) 


4.58xl0~ 01 ) 


(-9.932xl0~ 01 ) 


;2.12xl0" 01 ) 




2.500 


2.55x10+°°) 


;3.15xl0 +00 ) 


3.89xl0 +00 ) 


(1.147x10+°°) 


;2.12xl0" 01 ) 




3.000 


4.24xl0+ 01 ) 


1.53xl0+ 01 ) 


1.89xl0+ 01 ) 


(2.729x10+°°) 


;2.12xl0" 01 ) 




3.500 


4.18xl0+ 01 ) 


;5.17xl0+ 01 ) 


6.39xl0+ 01 ) 


(3.945x10+°°) 


;2.12xl0" 01 ) 




4.000 


1.10xl0+° 2 ) 


1.36xl0+° 2 ) 


1.68xl0+ 02 ) 


(4.913x10+°°) 


;2.12xl0- 01 ) 




5.000 


4.71xl0+° 2 ) 


;5.82xl0+ 02 ) 


7.19xl0+ 02 ) 


(6.366x10+°°) 


;2.12xl0- 01 ) 




6.000 


1.33xl0+° 3 ) 


1.64xl0 +03 ) 


2.03xl0+ 03 ) 


(7.405x10+°°) 


;2.12xl0" 01 ) 




7.000 


2.91xl0+° 3 ) 


;3.59xl0+ 03 ) 


4.44xl0+ 03 ) 


(8.186x10+°°) 


;2.12xl0" 01 ) 




8.000 


5.35xl0+° 3 ) 


;6.62xl0+ 03 ) 


8.18xl0+ 03 ) 


(8.798x10+°°) 


;2.12xl0- 01 ) 




9.000 


8.68xl0+° 3 ) 


;i.o7xio+ 04 ) 


1.33xl0+ 04 ) 


(9.281x10+°°) 


;2.12xl0" 01 ) 




10.000 


1.30xl0+° 4 ) 


;i.60xio+ 04 ) 


1.98xl0+ 04 ) 


(9.681x10+°°) 


;2.12xl0" 01 ) 





C.4. 22 Ne(a, 7 ) 26 Mg 



Comments: the Ne+a reaction rates have previously been published in lLongland et al. 



(|2012l ). Here, we also include the upd ate input file and addition al plots. These r ates d iffer from 
those presented in llliadis et al.l ( 2010a ) in the following ways (see lLongland et al.1 (20L2j) for more 
details) : (i) the inflated weighted average method has been employed to account for the ambiguities 
in experimental data sets for res onance stren g ths a nd energies; (ii) normalization of the relative 
spectroscopic factors reported in iGiesen et al.1 ( 19931 ) is now performed with respect to the Ej. = 
1434 keV resonance rather than the Ef . afe = 831 keV resonance used in previous work; and (iii) a 
previous ( 6 Li,d) transfer measurement (IUgaldell2008l) suggested a spectr oscopic factor for the state 
corresponding to a resonance at E' a& = 588 keV; in llliadis et al.l ( 2010a ). this preliminary value is 
treated as an upper limit; in the new evaluation, this upper limit has been removed and replaced 
by C 2 S=1. For T > 1.25 GK, rates were extrapolated using Hauser-Feshbach calculations. 



22Ne(a,g)26Mg 



2 

10 



4.003 

21.991 

1.009 

0.0 



Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 
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o.o 

0.5 

10614.78 
11093.08 
1.25 

2 



Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 






1.0 
5000 





Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures; =NT for rate output at selected temperatures 



Nonresonant Contribution 
S (keVb) S ' (b) S " (b/keV) 
0.0 0.0 0.0 
0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 



fracErr 


Cutoff Energy (keV) 


0.0 


0.0 


0.0 


0.0 



Ecm 


DEcm 


wg 


Dwg 


J 


Gl 


DG1 


LI 


G2 


DG2 


L2 


G3 


DG3 


L3 


Exf 


Int 


78.37 


1.7 








4 


1.5e-46 


1.2e-46 


4 


3.0 


1.5 


1 











0.0 





703.78 


2.11 








2 


7.2e-6 


4.4e-7 


2 


3.0 


1.5 


1 


2.5e2 


1.7e2 


1 


0.0 


1 


826 . 04 


0.19 








4 


3.78e-6 


4.44e-7 


4 


3.0 


1.5 


1 


1.47e3 


8.0el 


2 


0.0 


1 


850.44 


0.21 








5 


4.36e-6 


9.09e-7 


5 


3.0 


1.5 


1 


6.55e3 


9.0el 


3 


0.0 


1 


893.31 


0.90 








1 


1.17e-4 


2.0e-5 


1 


3.0 


1.5 


1 


1.27e4 


2.5e3 


1 


0.0 


1 


911.16 


1.69 








1 


2.77e-4 


2.33e-5 


1 


3.0 


1.5 


1 


1 . 80e3 


9.0e2 


1 


0.0 


1 


1015.22 


1.69 








1 


2.83e-3 


3.33e-4 


1 


3.0 


1.5 


1 


1 . 35e4 


1.7e3 


1 


0.0 


1 


1133.66 


8.46 








1 


2.0e-2 


3.0e-3 


1 


3.0 


1.5 


1 


6.35e4 


8.5e3 


1 


0.0 


1 


1171.74 


3.38 








1 


1.67e-2 


2.33e-3 


1 


3.0 


1.5 


1 


2.45e4 


2.4e3 


1 


0.0 


1 


1213.0 


2.0 








2 


1.84e-l 


1.03e-l 


2 


3.0 


1.5 


1 


1.10e3 


2.5e2 





0.0 


1 


1280.0 


4.0 


2.0e-3 


2.0e-4 


1 


























1 


0.0 





1297.0 


3.0 








1 


1.89 


7.88e-l 


1 


3.0 


1.5 


1 


5.0e3 


2.0e3 


1 


0.0 


1 


1338.0 


3.0 








3 


6.48e-l 


3.33e-l 


3 


3.0 


1.5 


1 


4.0e3 


2.0e3 





0.0 


1 


1437 . 


3.0 








3 


8.58e-l 


5.81e-l 


3 


3.0 


1.5 


1 


3.0e3 


2.0e3 





0.0 


1 


1525.0 


3.0 








1 


1.67 


4.01e-l 


1 


3.0 


1.5 


1 


1.5e4 


2.0e3 


1 


0.0 


1 


1569.0 


7.0 











1.21el 


2.86 





3.0 


1.5 


1 


3.3e4 


5.0e3 


2 


0.0 


1 


1658.0 


7.0 











1 . 63e2 


3.49el 





3.0 


1.5 


1 


5.5e4 


1.0e4 


2 


0.0 


1 


1728.0 


4.0 











6.30e2 


1.22e2 





3.0 


1.5 


1 


3.5e4 


5.0e3 


2 


0.0 


1 



Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 

Note: . . .PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 



Ecm 


DEcm 


Jr 


Gl 


DG1 


LI 


PT 


G2 


DG2 




L2 


PT 


G3 


DG3 




L3 


PT 


Exf 


Int 


191.08 


0.15 


1 


1.25e-23 





1 


0.01 


3.0 


1.5 




1 



















0.0 





334.31 


0.1 


1 


1.20e-9 





1 


0.01 


3.0 


1.5 




1 



















0.0 





328.21 


2.0 


7 


3.70e-23 





7 


0.01 


3.0 


1.5 




1 



















0.0 





497 . 38 


0.08 


2 


9.28e-12 





2 


0.01 


1.73 


3.1e- 


-2 


1 





2.58e3 


2.40el 








0.0 


1 


548 . 16 


0.10 


2 


8.74e-8 





2 


0.01 


4.56 


2.9e- 


-1 


1 





4.64e3 


1.00e2 


1 





0.0 


1 


556 . 28 


0.16 


2 


1.25e-7 





2 


0.01 


3.0 


1.5 




1 





1.44 


1.6e- 


-1 


2 





0.0 





568 . 27 


0.19 


1 


2.08e-7 





1 


0.01 


3.0 


1.5 




1 





5.4e-l 


8.8e- 


-2 


1 





0.0 
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628.43 


0.10 


2 


9.46e-7 





2 


0.01 


7.42 


6.06-1 


1 





4.51e3 


1.07e2 


1 





0.0 


1 


659 . 32 


0.12 


2 


9.97e-7 





2 


0.01 


3.24 


3.5e-l 


1 





5.4e2 


5.4el 








0.0 


1 


665.11 


0.11 


4 


9.16e-8 





4 


0.01 


5.9e-l 


2.4e-l 


1 





1.51e3 


3.4el 


1 





0.0 


1 


670.81 


0.13 


1 


1.69e-6 





1 


0.01 


7.9e-l 


4.6e-l 


1 





1.26e3 


1.0e2 


1 





0.0 


1 


671.59 


0.12 


2 


1.02e-6 





2 


0.01 


4.26 


6.0e-l 


1 





1.28el 


6.0 


2 





0.0 


1 


674.36 


0.25 


2 


1.02e-6 





2 


0.01 


3.0 


1.5 


1 





1.54 


4.6e-l 


1 





0.0 





681.21 


0.13 


3 


7.37e-7 





3 


0.01 


3.31 


7.3e-l 


1 





8.06e3 


1.2e2 


1 





0.0 


1 


695.95 


0.35 


1 


1.76e-6 





1 


0.01 


3.0 


1.5 


1 





1.12 


4.0e-l 


1 





0.0 





711.34 


0.54 


1 


1.80e-6 





1 


0.01 


3.0 


1.5 


1 





6.0e-l 


3.2e-l 


1 





0.0 





713.40 


0.14 


1 


1.81e-6 





1 


0.01 


3.63 


4.7e-l 


1 





4.2e2 


8.6el 


1 





0.0 


1 


714.34 


0.55 


1 


1.81e-6 





1 


0.01 


3.0 


1.5 


1 





2.8 


1.0 


1 





0.0 





722.12 


0.56 


1 


1 . 83e-6 





1 


0.01 


3.0 


1.5 


1 





1.42 


5.6e-l 


1 





0.0 





729.15 


0.15 


2 


l.lle-6 





2 


0.01 


1.18 


2.7e-l 


1 





1.53e2 


4.2el 


1 





0.0 


1 


730.03 


0.16 


4 


6.16e-7 





4 


0.01 


1.82 


3.8e-l 


1 





4.13e3 


1.9e2 


3 





0.0 


1 


777.78 


0.16 


5 


1.51e-7 





5 


0.01 


3.0 


1.5 


1 





2.9e2 


1.9el 


2 





0.0 


1 



0.0 


0.0 


0.0 


0.0 


0.0 








0.0 


0.0 


0.0 


0.0 


0.0 


0.0 








0.0 



Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 

!+- 

0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1 . Reaction Rates using recent FEL results 

2. The doublet state at "330 keV has been included twice as upper limits for the spectroscopic factor 

3. The 703 keV resonance is treated as the same as seen in 22Ne(a,n) 
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Table 12:: Total thermonuclear reaction rates for 
22 Ne(a,n) 25 Mg. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


0.010 


0.00 xl0+ 00 


0.00x10+°° 


0.00x10+°° 


-5.708 xl0+° 2 


1.86x10+°° 


6.14xl0+ 01 


0.011 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-5.202 xl0+° 2 


1.86x10+°° 


6.10xl0+ 01 


0.012 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-4.781 xl0+° 2 


1.85x10+°° 


6.06xl0+ 01 


0.013 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-4.424xl0+° 2 


1.85x10+°° 


6.06xl0+ 01 


0.014 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-4.118xl0+° 2 


1.85x10+°° 


6.04xl0+ 01 


0.015 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-3.853 xl0+° 2 


1.85x10+°° 


6.02xl0+ 01 


0.016 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-3.621 xl0+° 2 


1.85x10+°° 


6.00xl0+ 01 


0.018 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-3.235 xl0+° 2 


1.85x10+°° 


6.00xl0+ 01 


0.020 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-2.925 xl0+° 2 


1.86x10+°° 


6.03xl0+ 01 


0.025 


0.00x10+°° 


0.00x10+°° 


0.00x10+°° 


-2.365 xl0+° 2 


1.87x10+°° 


6.28xl0+ 01 


0.030 


5.12xl0~ 88 


5.08xl0~ 87 


2.25xl0" 86 


-1.991 xl0+° 2 


1.90x10+°° 


6.65xl0+ 01 


0.040 


1.46xl0~ 67 


1.49xl0~ 66 


6.64xl0~ 66 


-1.519xl0+° 2 


1.94x10+°° 


7.20xl0+ 01 


0.050 


2.99xl0~ 55 


3.05 xl0~ 54 


1.36xl0~ 53 


-1.236xl0+° 2 


1.95x10+°° 


7.29xl0+ 01 


0.060 


4.92xl0~ 47 


4.87xl0~ 46 


2.17xl0~ 45 


-1.047xl0+° 2 


1.92x10+°° 


6.81 xlO+ 01 


0.070 


3.70xl0" 41 


3.48xl0~ 4 ° 


1.55xl0~ 39 


-9.117xl0+ 01 


1.84x10+°° 


5.72xl0+ 01 


0.080 


1.03xl0~ 36 


8.44xl0~ 36 


3.73 xlO" 35 


-8.101 xl0+ 01 


1.74x10+°° 


4.36xl0+ 01 


0.090 


3.23xl0~ 33 


2.19xl0~ 32 


9.43 xlO" 32 


-7.309 xl0+ 01 


1.62x10+°° 


3.11xl0+ 01 


0.100 


2.17xl0~ 3 ° 


1.20xl0~ 29 


4.92xl0~ 29 


-6.673 xl0+ 01 


1.50x10+°° 


2.12xl0+ 01 


0.110 


4.65xl0~ 28 


2.12xl0" 27 


8.22xl0~ 27 


-6.151 xl0+ 01 


1.39x10+°° 


1.43xl0+ 01 


0.120 


4.24xl0~ 26 


1.62xl0" 25 


5.82xl0~ 25 


-5.714xl0+ 01 


1.29x10+°° 


9.86x10+°° 


0.130 


1.94xl0 -24 


6.61 xl0~ 24 


2.14xl0~ 23 


-5.342 xl0+ 01 


1.19x10+°° 


6.93x10+°° 


0.140 


5.27xl0~ 23 


1.64xl0" 22 


4.81 xlO" 22 


-5.020xl0+ 01 


1.08x10+°° 


4.64x10+°° 


0.150 


9.94xl0~ 22 


2.74xl0~ 21 


7.18xl0" 21 


-4.737xl0+ 01 


9.62xl0~ 01 


2.96x10+°° 


0.160 


1.43xl0~ 2 ° 


3.39xl0~ 2 ° 


7.89 xlO" 20 


-4.484 xl0+ 01 


8.29xl0~ 01 


3.27x10+°° 


0.180 


1.61xl0~ 18 


2.74xl0~ 18 


5.01xl0~ 18 


-4.040 xl0+ 01 


5.53xl0~ 01 


1.26xl0+ 01 


0.200 


9.14xl0~ 17 


1.24xl0~ 16 


1.79xl0~ 16 


-3.660 xl0+ 01 


3.43xl0~ 01 


1.12xl0+ 01 


0.250 


1.68xl0" 13 


2.06xl0~ 13 


2.53 xl0~ 13 


-2.921 xl0+ 01 


2.06X10" 01 


3.49xl0~ 01 


0.300 


2.74X10" 11 


3.36xl0- n 


4.15xl0~ n 


-2.411 xl0+ 01 


2.06X10" 01 


4.10X10" 01 


0.350 


1.05xl0~° 9 


1.29xl0~° 9 


1.59xl0"° 9 


-2.046 xl0+ 01 


2.05xl0~ 01 


4.68xl0~ 01 


0.400 


1.64xl0~° 8 


2.00xl0~° 8 


2.45 xlO" 08 


-1.773xl0+ 01 


1.99xl0~ 01 


5.25xl0~ 01 


0.450 


1.42xl0~° 7 


1.71xl0~° 7 


2.07xl0~° 7 


-1.558xl0+ 01 


1.88xl0~ 01 


6.99xl0~ 01 


0.500 


8.51xl0~° 7 


1.00xl0~° 6 


1.19xl0~° 6 


-1.381 xl0+ 01 


1.68xl0~ 01 


1.13x10+°° 


0.600 


1.74xl0~° 5 


1.92xl0~° 5 


2.15xl0~° 5 


-1.085xl0+ 01 


1.07xl0~ 01 


2.15x10+°° 


0.700 


2.36xl0~° 4 


2.51xl0~° 4 


2.69 xlO" 04 


-8.287x10+°° 


6.70xl0~° 2 


1.13x10+°° 


0.800 


2.15xl0~° 3 


2.27xl0~° 3 


2.42xl0~° 3 


-6.084x10+°° 


5.79xl0~° 2 


3.45x10+°° 


0.900 


1.36xl0~° 2 


1.43xl0" 02 


1.51xl0~° 2 


-4.246x10+°° 


5.33xl0^° 2 


4.07x10+°° 
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Table 12 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


1.000 


6.34xl0~° 2 


6.64xl0-° 2 


6.98xl0"° 2 


-2.711 xl0+ 00 


4.82xl0~° 2 


3.59x10+°° 


1.250 


1.18xl0 +00 


1.22x10+°° 


1.27xl0 +00 


1.998xl0- 01 


3.88xl0~° 2 


1.63x10+°° 


1.500 


;i.09xl0 +01 ) 


1.14xl0+ 01 ) 


1.18xl0+ 01 ) 


(2.431 xl0+ 00 ) 


;3.89xl0~° 2 ) 




1.750 


6.79xl0+ 01 ) 


7.06xl0+ 01 ) 


7.34xl0+ 01 ) 


(4.257x10+°°) 


;3.89xl0~° 2 ) 




2.000 


2.92xl0+° 2 ) 


;3.04xl0 +02 ) 


3.16xl0+ 02 ) 


(5.717x10+°°) 


;3.89xl0~° 2 ) 




2.500 


;2.74xl0 +03 ) 


;2.85xl0 +03 ) 


2.96xl0+ 03 ) 


(7.953x10+°°) 


;3.89xl0~ 02 ) 




3.000 


'1.41xl0+ 04 ) 


1.46 xl0+ 04 ) 


1.52xl0+ 04 ) 


(9.590x10+°°) 


;3.89xl0~° 2 ) 




3.500 


4.96xl0+° 4 ) 


;5.16xl0 +04 ) 


5.37xl0+ 04 ) 


(1.085xl0+ 01 ) 


;3.89xl0~° 2 ) 




4.000 


1.36xl0+° 5 ) 


;i.4ixio+ 05 ) 


1.47xl0+ 05 ) 


(1.186xl0+ 01 ) 


;3.89xl0"° 2 ) 




5.000 


6.10xl0+ 05 ) 


;6.34xl0+ 05 ) 


6.59xl0+ 05 ) 


(1.336xl0+ 01 ) 


;3.89xl0" 02 ) 




6.000 


1.80xl0+° 6 ) 


1.88xl0+° 6 ) 


1.95xl0+ 06 ) 


(1.444xl0+ 01 ) 


;3.89xl0~° 2 ) 




7.000 


4.07xl0+ 06 ) 


;4.23xl0+ 06 ) 


4.40 xl0+ 06 ) 


(1.526xl0+ 01 ) 


'3.89xl0" 02 ) 




8.000 


7.70xl0+ 06 ) 


;8.oixio+ 06 ) 


8.32 xl0+ 06 ) 


(1.590xl0+ 01 ) 


;3.89xl0" 02 ) 




9.000 


1.28xl0 +07 ) 


;i.33xl0+ 07 ) 


1.39xl0+ 07 ) 


(1.640xl0+ 01 ) 


;3.89xl0" 02 ) 




10.000 


1.97xl0+ 07 ) 


;2.04xl0+ 07 ) 


2.12xl0+ 07 ) 


(1.683xl0+ 01 ) 


;3.89xl0" 02 ) 





C.5. 22 Ne(a,n) 25 Mg 
Comments: see Appendix IC. 41 for details on the new nuclear data for this reaction. 



22Ne(a,n)25Mg 



2 

10 



4.003 

21.991 

1.009 

0.0 

0.0 

0.5 

10614.78 

11093.08 

1.25 

3 



Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 



1.0 ! Minimum energy for numerical integration (keV) 

5000 ! Number of random samples (>5000 for better statistics) 

! =0 for rate output at all temperatures; =NT for rate output at selected temperatures 

Nonresonant Contribution 
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S(keVb) 


S'Cb) 


S"(b/keV) 


fracErr 


Cutoff Energy (keV) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 



Ecm 


DEcm 


wg 




Dwg 




J 


Gl 




DG1 




LI 


G2 


DG2 


L2 


G3 


DG3 


L3 


Exf 


Int 


703.78 


2.11 












2 


2. 


36e- 


-5 


2. 


2e- 


-6 


2 


2. 


5e2 


1. 


7e2 


2 


3 


1. 


5 


1 


0.0 


1 


826 . 04 


0.19 












4 


3. 


78e- 


-6 


4. 


4e- 


-7 


4 


1. 


47e3 


8. 


Oel 


4 


3 


1. 


5 


1 


0.0 


1 


850.44 


0.21 












5 


4. 


36e- 


-6 


9. 


le- 


-7 


5 


6. 


55e3 


9. 


Oel 


5 


3 


1. 


5 


1 


0.0 


1 


893.31 


0.90 












1 


1. 


17e- 


-4 


2. 


Oe- 


-5 


1 


1. 


27e4 


2. 


5e3 


1 


3 


1, 


5 


1 


0.0 


1 


911.16 


1.69 












1 


2. 


77e- 


-4 


2. 


3e- 


-5 


1 


1. 


80e3 


9. 


0e2 


1 


3 


1. 


5 


1 


0.0 


1 


1015.22 


1.69 












1 


2. 


83e- 


-3 


3. 


3e- 


-4 


1 


1. 


35e4 


1. 


7e3 


1 


3 


1, 


5 


1 


0.0 


1 


1133.66 


8.46 












1 


2. 


Oe-2 


3. 


Oe- 


-3 


1 


6. 


35e4 


8. 


5e3 


1 


3 


1. 


5 


1 


0.0 


1 


1171.74 


3.38 












1 


1. 


67e- 


-2 


2. 


3e- 


-3 


1 


2. 


45e4 


2. 


4e3 


1 


3 


1. 


5 


1 


0.0 


1 


1213.19 


2.34 












2 


2. 


13e- 


-1 


8. 


4e- 


-3 


2 


1. 


10e3 


2. 


5e2 


2 


3 


1. 


5 


1 


0.0 


1 


1247.88 


2.54 












1 


1. 


5e-2 


1. 


Oe- 


-2 


1 


2. 


45e4 


3. 


4e3 


1 


3 


1. 


5 


1 


0.0 


1 


1264.80 


2.54 


3.9e 


-1 


5.7e- 


-2 


1 
















1 












1 










1 


0.0 





1275.80 


2.54 


5.6e 


-1 


6.0e- 


-2 


1 
















1 












1 










1 


0.0 





1295.25 


2.54 


1.5 




1.6e- 


-1 


1 
















1 












1 










1 


0.0 





1336.71 


2.54 


2.9 




3.0e- 


-1 


1 
















1 












1 










1 


0.0 





1437 . 38 


2.54 


6.0 




7.7e- 


-1 


1 
















1 












1 










1 


0.0 





1499.99 


4.23 


1.0 




2.4e- 


-1 


1 
















1 












1 










1 


0.0 





1526.22 


2.54 


3.0 




3.4e- 


-1 


1 
















1 












1 










1 


0.0 





1569.36 


6.77 


9.0e 


-1 


2.1e- 


-1 


1 
















1 












1 










1 


0.0 





1649.74 


8.46 


3.1e+l 


8.5 




1 
















1 












1 










1 


0.0 





1730.95 


6.77 


2.0e 


f2 


3.3e+l 


1 
















1 












1 










1 


0.0 





1820.63 


8.46 


2.8e+l 


7.0 




1 
















1 












1 










1 


0.0 





1936 . 54 


12.69 


1.2e- 


<-2 


4.5e+l 


1 
















1 












1 










1 


0.0 






Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 



Note: . 


. .PT=<B> 


for 


g-rays [e 


inter 


II 


upper _limit 0.0"] 


; for each resonance: # uppe 


r limits 


< # 


open 


channels ! 


Ecm 


DEcm 


Jr 


Gl 


DG1 


LI 


PT 


G2 


DG2 


L2 


PT 


G3 


DG3 


L3 


PT 


Exf 


Int 


497 . 38 


0.08 


2 


9.28e-12 





2 


0.01 


2 . 58e3 


2.4el 








1.73 


3e-2 


1 





0.0 


1 


548 . 16 


0.10 


2 


3.80e-8 





2 


0.01 


4.64e3 


1.0e2 


1 





4.56 


0.29 


1 





0.0 


1 


556 . 28 


0.16 


2 


1.50e-8 





2 


0.01 


1.44 


1.6e-l 


2 





3.0 


1.5 


1 





0.0 





568 . 27 


0.19 


1 


2.08e-7 





1 


0.01 


5.40e-l 


8.8e-2 


1 





3.0 


1.5 


1 





0.0 





628.43 


0.10 


2 


2.40e-8 





2 


0.01 


4.51e3 


l.le2 


1 





7.42 


0.60 


1 





0.0 


1 


659.32 


0.12 


2 


2.20e-8 





2 


0.01 


5.40e2 


5.4el 








3.24 


0.35 


1 





0.0 


1 


665.11 


0.11 


4 


1.44e-8 





4 


0.01 


1.51e3 


3.4el 


1 





5.9e-l 


2.4e-l 


1 





0.0 


1 


670.81 


0.13 


1 


2.57e-8 





1 


0.01 


1 . 26e3 


1.0e2 


1 





7.9e-l 


4.6e-l 


1 





0.0 


1 


671.59 


0.12 


2 


1 . 54e-8 





2 


0.01 


1.28el 


6.0 


2 





4.26 


0.60 


1 





0.0 


1 


674 . 36 


0.25 


2 


1 . 54e-8 





2 


0.01 


1.54 


0.46 


1 





3.0 


1.5 


1 





0.0 





681.21 


0.13 


3 


1.43e-8 





3 


0.01 


8 . 06e3 


1.2e2 


1 





3.31 


0.73 


1 





0.0 


1 


695.95 


0.35 


1 


5.34e-9 





1 


0.01 


1.12 


0.40 


1 





3.0 


1.5 


1 





0.0 





711.34 


0.54 


1 


4.11e-8 





1 


0.01 


6.0e-l 


3.2e-l 


1 





3.0 


1.5 


1 





0.0 





713.40 


0.14 


1 


1.67e-7 





1 


0.01 


4.24e2 


8.6el 


1 





3.63 


0.47 


1 





0.0 


1 


714.34 


0.55 


1 


4.12e-8 





1 


0.01 


2.8 


1.0 


1 





3.0 


1.5 


1 





0.0 
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722.12 


0.56 


1 


4.17e-8 





1 


0.01 


1.42 


0.56 


1 





3.0 


1.5 


1 





0.0 





729.15 


0.15 


2 


4.00e-8 





2 


0.01 


1.53e2 


4.2el 


1 





1.18 


0.27 


1 





0.0 


1 


730.03 


0.16 


4 


4.68e-9 





4 


0.01 


4.13e3 


1.9e2 


3 





1.82 


0.38 


1 





0.0 


1 


777 . 78 


0.16 


5 


3.34e-9 





5 


0.01 


2.90e2 


1.9el 


2 





3.0 


1.5 


1 





0.0 


1 



Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 

!+- 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0.0 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. a,n measurements including results from recent FEL run 

2. The 703 keV resonance is treated as the same as seen in 22Ne(a,g) 
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2 Ne(a, n) 25 Mg 



- 




T9 = 0.03 
A-D = 121 



- 




T9 = 0.06 
A-D = 120.7 



CO 

O 




- 










T9 = 0.3 
A-D = 0.2664 



2.0B-11 2.5G-11 3.08-11 3.5e-11 



- 






T9 = 0.6 
A-D = 1 .65 



l iJJIPI 






T9 = 1 
A-D = 0.5268 



= 0.075 0.080 0.085 0.090 0.095 0.100 



Reaction rate (cm 3 mol 1 s 1 ) 
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to 

<D 

DC 



to 

<D 



22 Ne(a, n) 25 Mg 




0.01 



0.10 



1.00 



10.00 



to 
a: 

<D 
"to 



to 
<D 

a: 



0.01 




0.10 1.00 

Temperature (GK) 



10.00 
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Table 13:: Total thermonuclear reaction rates for 
22 Na(>, 7) 23 Mg. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


0.010 


1.22xl0~ 3 ° 


2.52X10" 30 


5.16X10" 30 


-6.816xl0+ 01 


7.17xl0^ 01 


6.24xl0" 01 


0.011 


1.05xl0~ 28 


1.96xl0~ 28 


3.66 xlO" 28 


-6.380 xl0+ 01 


6.29xl0~ 01 


6.85xl0~ 01 


0.012 


4.20xl0~ 27 


7.33xl0~ 27 


1.27xl0" 26 


-6.018xl0+ 01 


5.59xl0~ 01 


6.84xl0~ 01 


0.013 


9.35xl0~ 26 


1.57xl0 -25 


2.56xl0" 25 


-5.712xl0+ 01 


5.07X10" 01 


5.62X10" 01 


0.014 


1.33xl0~ 24 


2.13xl0~ 24 


3.36 xlO" 24 


-5.451 xl0+ 01 


4.68X10" 01 


4.65X10" 01 


0.015 


1.30xl0~ 23 


2.03xl0~ 23 


3.11xl0" 23 


-5.226 xl0+ 01 


4.41 xlO~ 01 


3.96xl0~ 01 


0.016 


9.55xl0~ 23 


1.44xl0~ 22 


2.19xl0" 22 


-5.029xl0+ 01 


4.22xl0~ 01 


2.66xl0~ 01 


0.018 


2.55xl0~ 21 


3.76xl0~ 21 


5.64xl0" 21 


-4.702 xl0+ 01 


4.03X10" 01 


3.99xl0- 01 


0.020 


3.47xl0~ 2 ° 


5.10xl0~ 2 ° 


7.61xl0~ 2 ° 


-4.442 xl0+ 01 


3.98X10" 01 


3.48X10" 01 


0.025 


4.19xl0~ 18 


6.06xl0~ 18 


8.92xl0~ 18 


-3.965 xl0+ 01 


3.80xl0~ 01 


4.33xl0~ 01 


0.030 


1.54xl0~ 16 


2.15xl0~ 16 


2.97xl0~ 16 


-3.608 xl0+ 01 


3.31xl0~ 01 


3.12xl0~ 01 


0.040 


3.03 xl0~ 14 


4.34xl0~ 14 


6.24xl0~ 14 


-3.077xl0+ 01 


3.66X10" 01 


2.86X10" 01 


0.050 


8.55xl0~ 13 


1.33xl0" 12 


2.06xl0~ 12 


-2.735xl0+ 01 


4.45xl0~ 01 


8.19xl0~ 01 


0.060 


8.33xl0~ 12 


1.37X10 -11 


2.20X10" 11 


-2.502 xl0+ 01 


4.84X10" 01 


5.18xl0~ 01 


0.070 


1.47xl0~ 10 


1.95xl0~ 10 


2.60xl0~ 10 


-2.236 xl0+ 01 


2.79xl0- 01 


1.18x10+°° 


0.080 


5.29xl0~° 9 


6.79xl0~° 9 


8.88xl0~° 9 


-1.880xl0+ 01 


2.50xl0~ 01 


3.54x10+°° 


0.090 


1.14X10" 07 


1.43xl0-° 7 


1.82xl0~° 7 


-1.575xl0+ 01 


2.30X10" 01 


1.62x10+°° 


0.100 


1.35xl0~° 6 


1.66xl0~° 6 


2.06xl0~° 6 


-1.330xl0+ 01 


2.12xl0~ 01 


7.75xl0~ 01 


0.110 


l.OlxlO" 05 


1.23xl0"° 5 


1.50xl0~° 5 


-1.131xl0+ 01 


2.00xl0~ 01 


4.32xl0~ 01 


0.120 


5.34xl0~° 5 


6.45xl0~° 5 


7.76xl0~° 5 


-9.650 xl0+°° 


1.90xl0~ 01 


3.05xl0~ 01 


0.130 


2.17xl0~° 4 


2.61 xlO~° 4 


3.11xl0"° 4 


-8.254x10+°° 


1.83X10" 01 


2.74xl0~ 01 


0.140 


7.21xl0~° 4 


8.59xl0~° 4 


1.02xl0~° 3 


-7.062x10+°° 


1.76xl0~ 01 


2.79xl0~ 01 


0.150 


2.03xl0~° 3 


2.40xl0~° 3 


2.83 xlO" 03 


-6.032x10+°° 


1.71xl0~ 01 


3.29xl0~ 01 


0.160 


5.01xl0~° 3 


5.90xl0~° 3 


6.94xl0"° 3 


-5.133x10+°° 


1.66xl0~ 01 


3.93xl0~ 01 


0.180 


2.25xl0~° 2 


2.63xl0~° 2 


3.07xl0~° 2 


-3.639x10+°° 


1.57xl0~ 01 


5.19xl0~ 01 


0.200 


7.49 xl0~° 2 


8.67xl0~° 2 


l.OOxlO^ 01 


-2.445x10+°° 


1.50xl0~ 01 


6.04xl0~ 01 


0.250 


6.56xl0- 01 


7.46X10" 01 


8.53X10" 01 


-2.911xl0- 01 


1.35X10" 01 


6.39xl0~ 01 


0.300 


2.81x10+°° 


3.16x10+°° 


3.59x10+°° 


1.153x10+°° 


1.26X10" 01 


6.20X10" 01 


0.350 


7.93xl0 +00 


8.94x10+°° 


1.01 xlO+ 01 


2.191x10+°° 


1.24xl0~ 01 


4.25xl0~ 01 


0.400 


1.73xl0+ 01 


1.96xl0+ 01 


2.21xl0+ 01 


2.974x10+°° 


1.24xl0" 01 


4.83xl0~ 01 


0.450 


3.19xl0+ 01 


3.62xl0+ 01 


4.08 xlO+ 01 


3.587x10+°° 


1.24xl0~ 01 


6.08xl0~ 01 


0.500 


5.25xl0 +01 


5.94xl0+ 01 


6.69 xlO+ 01 


4.084x10+°° 


1.24xl0~ 01 


6.73X10" 01 


0.600 


1.14xl0+ 02 


1.28xl0+° 2 


1.44xl0+ 02 


4.852x10+°° 


1.20xl0~ 01 


8.94xl0~ 01 


0.700 


2.05xl0+° 2 


2.28xl0+° 2 


2.55xl0+ 02 


5.434x10+°° 


l.llxlO" 01 


1.24x10+°° 


0.800 


3.30xl0 +02 


3.65xl0+° 2 


4.04 xl0+ 02 


5.902x10+°° 


1.02xl0~ 01 


1.50x10+°° 


0.900 


4.92 xl0+ 02 


5.39xl0+° 2 


5.92 xl0+ 02 


6.292x10+°° 


9.32xl0" 02 


1.57x10+°° 
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Table 13 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a A-D 


1.000 


6.91 xl0 +02 


7.52xl0+ 02 


8.19xl0+ 02 


6.624x10+°° 


8.59xl0" 02 1.51x10+°° 


1.250 


;i.40xl0 +03 ) 


1.52xl0+° 3 ) 


1.65xl0 +03 ) 


(7.327x10+°°) 


;8.01xl0~° 2 ) 


1.500 


'2.41 xlO+ 03 ) 


2.61 xlO+ 03 ) 


2.82xl0+° 3 ) 


(7.865x10+°°) 


'8.01xl0~° 2 ) 


1.750 


3.57xl0+ 03 ) 


3.87xl0+° 3 ) 


4.20xl0+° 3 ) 


(8.262x10+°°) 


;8.01xl0~° 2 ) 


2.000 


4.85xl0+° 3 ) 


5.25xl0+° 3 ) 


5.69xl0+° 3 ) 


(8.566x10+°°) 


'8.01xl0~° 2 ) 


2.500 


7.53xl0+° 3 ) 


8.15xl0+° 3 ) 


8.83xl0 +03 ) 


(9.006x10+°°) 


;8.01xl0~° 2 ) 


3.000 


1.01xl0+ 04 ) 


1.09xl0 +04 ) 


1.18xl0+ 04 ) 


(9.297x10+°°) 


;8.oixio~° 2 ) 


3.500 


1.24xl0+ 04 ) 


1.35xl0+ 04 ) 


1.46xl0+° 4 ) 


(9.507x10+°°) 


;8.oixio~° 2 ) 


4.000 


T.45xl0+° 4 ) 


1.58xl0+° 4 ) 


1.71 xl0+ 04 ) 


(9.665x10+°°) 


'8.01xl0"° 2 ) 


5.000 


1.81xl0+ 04 ) 


1.97xl0+ 04 ) 


2.13xl0+ 04 ) 


(9.886x10+°°) 


;8.oixio-° 2 ) 


6.000 


2.10xl0+° 4 ) 


2.28xl0 +04 ) 


2.47xl0+° 4 ) 


(1.003xl0+ 01 ) 


'8.01xl0~° 2 ) 


7.000 


2.33xl0+° 4 ) 


2.52xl0+° 4 ) 


2.73 xl0+ 04 ) 


(1.013xl0+ 01 ) 


'8.01xl0"° 2 ) 


8.000 


2.51xl0+° 4 ) 


2.72 xl0+ 04 ) 


2.95xl0+° 4 ) 


(1.021xl0+ 01 ) 


'8.01xl0-° 2 ) 


9.000 


2.66xl0+° 4 ) 


2.88xl0+ 04 ) 


3.12xl0+ 04 ) 


(1.027xl0+ 01 ) 


;8.oixio-° 2 ) 


10.000 


2.82xl0+ 04 ) 


3.05 xl0+ 04 ) 


3.30xl0+ 04 ) 


(1.033xl0+ 01 ) 


;8.oixio-° 2 ) 



C.6. 22 Na(p, 7 ) 23 Mg 



values for resonance stren gths (ISallaska et al 



Comments: The input file from Iliadis et al. (2010b1 has been upd ated with new 



2011. 



201G ) and energies (ISallaska et al 



Saastamoinen et al.ll201ll ). A new Q value (7580.72 ± 0.71 keV) has also been adopted ( Wang et al 



20121 . This Q- value was used to calculate proton energi e s from the excitation energies of 



y measured 



2011 



2010 



Jenkins et al 



( 2004) 1 for resonances at 43, 66, and 1 89 ke V, of lEndtJ (J1998J1 for re sonances at 480, 494, 612, 706 



and 760 keV, and of iPerajarvi et al.1 (J2000J) and llaco 



) et al 



(200 



al for 



the resonance at 221 keV 



(weighted average). The energy analysis outlined in (jSallaska et al.ll201ll : ISallaskall2010i 1 was per- 
formed with the new Q value to extract proton energies for resonances at 204, 274, 434, and 583 keV, 
and the results changed by 0.1 keV or less. These values were combined in a weight e d average with 



proton energies calculated via exci tation energ i es fro m ISaastamoinen et al.l (|201ll ) ; iJenkins et al 



(|2004l 1 . The value calculated from llacob et al.l (|2006l 1 was also included for the resonance at 204 
keV. 

All resonance strengths an d proton partia l widths, inc luding upper limits, have been taken 
from the direct measurement of ISallaska et al.l (|201ll . |201Q| 1 , with the exception of resonances at 
43 and 66 keV. For two resonances measured by Sallaska et al. (204 and 583 keV), the reported 
strengths have an inflated upper uncertainty to account for possible unobserved contributions to 
the rate. To calculate the reaction rate with the present method, we simply take the symmetric, 
measured result. Because of the slight energy change for these resonances, the penetration factors 
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were recalculated, resulting in slightly different particle partial widths from llliadis et al.l (I2010bl) 



The r eso nance strengths not dir ectly measured by Sallaska et al. but measured by 



2011 



199C) or lStegmiiller et al.l (|1996l ) have been scaled by a factor of 2.5, as discussed in lSallaska et al 



Seutheetal 



This i nclud es resonances at 480, 706, and 760 keV. Although not quoted explicitly in 



Sallaska et al.l ()201ll ). this me thodology has also b e applied to potential resonances at 493 and 612 
keV, originally investigated in lSeuthe et al.l ()1990l ). For T > 1.0 GK, rates were extrapolated using 
Hauser-Feshbach calculations. 



The direct capture component has roughly been estimated in ISeuthe et al.l (|1990l ). As a test, 
S(0) was increased and decreased by a factor of 5, which changed the total rate by no more than 
3%. The conclusion is the direct capture process is insignificant. 



22Na(p,g)23Mg 



1 ! Zproj 
11 ! Ztarget 

! Zexitparticle (=0 when only 2 channels open) 

1.0078 ! Aproj 

21.9944 ! Atarget 

! Aexitparticle (=0 when only 2 channels open) 

0.5 ! Jpro j 

3.0 ! Jtarget 

0.0 ! Jexitparticle (=0 when only 2 channels open) 

7580.72 ! projectile separation energy 

0.0 ! exit particle separation energy (=0 when only 2 channels open) 

1.25 ! Radius parameter R0 (fm) 

2 ! Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 



1.0 
5000 





Minimum energy for numerical integration (keV) 

Number of random samples (>5000 for better statistics) 

=0 for rate output at all temperatures ; =NT for rate output at selected temperatures 



Non-resonant contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

1.8el 0.0 0.0 0.4 1500.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 

Ecm DEcm wg Dwg Jr Gl DG1 LI G2 DG2 L2 G3 DG3 L3 Exf Int 

7.8e-17 3.2e-17 2 

1.6e-12 6.6e-13 2 













42 . 7 1.1 

66.2 2.7 

204.4 0.9 

273 . 8 0.9 

434.4 0.7 
480 . 3 2.1 

582 . 5 0.7 



4.5 

1.5 

5.7e-3 0.9e-3 

3.9e-2 0.8e-2 

1.7e-l 0.2e-l 

9.3e-2 3.6e-2 

5.9e-l 0.7e-l 



1.6e-l 





8e- 


-1 


1 











0.0 


1 


2.0e-2 


1 


Oe- 


-2 


1 











0.0 



























0.0 



























0.0 



























0.0 



























0.0 



























0.0 






221.3 


1.6 


2.5 


0.0016 





493.3 


6.0 


2.5 


. 044 





612.3 


8.0 


2.5 


0.052 


















2715.0 

















0.0 

















0.0 

















0.0 
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706.3 2.1 9.1e-l 1.76-1 00000. 00 

760.3 2.1 2.4e-l 0.8e-l 0000 0. 

Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 
Note: ...PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 
Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 

188.5 1.2 2.5 0.0012 0.0045 0.33 0.16 1 

0.0045 0.2 0.1 1 

0.0045 0.2 0.1 1 

0.0045 0.2 0.1 1 

Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf 

!+- 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Gg=2.0e-2 eV for Er=66.6 keV resonance is a guess (not important since Gp«Gg) . 

2. Gg=2.0e-1 eV for Er=222, 494 and 614 keV resonances is a guess; we assume for these undetected resonances 
that Gp«Gg (assumption most likely inconsequential for total rates) . 

3. Spin and parity of Er=494 and 614 keV resonances unknown; for upper limit contributions we assume s-waves 
(Jp=5/2+) . 

4. Strengths taken directly from Sallaska et al . (2010, 2011), except at Er=43 and 66 keV; strengths at 480, 
493, 612, 706, and 760 have been scaled up from Seuthe et al . (1990) by 2.5, per Sallaska et al. (2011). 

5. Energies are from Sallaska et al. (2010, 2011), Perajarvi et al . (2000), Iacob et al. (2006), 
Saastamoinen et al . (2011), Jenkins (2004), and Endt (1998), with a Q-value from Audi & Meng (2011). 
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2 Na(p, y) 23 Mg 



to 


T9 = 0.03 


o - 
d 




A-D = 0.3122 


o 

o - 


l 




o 


II 


CM 

o 

o - 


if 


o 


I \ 


o 

o 
o - 


m mi in ii I I I 



4e-16 6e-16 8e-16 



to 



CB 
_Q 
O 



L*T 


T9 = 0.1 
A-D = 0.7755 

iHiky! i 



2.0e-06 3.0e-06 




iuJIP! 




T9 = 0.6 
A-D = 0.8939 







T9=1 


o - 




A-D = 1 .507 


o 






o 
o - 






o 






o 
a - 






o 




o 
o - 




o 




Q 

o - 




o 




o 

o 





Reaction rate (cm 3 mol 1 s 1 ) 
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Na(p,7) 23 Mg 



o 
a: 






2.5- 


\ 
\ 


/ v 


2.0- 


\ 
\\ 
\\ 


1 x 


1.5- 


^^^-^v 


^S \ S. 


1.0- 








0.5- 


~-<_-- " 






.'J_i :j d : 



10 



10 



10" 



10' 




Temperature (GK) 
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Table 14:: Total thermonuclear reaction rates for 
29 P(p, 7 ) 30 S. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


0.010 


4.57xl0" 42 


6.75xl0~ 42 


9.90xl0~ 42 


-9.480 xl0+ 01 


3.87xl0" 01 


2.82xl0" 01 


0.011 


1.79xl0~ 4 ° 


2.60xl0~ 4 ° 


3.85 xl0" 4 ° 


-9.114xl0+ 01 


3.80xl0~ 01 


4.17xl0~ 01 


0.012 


4.50xl0~ 39 


6.70xl0~ 39 


9.77xl0" 39 


-8.790 xl0+ 01 


3.93X10" 01 


3.56xl0~ 01 


0.013 


8.19xl0~ 38 


1.21xl0~ 37 


1.74xl0" 37 


-8.501 xl0+ 01 


3.82X10" 01 


7.03xl0- 01 


0.014 


1.09xl0~ 36 


1.64xl0~ 36 


2.39xl0~ 36 


-8.241 xl0+ 01 


3.89X10" 01 


4.17xl0- 01 


0.015 


1.19xl0~ 35 


1.73xl0~ 35 


2.54xl0" 35 


-8.004xl0+ 01 


3.84xl0~ 01 


5.09xl0~ 01 


0.016 


1.03xl0~ 34 


1.52xl0~ 34 


2.22xl0" 34 


-7.787xl0+ 01 


3.86xl0~ 01 


2.10xl0~ 01 


0.018 


4.79xl0~ 33 


7.01xl0 -33 


1.03xl0" 32 


-7.404xl0+ 01 


3.83xl0- 01 


3.63xl0- 01 


0.020 


1.30xl0~ 31 


1.90xl0" 31 


2.79xl0~ 31 


-7.074xl0+ 01 


3.86xl0~ 01 


1.88xl0- 01 


0.025 


9.70xl0~ 29 


1.44xl0~ 28 


2.10xl0" 28 


-6.412xl0+ 01 


3.87xl0~ 01 


4.47xl0~ 01 


0.030 


1.51xl0~ 26 


2.22xl0~ 26 


3.25xl0" 26 


-5.907xl0+ 01 


3.82xl0~ 01 


1.81xl0~ 01 


0.040 


2.34xl0~ 23 


3.39xl0~ 23 


4.98 xlO" 23 


-5.173xl0+ 01 


3.83xl0- 01 


2.73xl0~ 01 


0.050 


4.28xl0~ 21 


6.25xl0" 21 


9.16xl0~ 21 


-4.652 xl0+ 01 


3.87xl0~ 01 


7.71 xlO~ 01 


0.060 


2.22xl0~ 19 


3.29xl0~ 19 


4.82xl0~ 19 


-4.256 xl0+ 01 


3.84xl0~ 01 


4.39xl0~ 01 


0.070 


5.70xl0~ 18 


8.23xl0~ 18 


1.18xl0~ 17 


-3.934xl0+ 01 


3.72xl0~ 01 


2.75xl0~ 01 


0.080 


1.50xl0~ 16 


2.04xl0~ 16 


2.78xl0~ 16 


-3.613xl0+ 01 


3.13xl0~ 01 


6.52xl0~ 01 


0.090 


5.96xl0~ 15 


8.77xl0~ 15 


1.33xl0" 14 


-3.235xl0+ 01 


4.08xl0~ 01 


1.93x10+°° 


0.100 


1.87xl0~ 13 


2.82xl0~ 13 


4.35 xlO" 13 


-2.889 xl0+ 01 


4.30xl0~ 01 


6.69xl0~ 01 


0.110 


3.39xl0~ 12 


5.09xl0~ 12 


7.79xl0" 12 


-2.600 xl0+ 01 


4.23xl0~ 01 


5.10xl0" 01 


0.120 


3.80X10" 11 


5.65xl0~ n 


8.59xl0~ n 


-2.359xl0+ 01 


4.15xl0~ 01 


4.87xl0~ 01 


0.130 


2.91xl0~ 10 


4.32xl0~ 10 


6.50xl0~ 10 


-2.156xl0+ 01 


4.07xl0~ 01 


5.06xl0~ 01 


0.140 


1.66xl0~° 9 


2.45xl0~° 9 


3.67xl0~° 9 


-1.982xl0+ 01 


4.01 xlO~ 01 


5.28xl0~ 01 


0.150 


7.49xl0~° 9 


1.10xl0~° 8 


1.64xl0"° 8 


-1.832xl0+ 01 


3.95xl0~ 01 


5.77xl0~ 01 


0.160 


2.78xl0~° 8 


4.08xl0~° 8 


6.04xl0"° 8 


-1.701xl0+ 01 


3.90xl0~ 01 


6.68xl0~ 01 


0.180 


2.49xl0~° 7 


3.58xl0~° 7 


5.29xl0~° 7 


-1.483xl0+ 01 


3.77xl0~ 01 


1.08x10+°° 


0.200 


1.46X10" 06 


2.07xl0~° 6 


3.00 xlO" 06 


-1.308xl0+ 01 


3.59xl0~ 01 


1.80x10+°° 


0.250 


3.98xl0~° 5 


5.31 xl0~° 5 


7.22xl0~° 5 


-9.833 xl0+°° 


3.01xl0~ 01 


3.74x10+°° 


0.300 


4.16xl0~° 4 


5.33xl0-° 4 


6.88 xlO" 04 


-7.532 xl0+°° 


2.57xl0~ 01 


1.82x10+°° 


0.350 


2.39xl0~° 3 


3.03xl0~° 3 


3.86 xlO" 03 


-5.798 xl0+°° 


2.44xl0~ 01 


4.95xl0~ 01 


0.400 


9.07xl0~° 3 


1.16xl0~° 2 


1.48xl0~° 2 


-4.458 xl0+°° 


2.48xl0~ 01 


4.46xl0~ 01 


0.450 


2.60xl0~° 2 


3.31 xl0~° 2 


4.30 xlO" 02 


-3.401 xl0+°° 


2.56xl0~ 01 


6.09xl0~ 01 


0.500 


6.03xl0~° 2 


7.73xl0~° 2 


1.02xl0~ 01 


-2.553x10+°° 


2.63xl0~ 01 


6.65xl0~ 01 


0.600 


2.11X10" 01 


2.75X10" 01 


3.62X10" 01 


-1.287x10+°° 


2.71xl0~ 01 


5.70xl0~ 01 


0.700 


5.17xl0~ 01 


6.73X10" 01 


8.88xl0~ 01 


-3.918xl0~ 01 


2.73xl0~ 01 


4.91 xlO~ 01 


0.800 


1.01xl0 +0 ° 


1.31x10+°° 


1.73x10+°° 


2.775xl0~ 01 


2.69xl0~ 01 


4.67xl0~ 01 


0.900 


1.73 xl0+ 00 


2.22x10+°° 


2.90x10+°° 


8.021 xlO" 01 


2.61xl0~ 01 


4.17xl0^ 01 
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Table 14 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a 


A-D 


1.000 


2.66xl0 +00 


3.41x10+°° 


4.41x10+°° 


1.230x10+°° 


2.52xl0~ 01 


3.43xl0~ 01 


1.250 


6.09xl0 +00 


7.65x10+°° 


9.62x10+°° 


2.037x10+°° 


2.32xl0" 01 


2.60xl0~ 01 


1.500 


1.10xl0 +01 


1.36xl0+ 01 


1.70xl0+ 01 


2.616x10+°° 


2.20xl0~ 01 


3.78xlQ- 01 


1.750 


1.70xl0+ 01 ) 


(2.11xl0+ 01 ) 


2.62xl0+ 01 ) 


(3.048x10+°°) 


;2.17xl0~ 01 ) 




2.000 


'2.48xl0+ 01 ) 


(3.08xl0+ 01 ) 


3.82xl0+ 01 ) 


(3.426x10+°°) 


;2.17xl0~ 01 ) 




2.500 


;4.29xl0+ 01 ) 


(5.33xl0+ 01 ) 


6.62xl0+ 01 ) 


(3.976x10+°°) 


;2.17xl0~ 01 ) 




3.000 


6.33xl0+ 01 ) 


(7.86xl0+ 01 ) 


9.75xl0+ 01 ) 


(4.364x10+°°) 


;2.17xl0~ 01 ) 




3.500 


;8.38xl0+ 01 ) 


(1.04xl0+° 2 ) 


1.29xl0+° 2 ) 


(4.645x10+°°) 


;2.17xl0~ 01 ) 




4.000 


'1.05xl0+ 02 ) 


(1.31xl0+° 2 ) 


1.63xl0+° 2 ) 


(4.875x10+°°) 


;2.17xl0- 01 ) 




5.000 


;i.47xl0+ 02 ) 


(1.83xl0+° 2 ) 


2.27xl0+° 2 ) 


(5.208x10+°°) 


;2.17xl0" 01 ) 




6.000 


1.88xl0+ 02 ) 


(2.34xl0+° 2 ) 


2.90xl0+° 2 ) 


(5.454x10+°°) 


;2.17xl0~ 01 ) 




7.000 


2.27xl0+ 02 ) 


(2.81xl0+° 2 ) 


3.50xl0+° 2 ) 


(5.640x10+°°) 


;2.17xl0" 01 ) 




8.000 


;2.63xl0+ 02 ) 


(3.27xl0+° 2 ) 


4.06xl0+° 2 ) 


(5.789x10+°°) 


;2.17xl0~ 01 ) 




9.000 


2.97xl0+ 02 ) 


(3.69xl0+° 2 ) 


4.58xl0+° 2 ) 


(5.911x10+°°) 


;2.17xl0~ 01 ) 




10.000 


3.36xl0+ 02 ) 


(4.17xl0+° 2 ) 


5.18xl0+° 2 ) 


(6.034x10+°°) 


;2.17xl0" 01 ) 





C.7. 29 P(p, 7 ) JU S 



30 s 



Comments: The input file from llliadis et al.l (l2010bh has been updated w i th ne w information on 
the en ergies of astrophysically important resonances ISetoodehnia et al.l (J2010l . l201ll ). ISetoodehnia et al 
(12010]') reaffirms the exist ence of a level predicted to be 3+ (E x = 4693(5) keV), previously observed 



in 



Bardayan et al.1 (120071 ). and measured the energy of the potential 2+ level (E x = 481 4(3) keV) for 



the firs t time, decreasing the uncertainty in the energy estimated u sing the IMME from llliadis et al 



(hoiObh bv a factor of 10. In our previous work (jlliadis et alJboiObh . the E x = 5288 keV (3~) l evel of 
Yokota et all |l98J) had been associated with the E x = 5217 keV level of lFvnbo. et al.l (J2000). The 
new measurements of lSetoodehnia et al.l (|2010l ) indicate that these are separate l evels, with energies 



of 52 26(3) and 5318(4) keV. Because of this development, it has been suggested (ISetoodehnia et al 
201lh that the doublet level (0+ and 4+) of the E x = 5168 keV of bardavan et al.1 (120071) is two sep- 



arate levels with different energies. The proposed assignments in ISetoodehnia et al.l (|201ll ) leave a 
3+ level missing; we therefore chose to associate the 0+ level with E x = 5127 keV, leaving E x = 5168 
keV with a 4+ assignment (each with C 2 S < 0.01). This allows the level at E x = 5318 keV to 
be designated as 3~ ( with C 2 S = . 36). O ther values for spins remain unchanged. Spectroscopic 



factors are taken from iMackh et al.l ( 19731 ). 



Along with a new Q-value (4398.72 ± 3.06 keV from IWang et al.l |2012| ). all resonance en- 
ergies were reevalua ted. Weighted averages of excitation energie s were calculated for th e res- 



onances at 289 keV (ISetoodehnia et al.ll2010l : iBardayan et al.l 120071 : ISetoodehnia et al.ll201ll ). 412 



80 



keV dSetoodehnia et all201ol.l2011 



1982 



2011), 819 keV (Setoodehnia et al 



Setoodehnia et al 



Yokota et al. 



2010 



Paddock 



. 734 keV (lYokota et al.lll982l:lKuh 



2010l:lFvnbo. et al 



2000 



1972 1. 993 keV (jFvnbo. et al. 



Paddock 



1972V 915 keV (lYokota et a 



maim et all 19 73; Setoodehnia et al 



200d: ISetoodehnia et al.ll2010l; 



19821 ; |Paddocklll972 ) . This was unnecessary for resonances at 769 keV ( Bardavan et al. 



20071 ) and 1443 keV ( Fvnbo. et al.ll2000l ). All penetration factors were recalcul ated with the new res - 
onance energies, and subsequent particle partial widths changed slightly from llliadis et al.1 ( 2010bl ). 
Gamma-ray partial widths are also unchanged, although which level each is associated with changed 
for consistency with the above spin assignment changes. For T > 1.5 GK, rates were extrapolated 
using Hauser-Feshbach calculations. 

The significance of the direct capture component was estimated through simulation. Above 
0.1 GK, the effect of increasing or decreasing S(0) by a factor of 5 is negligible; however, between 
0.01 and 0.1 GK, the rate changes linearly with this factor. 



29P(p,g)30S 






1 

15 



1.0078 

28.982 



0.5 

0.5 

0.0 

4398.72 

0.0 

1.25 

2 



Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 






1.0 
5000 





! Minimum energy for numerical integration (keV) 

! Number of random samples (>5000 for better statistics) 

! =0 for rate output at all temperatures; =NT for rate output at selected temperatures 

Non-resonant contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

7.3el -1.25e-2 0.0 0.4 1000.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 

Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 

Ecm DEcm wg Dwg J Gl DG1 LI G2 DG2 L2 G3 DG3 L3 Exf Int 

3 1.3e-5 5.3e-6 2 4.9e-3 2.5e-3 1 

2 3.3e-3 1.3e-3 2 4.2e-3 2.1e-3 1 

3 2.2e-l 8.8e-2 2 8.2e-3 4.1e-3 1 

3 l.leO 4.4e-l 3 9.4e-3 4.7e-3 1 



289.4 


3.1 





411.8 


3.1 





734.1 


3.1 





915.0 


4.8 















0.0 


1 











0.0 


1 











0.0 


1 











0.0 


1 
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992.5 3.5 2 6.4e0 2.6e0 2 0.01 5.0e-3 1 0.0 

Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 

Note: ...PT=<B> for g-rays [enter: "upper_limit 0.0"]; for each resonance: # upper limits < # open channels! 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf Int 

769.3 6.7 4 3.5e-4 0.0 4 0.0045 5.5e-3 2.8e-3 10 00 0.01 

819.1 3.1 1.8el 0.0 0.0045 7.7e-3 3.9e-3 10 00 0.01 

1443.3 5.0 4 1.0e-l 0.0 4 0.0045 2.7e-2 1.4e-2 10 00 0.01 

Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf 

!+- 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Resonance energies calculated from excitation energies (Endt 1990, Tab. I of Bardayan et al . 2007, 
Setoodehnia et al . 2010, and Setoodehnia et al., 2011) and Q-value, except for Er=489 keV, which is based on 
the IMME (see Iliadis et al. 2001). 

2. The spin-parity assignments are not unambiguous; they are based on experimental Jp restrictions (Endt 1990, 
Bardayan et al . 2007) and application of the IMME (Iliadis et al . 2001). 

3. Proton partial widths are computed using C2S of mirror states from (d,p) transfer experiment (Mackh et al. 
1973) . 

4. Gamma-ray partial widths are computed from measured lifetimes of 30Si mirror states, except for Er=990 keV 
for which Gg=0.012 eV is a rough estimate (application of RULs to the known g-ray branchings of the mirror 
state yields Gg<0 . 1 eV for this level). 
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9 P(P,Y) 30 S 



j 


\ 


v 


T9 - 0.03 
A-D = 0.1812 




T9 = 0.06 
A-D = 0.4395 



2e-26 4e-26 6e-26 8e-26 



CO 



CB 

o 




J 








T9 = 0.3 
A-D = 1.818 



0e-13 1.06-12 



0.0006 0.0010 0.0014 



-4™- 






T9 = 0.6 
A-D = 0.5699 




Reaction rate (cm 3 mol 1 s 1 ) 
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9 P(p,7) 30 S 





Temperature (GK) 
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Table 15:: Total thermonuclear reaction rates for 
38 Ar(p, 7 ) 39 K. 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal fj, 


lognormal a 


A-D 


0.010 


1.89 xlO" 46 


3.18xl0" 45 


2.20xl0" 44 


-1.028xl0+° 2 


2.28x10+°° 


3.11xl0+ 01 


0.011 


3.97xl0~ 45 


5.68xl0~ 44 


4.43 xlO" 43 


-9.975 xl0+ 01 


2.16x10+°° 


1.64xl0+ 01 


0.012 


1.53xl0~ 43 


9.00xl0~ 43 


5.92xl0" 42 


-9.676 xl0+ 01 


1.76x10+°° 


5.80x10+°° 


0.013 


7.87xl0~ 42 


3.61xl0~ 41 


1.53xl0" 4 ° 


-9.316xl0+ 01 


1.46x10+°° 


2.19x10+°° 


0.014 


4.53xl0~ 4 ° 


2.37xl0~ 39 


1.02xl0" 38 


-8.903 xl0+ 01 


1.55x10+°° 


3.75x10+°° 


0.015 


2.42xl0~ 38 


1.39xl0~ 37 


6.04xl0" 37 


-8.499 xl0+ 01 


1.61x10+°° 


9.28x10+°° 


0.016 


9.80xl0~ 37 


5.50xl0~ 36 


2.21 xlO" 35 


-8.135xl0+ 01 


1.58x10+°° 


1.55xl0+ 01 


0.018 


5.20xl0~ 34 


2.56xl0~ 33 


8.77xl0" 33 


-7.522xl0+ 01 


1.48x10+°° 


2.43xl0+ 01 


0.020 


7.76 xl0~ 32 


3.56xl0~ 31 


l.llxlO" 30 


-7.030 xl0+ 01 


1.40x10+°° 


3.31 xlO+ 01 


0.025 


6.33xl0~ 28 


2.53xl0~ 27 


7.08 xlO" 27 


-6.144xl0+ 01 


1.30x10+°° 


5.12xl0+ 01 


0.030 


2.40xl0~ 25 


9.39xl0~ 25 


2.56xl0" 24 


-5.552xl0+ 01 


1.27x10+°° 


5.48xl0+ 01 


0.040 


4.52xl0~ 22 


1.53xl0" 21 


4.28xl0~ 21 


-4.803 xl0+ 01 


1.11x10+°° 


1.77xl0+ 01 


0.050 


1.13xl0~ 19 


2.99xl0~ 19 


6.92xl0~ 19 


-4.272 xl0+ 01 


9.03xl0" 01 


1.13xl0+ 01 


0.060 


7.26xl0~ 18 


2.46xl0~ 17 


8.73xl0~ 17 


-3.825xl0+ 01 


1.19x10+°° 


5.51x10+°° 


0.070 


1.95xl0~ 16 


9.25xl0~ 16 


3.80 xl0~ 15 


-3.467xl0+ 01 


1.36x10+°° 


1.60xl0+ 01 


0.080 


2.97xl0~ 15 


1.54xl0 -14 


6.45xl0" 14 


-3.188xl0+ 01 


1.40x10+°° 


2.02xl0+ 01 


0.090 


2.77xl0~ 14 


1.38xl0 -13 


5.78xl0~ 13 


-2.966 xl0+ 01 


1.37x10+°° 


2.09xl0+ 01 


0.100 


1.92xl0~ 13 


8.13xl0~ 13 


3.33xl0~ 12 


-2.784xl0+ 01 


1.28x10+°° 


2.04xl0+ 01 


0.110 


1.19xl0~ 12 


3.88xl0~ 12 


1.41xl0 -11 


-2.624xl0+ 01 


1.14x10+°° 


1.21xl0+ 01 


0.120 


6.46xl0~ 12 


1.85xl0 -11 


5.18xl0~ n 


-2.473 xl0+ 01 


1.00x10+°° 


4.60x10+°° 


0.130 


2.87xl0" n 


8.26xl0~ n 


2.02xl0~ 10 


-2.329xl0+ 01 


9.56xl0- 01 


1.30xl0+ 01 


0.140 


l.lOxlO" 10 


3.31xl0~ 10 


8.42xl0" 10 


-2.190xl0+ 01 


1.00x10+°° 


1.07xl0+ 01 


0.150 


3.71xl0~ 10 


1.18xl0~° 9 


3.39 xlO" 09 


-2.060 xl0+ 01 


1.09x10+°° 


5.29x10+°° 


0.160 


1.16xl0~° 9 


3.92xl0~° 9 


1.27xl0"° 8 


-1.939xl0+ 01 


1.16x10+°° 


5.70x10+°° 


0.180 


9.13xl0~° 9 


3.41xl0~° 8 


1.20xl0~° 7 


-1.722xl0+ 01 


1.22x10+°° 


1.02xl0+ 01 


0.200 


6.11xl0-° 8 


2.15xl0~° 7 


7.38xl0~° 7 


-1.537xl0+ 01 


1.20x10+°° 


7.91x10+°° 


0.250 


2.58xl0~° 6 


8.16xl0~° 6 


2.17xl0~° 5 


-1.179xl0+ 01 


1.04x10+°° 


9.50x10+°° 


0.300 


3.67xl0~° 5 


LlOxlO" 04 


2.62xl0"° 4 


-9.208 xl0+ 00 


9.56xl0- 01 


1.67xl0+ 01 


0.350 


2.79xl0~° 4 


7.60xl0~° 4 


1.83xl0"° 3 


-7.232 xl0+ 00 


9.01xl0~ 01 


7.88x10+°° 


0.400 


1.49 xl0~° 3 


3.52xl0~° 3 


8.57xl0"° 3 


-5.636 xl0+ 00 


8.37xl0- 01 


4.05x10+°° 


0.450 


6.28xl0~° 3 


1.30xl0~° 2 


3.03xl0~° 2 


-4.299 xl0+ 00 


7.62xl0~ 01 


7.17x10+°° 


0.500 


2.18xl0~° 2 


4.04xl0~° 2 


8.68xl0~° 2 


-3.156x10+°° 


6.85xl0~ 01 


7.68x10+°° 


0.600 


1.59X10" 01 


2.69X10" 01 


4.67xl0~ 01 


-1.309x10+°° 


5.55xl0~ 01 


2.12x10+°° 


0.700 


7.13xl0~ 01 


1.13x10+°° 


1.75x10+°° 


1.060xl0~ 01 


4.67xl0~ 01 


2.52x10+°° 


0.800 


2.26xl0 +00 


3.47x10+°° 


5.04x10+°° 


1.218x10+°° 


4.09xl0~ 01 


6.04x10+°° 


0.900 


5.76x10+°° 


8.46x10+°° 


1.19xl0+ 01 


2.112x10+°° 


3.70xl0" 01 


7.80x10+°° 
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Table 15 - continued 



T (GK) 


Low rate 


Median rate 


High rate 


lognormal \x 


lognormal a 


A-D 


1.000 


1.24xl0+ 01 


1.76xl0+ 01 


2.41xl0+ 01 


2.845 xl0+ 00 


3.40xl0~ 01 


7.73x10+°° 


1.250 


5.10xl0+ 01 


6.86xl0 +01 


8.91 xl0+ 01 


4.211x10+°° 


2.84xl0" 01 


5.01x10+°° 


1.500 


1.38xl0 +02 


1.77xl0 +02 


2.22xl0 +02 


5.164x10+°° 


2.39xl0" 01 


2.78x10+°° 


1.750 


2.91 xl0+ 02 


3.59xl0+ 02 


4.36xl0+° 2 


5.876x10+°° 


2.03xl0~ 01 


1.47x10+°° 


2.000 


5.22xl0 +02 


6.24xl0+° 2 


7.37xl0 +02 


6.433x10+°° 


1.72X10" 01 


7.61 xl0~ 01 


2.500 


1.24xl0 +03 


1.42 xl0+ 03 


1.61 xl0 +03 


7.258x10+°° 


1.29xl0~ 01 


2.78xl0~ 01 


3.000 


2.30xl0 +03 


2.55xl0 +03 


2.81 xl0 +03 


7.842x10+°° 


1.01xl0~ 01 


1.89xl0~ 01 


3.500 


(3.59xl0 +03 ) 


(3.91xl0+° 3 ) 


(4.26xl0+° 3 ) 


(8.271x10+°°) 


(8.56xl0~° 2 ) 




4.000 


(4.98xl0 +03 ) 


(5.42xl0+° 3 ) 


(5.91 xl0 +03 ) 


(8.598x10+°°) 


(8.56xl0"° 2 ) 




5.000 


(8.18xl0+ 03 ) 


(8.92xl0+° 3 ) 


(9.71 xl0+ 03 ) 


(9.096x10+°°) 


(8.56xl0"° 2 ) 




6.000 


(1.17xl0+ 04 ) 


(1.28xl0+ 04 ) 


(1.39xl0+ 04 ) 


(9.455x10+°°) 


(8.56xl0~° 2 ) 




7.000 


(1.50xl0+ 04 ) 


(1.64xl0+ 04 ) 


(1.78xl0+ 04 ) 


(9.704x10+°°) 


(8.56xl0"° 2 ) 




8.000 


(1.70xl0+ 04 ) 


(1.85xl0+ 04 ) 


(2.02xl0+ 04 ) 


(9.827x10+°°) 


(8.56xl0"° 2 ) 




9.000 


(1.66xl0+ 04 ) 


(1.80xl0+ 04 ) 


(1.97xl0+° 4 ) 


(9.801x10+°°) 


(8.56xl0"° 2 ) 




10.000 


(1.61 xl0+ 04 ) 


(1.76xl0+ 04 ) 


(1.91 xl0+ 04 ) 


(9.774x10+°°) 


(8.56xl0" 02 ) 





C.8. 38 Ar(p, 7 ) 39 K 

Comments: In total, 118 resonances are taken into account. For 99 resonances i n the energy 
regi on of E r = 873 — 2338 keV, the measured energies and strengths are adopted from lEndtl (Il990l ) 



and iHanninenl (119841 ) . respectively. In the latter work, all the resonance strengths were normalized 
to E r = 1 358 k eV, <^7=1.30±0.25 eV. This value is in agreement with the earlier measurement of 
Maripuul (|l970l ). For the resonance at E r = 165 keV (J 77 ; T = 7/2"; 3/2), the strength is estimated 
us ing the proton spectr oscopic factor, C 2 S = 0.25, measured in the 38 Ar( 3 He,d) 39 K transfer study 



of iKnopfle et al.l (J1974I ). For 18 resonances between the lowest-lying directly measured resonance 
(E r = 873 keV) and the proton threshold, spins, parities and mean lifetimes are adopted from 



Endtl (jl990l ) or the Evaluated Nuclear Structure Data File (ENSDF). For these resonances, the 
dimensionless reduced proton widths are sampled randomly using a Porter-Thomas distribution 
with a cutoff at C 2 S = 1. The direct capture rate contribution is estimated by weighing the 
partial S-factors, obtai ned from a single-pa rticle potential model, with the experimental proton 
spectroscopic factors of IKnopfle et al.1 (119741 ) . At 3.4 GK, the rates need to be matched to Hauser- 
Feshbach results, for extrapolation of the laboratory rates to higher temperatures. In t he bottom 
panel of the comparison plot of the new Monte Carlo-based rates to the previous rates of lHanninen 
(| 19841 ). it can be seen that our results are much higher at temperatures between 0.5 GK and 2.5 GK. 
The reason is that the work of lHanninenl (]1984l ) disregarded any rate contributions from resonances 
below E r = 873 keV. 
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38Ar(p,g)39K 



1 

18 



1.0078 

37.962732 



0.5 

0.0 

0.0 

6381.43 

0.0 

1.25 

2 



Zproj 

Ztarget 

Zexitparticle (=0 when only 2 channels open) 

Aproj 

Atarget 

Aexitparticle (=0 when only 2 channels open) 

Jproj 

Jtarget 

Jexitparticle (=0 when only 2 channels open) 

projectile separation energy (keV) 

exit particle separation energy (=0 when only 2 channels open) 

Radius parameter R0 (fm) 

Gamma-ray channel number (=2 if ejectile is a g-ray; =3 otherwise) 



1.0 ! Minimum energy for numerical integration (keV) 

5000 ! Number of random samples (>5000 for better statistics) 

! =0 for rate output at all temperatures; =NT for rate output at selected temperatures 

Non-resonant contribution 

S(keVb) S'(b) S"(b/keV) fracErr Cutoff Energy (keV) 

741.0 -0.2598 1.208e-4 0.5 1600.0 

0.0 0.0 0.0 0.0 0.0 

Resonant Contribution 

Note: Gl = entrance channel, G2 = exit channel, G3 = spectator channel ! ! Ecm, Exf in (keV) ; wg, Gx in (eV) ! ! 
Note: if Er<0, theta~2=C2S*theta_sp~2 must be entered instead of entrance channel partial width 

DG1 LI G2 DG2 L2 G3 DG3 L3 
ie-13 3.8e-13 











































Ecm 


DEcm 


wg 


Dwg 


Jr 


Gl 


164.6 


2.0 


0.0 


0.0 


3. 


5 7 


873.2 


0.8 


0.085 


0.025 








897.3 


0.6 


0.165 


0.05 








905.2 


0.8 


0.02 


0.01 








955.0 


0.8 


0.03 


0.01 








1000.0 


0.6 


0.05 


0.02 








1057.6 


0.7 


0.175 


0.05 








1067.0 


0.8 


0.03 


0.01 








1080.0 


0.7 


0.37 


0.095 








1100.4 


0.7 


0.115 


0.03 








1154.3 


0.8 


0.5 


0.2 








1159.1 


0.8 


0.06 


0.02 








1170.4 


0.8 


0.04 


0.02 








1178.9 


0.8 


0.02 


0.01 








1220.9 


0.8 


0.06 


0.02 








1223.5 


0.8 


0.07 


0.025 








1251.8 


0.6 


0.27 


0.085 








1318.2 


0.8 


0.06 


0.02 








1332.9 


0.8 


0.145 


0.05 








1358.0 


0.6 


1.3 


0.25 








1374.3 


0.8 


0.135 


0.03 











G2 


DG2 


L2 


G3 


DG3 


3 


0. 


1 0.001 


1 



















































































































































































































































































































































































Exf 


Int 
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1384.6 


0.8 


0.05 


0.02 






































1391.5 


0.8 


0.145 


0.04 






































1403 . 2 


0.8 


0.08 


0.03 






































1415.9 


0.8 


0.435 


0.115 






































1420 . 5 


0.8 


0.215 


0.05 






































1423.3 


0.8 


0.1 


0.03 






































1439 . 2 


0.8 


0.205 


0.05 






































1465 . 3 


1.2 


0.3 


0.095 






































1486 . 6 


1.0 


0.38 


0.10 






































1576.3 


0.8 


0.39 


0.10 






































1597.1 


0.8 


0.25 


0.06 






































1600.1 


1.0 


0.07 


0.03 






































1601.8 


0.8 


0.5 


0.2 






































1604.8 


0.8 


0.06 


0.02 






































1610.8 


0.8 


0.465 


0.105 






































1613.9 


0.8 


0.07 


0.02 






































1616.9 


0.8 


0.06 


0.02 






































1649.5 


0.8 


0.25 


0.06 






































1652.4 


1.0 


0.105 


0.03 






































1656.9 


0.8 


0.65 


0.20 






































1698.1 


1.2 


0.175 


0.070 






































1699.9 


1.2 


0.260 


0.125 






































1705.5 


1.0 


0.03 


0.01 






































1711.7 


1.0 


0.27 


0.07 






































1717.6 


1.0 


0.175 


0.06 






































1726.2 


1.0 


0.135 


0.070 






































1736.4 


0.8 


0.435 


0.105 






































1747.1 


1.0 


0.09 


0.03 






































1756.8 


1.0 


0.095 


0.030 






































1788.7 


1.0 


0.155 


0.060 






































1803.4 


1.0 


0.280 


0.095 






































1807.6 


1.0 


0.260 


0.085 






































1809.7 


1.0 


0.280 


0.095 






































1816.7 


1.0 


0.6 


0.2 






































1821.6 


1.2 


0.10 


0.03 






































1871.2 


0.8 


1.35 


0.45 






































1880.8 


1.0 


0.70 


0.25 






































1889.4 


1.0 


0.125 


0.04 






































1898.1 


1.0 


0.185 


0.060 






































1903.6 


1.0 


0.095 


0.040 






































1912.6 


1.0 


0.235 


0.080 






































1923.4 


1.0 


0.465 


0.155 






































1933.1 


1.5 


0.145 


0.050 






































1944.4 


1.0 


0.07 


0.03 






































1959.3 


1.0 


0.310 


0.105 






































1965.4 


1.0 


0.185 


0.060 






































1998.2 


1.0 


0.195 


0.060 






































2004.8 


1.0 


0.55 


0.20 






































2013.7 


1.5 


0.260 


0.085 






































2032.4 


1.2 


0.6 


0.2 







































2045.4 


1.2 


0.185 


0.060 






































2048 . 1 


1.2 


0.475 


0.155 






































2083 . 9 


1.2 


0.235 


0.085 






































2094 . 6 


1.5 


0.270 


0.095 






































2102.7 


1.2 


0.50 


0.20 






































2127.8 


1.2 


0.135 


0.040 






































2132.2 


1.2 


0.185 


0.060 






































2142.8 


1.2 


0.195 


0.060 






































2149.0 


1.2 


0.55 


0.20 






































2166.1 


1.3 


0.330 


0.115 






































2177.9 


1.5 


0.105 


0.030 






































2185.7 


1.3 


0.135 


0.040 






































2202 . 2 


1.3 


0.80 


0.25 






































2210.8 


1.5 


0.205 


0.070 






































2216.6 


1.5 


0.105 


0.050 






































2243.4 


1.5 


0.310 


0.105 






































2245 . 5 


1.5 


0.85 


0.30 






































2256.9 


1.5 


0.29 


0.10 






































2273.5 


1.5 


0.205 


0.075 






































2285 . 3 


1.5 


0.125 


0.040 






































2290 . 8 


1.5 


0.185 


0.060 






































2292 . 8 


1.5 


0.155 


0.050 






































2301.9 


1.5 


0.415 


0.135 






































2306.4 


1.5 


0.415 


0.135 






































2311.7 


1.5 


0.260 


0.085 






































2322.4 


1.5 


0.310 


0.105 






































2327.6 


1.5 


0.6 


0.2 






































2332 . 7 


1.5 


0.310 


0.105 






































2338.2 


1.5 


0.290 


0.105 







































Upper Limits of Resonances 

Note: enter partial width upper limit by chosing non-zero value for PT, where PT=<theta~2> for particles and... 



Note : . 


. .PT= 


<B> 


for g 


-rays [en 


ter: " 


upper 


.limit 0.0 


"] ; for 


each resonance: 


# upper limits 


< # 


open 


channels ! 


Ecm 


DEcm 


Jr 


Gl 


DG1 LI 




PT 


G2 




DG2 L2 


PT 




G3 


DG3 L3 


PT 




Exf 


Int 


14.6 


2.0 




0.5 


3.6e-53 


0.0 







. 0045 


0.1 




0.001 


1 























28.6 


2.0 




1.5 


3.9e-37 


0.0 


2 




. 0045 


0.1 




0.001 


1 























76.6 


2.0 




0.5 


1.9e-17 


0.0 







. 0045 


0.1 




0.001 


1 























83.6 


2.0 




0.5 


2.9e-16 


0.0 







0.0045 


0.1 




0.001 


1 























119.6 


2.0 




1.5 


5.9e-14 


0.0 


2 




0.0045 


7.7e- 


■3 


3.8e-3 


1 




















1 


146.6 


2.0 




0.5 


8.6e-10 


0.0 







0.0045 


5.0e- 


■3 


2.5e-3 


1 




















1 


271.6 


2.0 




1.5 


1.7e-6 


0.0 


2 




0.0045 


1.9e- 


■2 


0.9e-2 


1 




















1 


304.6 


2.0 




0.5 


l.le-3 


0.0 







. 0045 


0.10 




0.05 


1 























358.6 


2.0 




1.5 


1.56-4 


0.0 


2 




. 0045 


0.19 




0.09 


1 




















1 


436.6 


3.0 




1.5 


2.3e-3 


0.0 


2 




. 0045 


0.14 




0.07 


1 




















1 


446.6 


2.0 




0.5 


0.26 


0.0 







. 0045 


0.10 




0.05 


1 























534.6 


2.0 




1.5 


0.032 


0.0 


2 




. 0045 


0.33 




0.17 


1 




















1 


561.6 


2.0 




1.5 


0.057 


0.0 


2 




0.0045 


0.10 




0.05 


1 























639.6 


2.0 




0.5 


18.4 


0.0 







0.0045 


0.10 




0.05 


1 























669.6 


2.0 




1.5 


10.3 


0.0 


1 




0.0045 


0.10 




0.05 


1 























740.6 


5.0 




1.5 


29.2 


0.0 


1 




. 0045 


0.10 




0.05 


1 
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788.6 2.0 0.5 152.0 0.0 0.0045 0.10 0.05 10 00 

818.6 2.0 1.5 78.0 0.0 1 0.0045 0.10 0.05 10 00 

Interference between Resonances [numerical integration only] 

Note: + for positive, - for negative interference; +- if interference sign is unknown 

Ecm DEcm Jr Gl DG1 LI PT G2 DG2 L2 PT G3 DG3 L3 PT Exf 

!+- 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

Reaction Rate and PDF at NT selected temperatures only 

Note: default values are used for reaction rate range if Min=Max=0 . 

T9 Min Max 

0.01 0.0 0.0 

0.1 0.0 0.0 

Comments : 

1. Excitation energies, spin-parities and lifetimes from Endt (1990) and ENSDF, unless noted otherwise. 

2. Resonance energies and strengths adopted from Haenninen (1984). In that work, all strengths were normalized 
to Erlab=1394 keV, wg=l .30+-0. 25 eV. This strength is in agreement with the earlier measurement of Maripuu 
(1970) . 

3. Direct capture S-f actor calcuted using measured spectroscopic factors from (3He,d) work of Knoepfle et al. 
(1974) . 

4. For Ercm=165 keV, Gp is calculated from measured C2S=0.25 of Knoepfle et al . (1974); assumption of Gg=0.1 
eV is inconsequential. 

5. For resonance upper limit calculation: (i) smalles possible orbital angular momentum ell assumed; (ii) high 
-spin states at 6434, 6475, 7092, 7142 keV have been disregarded since ell>4; (iii) for 6546 kev (T=3/2) 
level, C2S=0.245 adopted from (3He,d) study of Knoepfle et al. (1974); (iv) Gg=0 . 1 eV assumption for 
Ercm=15, 29, 77, 84 keV is inconsequential since Gg>>Gp; 

6. Gg=0.1 eV for Ercm=305, 447, 562, 640, 670, 741, 789, 819 keV is a crude guess. 

7. For Ercm=120, 147, 272, 359, 437 keV spin-parity is ambiguous; we assumed smallest possible spin for upper 
limit calculation. 
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